Search

The Caribbean and South American tectonic plates bound the north-eastwards expulsion of the North Andean Block in western Venezuela. This complex geodynamic setting resulted in the formation of major strike-slip fault systems and sizeable mountain chains. The 100-km-wide Mérida Andes extend from the Colombian/Venezuelan border to the Caribbean coast. To the north and south, the Mérida Andes are bound by hydrocarbon-rich sedimentary basins. Knowledge of lithospheric structures, related to the formation of the Mérida Andes, is limited though, due to a lack of deep geophysical data. In this study, we present results of the first broad-band magnetotelluric profile crossing the Mérida Andes and the Maracaibo and Barinas–Apure foreland basins on a length of 240 km. Geoelectrical strike and dimensionality analysis are consistent with 1-D or 2-D subsurface structures for the sedimentary basins but also indicate a strong 3-D setting for the Mérida Andes. Using a combination of 2-D and 3-D modelling we systematically examined the influence of 3-D structures on 2-D inversions. Synthetic data sets derived from 3-D modelling allow identification and quantification of spurious off-profile features as well as smoothing artefact due to limited areal station coverage of data collected along a profile. The 2-D inversion models show electrically conductive basins with depths of 2–5 km for the Barinas-Apure and 2–7 km for the Maracaibo basins. A number of resistive bodies within the Maracaibo basin could be related to active deformation causing juxtaposition of older geological formations and younger basin sediments. The most important fault systems of the area, the Boconó and Valera Faults, cross-cut the Mérida Andes in NE–SW direction along its strike on a length 400 km and N–S direction at its centre on a length 60 km, respectively. Both faults are associated with subvertical zones of high electrical conductivity and sensitivity tests suggest that they reach depths of up to 12 km. A sizeable conductor at 50 km depth, which appears consistently in the 2-D sections, could be identified as an inversion artefact caused by a conductor east of the profile. We speculate the high conductivity associated with the off-profile conductor may be related to the detachment of the Trujillo Block. Our results partially support the 'floating orogen hypothesis' developed to explain the geodynamic evolution of western Venezuela and they highlight the relevance of the Trujillo Block in this process. [ABSTRACT FROM AUTHOR]

We investigate seismoelectric (SE) signals accompanying seismic waves radiated from earthquake sources. SE signals are mostly generated from compressional portions of seismic waves by electrokinetic coupling. They contain coseismic electric fields travelling with seismic wave velocity and interface response (IR) waves, which originate at layer interfaces and travel with electromagnetic wave speed. IR wave amplitudes are sensitive to contrasts in poroelastic and electric rock parameters. We introduce SE spectral ratios (SESRs) as a tool to evaluate the influence of IRs on the overall SE signal independently of the earthquake source–time function. Based on data from Northern Chile we show that SESRs show a site specific frequency dependence with a trend of decreasing amplitudes towards increasing frequency. Modelling results reveal that the specific frequency dependence of the SESRs is caused by IRs excited at depths of some hundred metres underneath the recording stations. We analyse the SESR sensitivity towards porosity, permeability, fluid salinity and the depth of the interfaces. We verify that observed SESRs can be reproduced through forward modelling and linearized inversion based on realistic subsurface parameters. [ABSTRACT FROM AUTHOR]

We present the crustal resistivity structure of the Pamir and Southern Tian Shan orogenic belts at the northwestern promontory of the India–Asia collision zone. The magnetotelluric (MT) data were recorded along a roughly north–south trending, 350 km long corridor from the Pamir Plateau in southern Tajikistan across the Pamir frontal ranges, the Alai Valley and the southwestern Tian Shan to Osh in the Kyrgyz part of the Fergana Basin. In total, we measured at 178 sites, whereof 26 combine broad band and long period recordings. One of the most intriguing features of the 2-D and 3-D inversion results is a laterally extended zone of high electrical conductivity below the Pamir Plateau, with resistivities below 1 Ωm, starting at a depth of ∼10–15 km. The high conductivity can be explained with the presence of partially molten rocks at middle to lower crustal levels, possibly related to ongoing migmatization and/or middle/lower crustal flow underneath the Southern Pamir. This interpretation is consistent with a low velocity zone found from local earthquake tomography, relatively high vp/vs ratios, elevated surface heat flow, and thermomechanical modelling suggesting that melting temperatures are reached in the felsic middle crust. In the upper crust of the Pamir and Tian Shan, the Palaeozoic–Mesozoic suture zones appear as electrically conductive, whereas the compact metamorphic rocks of the Muskol-Shatput Dome of the Central Pamir are highly resistive. The intra-montane basin of the Alai Valley—sandwiched between the Pamir and Tian Shan—exhibits a generally conductive upper crust that bifurcates into two conductors at depth. One of them connects to the active Main Pamir Thrust, which is absorbing most of today's convergence between the Pamir and the Tian Shan. Several deeper zones of high conductivity in the middle and lower crust of Central and Northern Pamir likely record fluid release due to metamorphism associated with active continental subduction/delamination. [ABSTRACT FROM PUBLISHER]

