WALDIM: A code for the dimensionality analysis of magnetotelluric data using the Rotational Invariants of the Magnetotelluric Tensor

In the magnetotelluric (MT) method, the analysis of geoelectric dimensionality has acquired special importance in the last years, because numerical codes have made it possible to model and invert data using either one-dimensional (1D), two-dimensional (2D) or three-dimensional (3D) approaches. We present a FORTRAN code termed WALDIM to perform the dimensionality analysis of a set of MT data, according to the WAL rotational invariants criteria. These criteria are based on the possible annulment of the invariants of the MT tensor, which allow retrieval of as much information as possible from this tensor, without taking any a priori dimensionality assumption. When determining the dimensionality of real and therefore noisy data, two problems arise. The first is due to the data errors, which propagate into the invariants values, and therefore, to the determination of the dimensionality. The second is the fact that the invariants are rarely precisely zero, and the definition of a threshold is necessary. To solve these problems, WALDIM takes into account the data errors. Additionally, the dimensionality results can be grouped into frequency bands. Thus, we provide a software utility that allows providing a robust description of the dimensionality, and the parameters necessary for data correction prior to modeling. Given its completeness at analyzing the MT tensor for both individual and bands of frequencies, this code is meant to be a practical tool for MT data analysis

Capability of cross-hole electrical configurations for monitoring rapid plume migration experiments

Cross-hole electrical resistivity tomography is a useful tool in geotechnical, hydrogeological or fluid/gas plume migration studies. It allows better characterization of deep subsurface structures and monitoring of the involved processes. However, due to the large amount of possible four-electrode combinations between boreholes, the choice of the most efficient ones for rapid plume migration experiments (real-time monitoring), becomes a challenge. In this work, a numerical simulation to assess the capabilities and constraints of the most common cross-hole configurations for real-time monitoring is presented. Four-electrode configurations, sensitivity, dependence on the body location and amount of data were taken into account. The analysis of anomaly detection and the symmetry of the sensitivity pattern of cross-hole configurations allowed significant reduction of the amount of data and maintaining the maximum potential resolution of each configuration for real-time monitoring. The obtained results also highlighted the benefit of using the cross-hole AB-MN configuration (with both current - or potential - electrodes located in the same borehole) combined with other configurations with complementary sensitivity pattern

Land CSEM Simulations and Experimental Test Using Metallic Casing in a Geothermal Exploration Context: Vallès Basin (NE Spain) Case Study

Controlled-source electromagnetic (CSEM) measurements are complementary data for magnetotelluric (MT) characterization although its methodology on land is not sufficiently developed and tested as in marine environments. Acquiring expertise in CSEM is crucial for surveys in places where MT cannot be performed due to high levels of cultural noise. To acquire that expertise, we perform CSEM experiments in the Vallès fault [Northeast (NE), Spain], where MT results have been satisfactory and allow us to verify the CSEM results. The Vallès basin is relevant for potential heat generation because of the presence of several geothermal anomalies and its nearby location in urban areas. In this article, we present the experimental setup for that region, a 2-D joint MT+CSEM inverse model, several 3-D CSEM simulations in the presence of metallic casing, and its comparison with real data measurements. We employ a parallel and high-order vector finite element algorithm to discretize the governing equations. By using an adapted meshing strategy, different scenarios are simulated to study the influence of the source position/direction and the conductivity model in a metallic casing presence. An excellent agreement between the simulated data and analytical/real field data demonstrates the feasibility of study metallic structures in realistic configurations. Our numerical results confirm that metallic casing strongly influences electromagnetic (EM) responses, making surface measurements more sensitive to resistivity variations near the metallic structure. It could be beneficial getting higher signal-to-noise ratios and sensitivity to deep targets. However, such a casing effect depends on the input model (e.g., conductivity contrasts, frequency, and geometry)

Dimensionality imprint of electrical anisotropy in magnetotelluric responses.

Dimensionality analysis of magnetotelluric data is a common procedure for inferring the main properties of the geoelectric structures of the subsurface such as the strike direction or the presence of superficial distorting bodies, and enables the most appropriate modeling approach (1D, 2D or 3D) to be determined. Most of the methods currently used assume that the electrical conductivity of individual parts of a structure is isotropic, although some traces of anisotropy in data responses can be recognized. In this paper we investigate the imprints of anisotropic media responses in dimensionality analysis using rotational invariants of the magnetotelluric tensor. We show results for responses generated from 2D synthetic anisotropic models and for field data that have been interpreted as showing the effects of electrical anisotropy in parts of the subsurface structure. As a result of this study we extend the WAL dimensionality criteria to include extra conditions that allow anisotropic media to be distinguished from 2D isotropic ones. The new conditions require the analysis of the strike directions obtained and take into account the overall behavior of different sites in a survey

Improving Bahr's invariant parameters using the WAL approach.

In the magnetotelluric technique, several methods exist to perform dimensionality analysis of the measured data using rotational invariants of the impedance tensor. Among these methods there is some dilemma on the different criteria established, which sometimes lead to non-equivalent interpretation

On the detectability of Teide volcano magma chambers (Tenerife, Canary Islands) with magnetotelluric data