For 3-D inversion of controlled-source electromagnetic (CSEM) data, increasing availability of high-performance computers enables us to apply inversion techniques that are theoretically favourable, yet have previously been considered to be computationally too demanding. We present a newly developed parallel distributed 3-D inversion algorithm for interpreting CSEM data in the frequency domain. Our scheme is based on a direct forward solver and uses Gauss–Newton minimization with explicit formation of the Jacobian. This combination is advantageous, because Gauss–Newton minimization converges rapidly, limiting the number of expensive forward modelling cycles. Explicit calculation of the Jacobian allows us to (i) precondition the Gauss–Newton system, which further accelerates convergence, (ii) determine suitable regularization parameters by comparing matrix norms of data- and model-dependent terms in the objective function and (iii) thoroughly analyse data sensitivities and interdependencies. We show that explicit Jacobian formation in combination with direct solvers is likely to require less memory than combinations of direct solvers and implicit Jacobian usage for many moderate-scale CSEM surveys. We demonstrate the excellent convergence properties of the new inversion scheme for several synthetic models. We compare model updates determined by solving either a system of normal equations or, alternatively, a linear least-squares system. We assess the behaviour of three different stabilizing functionals in the framework of our inversion scheme, and demonstrate that implicit regularization resulting from incomplete iterative solution of the model update equations helps stabilize the inversion. We show inversions of models with up to two million unknowns in the forward solution, which clearly demonstrates applicability of our approach to real-world problems. [ABSTRACT FROM AUTHOR]

Electrical conductivity is a key parameter for the exploration and characterization of geothermal reservoirs as hot mineralized formation water of active geothermal areas usually exhibits significantly higher conductivity than the surrounding host rock. Here we present results of a magnetotelluric (MT) exploration experiment carried out in the vicinity of the Groß Schönebeck geothermal test site in Northern Germany, where a doublet system of two 4.3-km deep boreholes was drilled to establish an in situ laboratory to investigate the potential for geothermal energy production. Classical 2-D smooth inversion of the MT data, recorded along two profiles, reveals a shallow conductive structure in good agreement with information from regional geology and seismic images. However, at the northernmost part of the profiles, the conductivity models reveal deep-reaching conductive structures, which appear uncorrelated with existing (geophysical or geological) data. Incorporating information from seismics as independent constraints for MT inversions allows us to examine the model space rigorously but target oriented. Employing so-called tear-zone inversions we can effectively derive an alternative class of models, which are consistent with the MT observations but also with the other data sets. We speculate that the zones of high conductivity imaged at reservoir depth are related areas of reduced thickness of the overlaying evapourite layer. The enhanced conductivity can be explained by a higher fracture density in anhydritic layers and/or generally lower resistivity of the pore fluid. [ABSTRACT FROM AUTHOR]

Magnetotelluric (MT) data from 66 sites along a 45-km-long profile across the San Andreas Fault (SAF) were inverted to obtain the 2-D electrical resistivity structure of the crust near the San Andreas Fault Observatory at Depth (SAFOD). The most intriguing feature of the resistivity model is a steeply dipping upper crustal high-conductivity zone flanking the seismically defined SAF to the NE, that widens into the lower crust and appears to be connected to a broad conductivity anomaly in the upper mantle. Hypothesis tests of the inversion model suggest that upper and lower crustal and upper-mantle anomalies may be interconnected. We speculate that the high conductivities are caused by fluids and may represent a deep-rooted channel for crustal and/or mantle fluid ascent. Based on the chemical analysis of well waters, it was previously suggested that fluids can enter the brittle regime of the SAF system from the lower crust and mantle. At high pressures, these fluids can contribute to fault-weakening at seismogenic depths. These geochemical studies predicted the existence of a deep fluid source and a permeable pathway through the crust. Our resistivity model images a conductive pathway, which penetrates the entire crust, in agreement with the geochemical interpretation. However, the resistivity model also shows that the upper crustal branch of the high-conductivity zone is located NE of the seismically defined SAF, suggesting that the SAF does not itself act as a major fluid pathway. This interpretation is supported by both, the location of the upper crustal high-conductivity zone and recent studies within the SAFOD main hole, which indicate that pore pressures within the core of the SAF zone are not anomalously high, that mantle-derived fluids are minor constituents to the fault-zone fluid composition and that both the volume of mantle fluids and the fluid pressure increase to the NE of the SAF. We further infer from the MT model that the resistive Salinian block basement to the SW of the SAFOD represents an isolated body, being wide and reaching to depths , in agreement with aeromagnetic data. This body is separated from a massive block of Salinian crust farther to the SW. The NE terminus of resistive Salinian crust has a spatial relationship with a near-vertical zone of increased seismic reflectivity SW of the SAF and likely represents a deep-reaching fault zone. [ABSTRACT FROM AUTHOR]

Magnetotelluric (MT) responses in complex, 3-D terrains are in general characterized by (i) elliptical polarization states of horizontal electric and magnetic fields; (ii) the non-orthogonality of electric and magnetic fields and (iii) a coupling of the anomalous tangential-electric (TE) and tangential-magnetic (TM) modes, giving rise to a mode-mixed anomalous electric field at the surface. These 3-D effects are propagated into the MT impedance tensor, which is derived from horizontal electric and magnetic fields, recorded at the earth's surface. The impedance tensor is in general fully occupied, and each of its elements is a mode-mixed quantity. To study 3-D effects of MT (impedance) data, the TE and TM mode contributions must be separated. This becomes possible with the inclusion of the single-mode vertical magnetic transfer function (the ratio of vertical to horizontal magnetic fields). Then, the individual modes can be resolved without prior knowledge of the underlying 3-D conductivity structure. For this purpose, we consider (i) the spatial relations between electromagnetic field components recorded in an array of sites (Faraday's law) and (ii) that the magnetic TE mode and electric TM mode fields are potential fields within the insulating air half-space above the earth's surface. Based on these two dependencies, it is possible to reconstruct the entire electromagnetic field from (measured) mixed-mode impedances and vertical magnetic transfer functions and to separate it into TE and TM modes, and into normal and anomalous parts. Hereby, we cannot only study the contribution of the two modes on the observed impedance tensor but also quantify the influence of 3-D effects at each location and frequency of a particular data set. Results of a modelling study suggest, that (i) none of the elements of a 3-D impedances tensor can be regarded as favourable for a 2-D interpretation (only 3-D models can explain 3-D data), (ii) a heterogeneous crust can strongly obscure identification of responses originating from the lower crust or upper mantle even at very long periods and (iii) the TE mode magnetic transfer functions are most important for sensing deep anomalies. [ABSTRACT FROM AUTHOR]

Magnetotelluric and seismic methods provide complementary information about the resistivity and velocity structure of the subsurface on similar scales and resolutions. No global relation, however, exists between these parameters, and correlations are often valid for only a limited target area. Independently derived inverse models from these methods can be combined using a classification approach to map geologic structure. The method employed is based solely on the statistical correlation of physical properties in a joint parameter space and is independent of theoretical or empirical relations linking electrical and seismic parameters. Regions of high correlation (classes) between resistivity and velocity can in turn be mapped back and re-examined in depth section. The spatial distribution of these classes, and the boundaries between them, provide structural information not evident in the individual models. This method is applied to a 10 km long profile crossing the Dead Sea Transform in Jordan. Several prominent classes are identified with specific lithologies in accordance with local geology. An abrupt change in lithology across the fault, together with vertical uplift of the basement suggest the fault is sub-vertical within the upper crust. [ABSTRACT FROM AUTHOR]

Magnetotelluric (MT) response function estimates can be severely disturbed by the effects of cultural noise. Methods to isolate and remove these disturbances are typically based on time-series editing, robust statistics, remote reference processing, or some combination of the above. Robust remote reference processing can improve the data quality at a local site, but only if synchronous recordings of at least one additional site are available and if electromagnetic noise between these sites is uncorrelated. If these prerequisites are not met, we suggest an alternative approach for noise removal, based on a combination of frequency domain editing with subsequent single site robust processing. The data pre-selection relies on a thorough visual inspection of a variety of statistical parameters such as spectral power densities, coherences, the distribution of response functions and their errors, etc. Extreme outliers and particularly noisy data segments are excluded from further data processing by setting threshold values for individual parameters. Examples from Namibia and Jordan illustrate that this scheme can improve data quality significantly. However, the examples also suggest that it is not possible to establish generally valid rules for selection as they depend strongly on the local noise conditions. High coherence, for example, can indicate a good signal-to-noise ratio or strongly correlated noise. However, we found that strong polarization of the magnetic field channels and the distribution of response function errors are two important parameters for noise detection. [ABSTRACT FROM AUTHOR]