Tenerife has been the subject of numerous studies covering a wide range of fields. Many studies have been focused on characterising the magmatic plumbing system. Even so, a controversy still exists regarding the location and size of the current magma chambers. Several magnetotelluric (MT) surveys have been carried out in the island, but no conductivity anomalies associated with the chambers have been detected. We report the results of a set of tests conducted against the 3-D resistivity model of the island, to determine the characteristics of the detectable chambers with the MT data. The most remarkable results indicate that the MT dataset is incompatible with a large-scale mafic reservoir located at shallower depths than 8 km b.s.l. However, shallower phonolitic chambers smaller than 3 x 3 x 1 km(3) could be undetected by the existing MT sites and new data should be acquired to confirm or not their existence. This new information is essential in volcanic islands like Tenerife, since many volcanic hazards are related to the size and depth of the sources of magma. Additionally, a joint interpretation of the obtained results together with other information is summarised in a hypothetical model, allowing us to better understand the internal structure of the island

El potencial espontáneo como herramienta para la caracterización de la afección por hidrocarburos

El objeto del presente documento es poner en conocimiento de la sociedad la disponibilidad de una técnica de prospección geofísica no destructiva que permite detectar y localizar la afección por hidrocarburos en el subsuelo de instalaciones industriales.La memoria recoge la descripción del fenómeno físico y los antecedentes de dicha técnica, citándose también novedades respecto a su aplicación en instalaciones industriales.El documento puede ser de utilidad para los agentes que intervienen en actividades potencialmente contaminantes del suelo, así como de cualquier entidad que pueda acometer trabajos de caracterización de suelos con técnicas que complementen la ejecución de estas tareas

A layer stripping approach for monitoring resistivity variations using surface magnetotelluric responses

The resolution of surface-acquired magnetotelluric data is typically not sufficiently high enough in monitoring surveys to detect and quantify small resistivity variations produced within an anomalous structure at a given depth within the subsurface. To address this deficiency we present an approach, called "layer stripping", based on the analytical solution of the one-dimensional magnetotelluric problem to enhance the sensitivity of surface magnetotelluric responses to such subtle subsurface temporal variations in resistivity within e.g. reservoirs. Given a well-known geoelectrical baseline model of a reservoir site, the layer stripping approach aims to remove the effect of the upper, unchanging structures in order to simulate the time-varying magnetotelluric responses at depth. This methodology is suggested for monitoring all kinds of reservoirs, e.g. hydrocarbons, gas, geothermal, compress air storage, etc., but here we focus on CO2 geological storage. We study one-dimensional and three-dimensional resistivity variations in the reservoir layer and the feasibility of the method is appraised by evaluating the error of the approach and defining different detectability parameters. The geoelectrical baseline model of the Hontomín site (Spain) for CO2 geological storage in a deep saline aquifer is taken as our exemplar for studying the validity of the 1D assumption in a real scenario. We conclude that layer stripping could help detect resistivity variations and locate them in the space, showing potential to also sense unforeseen resistivity variations at all depths. The proposed approach constitutes an innovative contribution to take greater advantage of surface magnetotelluric data and to use the method as a cost-effective permanent monitoring technique in suitable geoelectrical scenarios

Improving the modeling of geomagnetically induced currents in Spain

Vulnerability assessments of the risk posed by geomagnetically induced currents (GICs) to power transmission grids benefit from accurate knowledge of the geomagnetic field variations at each node of the grid, the Earth's geoelectrical structures beneath them, and the topology and relative resistances of the grid elements in the precise instant of a storm. The results of previous analyses on the threat posed by GICs to the Spanish 400 kV grid are improved in this study by resorting to different strategies to progress in the three aspects identified above. First, although at midlatitude regions the source fields are rather uniform, we have investigated the effect of their spatial changes by interpolating the field from the records of several close observatories with different techniques. Second, we have performed a magnetotelluric (MT) sounding in the vicinity of one of the transformers where GICs are measured to determine the geoelectrical structure of the Earth, and we have identified the importance of estimating the MT impedance tensor when predicting GIC, especially where the effect of lateral heterogeneities is important. Finally, a sensitivity analysis to network changes has allowed us to assess the reliability of both the information about the network topology and resistances, and the assumptions made when all the details or the network status are not available. In our case, the most essential issue to improve the coincidence between model predictions and actual observations came from the use of realistic geoelectric information involving local MT measurements

Geodynamic implications for the formation of the Betic-Rif orogen from magnetotelluric studies

Magnetotelluric data from the central Betics mountains (Spain) have been used to determine the electrical resistivity of the crust after a three‐dimensional (3D) interpretive approach. At shallow levels (3 km), the most striking and well‐resolved feature of the model is an upper‐middle crust conductive body, located at the core of the Internal Betics antiform. This approximately 14‐km‐thick body is interpreted as basic or ultrabasic rocks containing a conducting mineral phase. Its structural location above the sole thrust of the Betic orogen and beneath the Nevado‐Filábride complex confirms the presence of a major suture zone between this complex and the autochthonous Iberian plate. This suture may correspond to an ancient oceanic or transitional domain developed between Iberia and the Alboran Domain during the opening of the Tethys Ocean, partially subducted and closed during the development of the Betic orogen. The possible geodynamic scenarios for the Betics have been reconsidered, taking into account this new constraint