Imaging the magmatic system beneath the Krafla geothermal field, Iceland: A new 3-D electrical resistivity model from inversion of magnetotelluric data

2019 ◽  
Vol 220 (1) ◽  
pp. 541-567 ◽  
Author(s):  
Benjamin Lee ◽  
Martyn Unsworth ◽  
Knútur Árnason ◽  
Darcy Cordell

SUMMARY Krafla is an active volcanic field and a high-temperature geothermal system in northeast Iceland. As part of a program to produce more energy from higher temperature wells, the IDDP-1 well was drilled in 2009 to reach supercritical fluid conditions below the Krafla geothermal field. However, drilling ended prematurely when the well unexpectedly encountered rhyolite magma at a depth of 2.1 km. In this paper we re-examine the magnetotelluric (MT) data that were used to model the electrical resistivity structure at Krafla. We present a new 3-D resistivity model that differs from previous inversions due to (1) using the full impedance tensor data and (2) a finely discretized mesh with horizontal cell dimensions of 100 m by 100 m. We obtained similar resistivity models from using two different prior models: a uniform half-space, and a previously published 1-D resistivity model. Our model contains a near-surface resistive layer of unaltered basalt and a low resistivity layer of hydrothermal alteration (C1). A resistive region (R1) at 1 to 2 km depth corresponds to chlorite-epidote alteration minerals that are stable at temperatures of about 220 to 500 °C. A low resistivity feature (C2) coincides with the Hveragil fault system, a zone of increased permeability allowing interaction of aquifer fluids with magmatic fluids and gases. Our model contains a large, low resistivity zone (C3) below the northern half of the Krafla volcanic field that domes upward to a depth of about 1.6 km b.s.l. C3 is partially coincident with reported low S-wave velocity zones which could be due to partial melt or aqueous fluids. The low resistivity could also be attributed to dehydration and decomposition of chlorite and epidote that occurs above 500 °C. As opposed to previously published resistivity models, our resistivity model shows that IDDP-1 encountered rhyolite magma near the upper edge of C3, where it intersects C2. In order to assess the sensitivity of the MT data to melt at the bottom of IDDP-1, we added hypothetical magma bodies with resistivities of 0.1 to 30 Ωm to our resistivity model and compared the synthetic MT data to the original inversion response. We used two methods to compare the MT data fit: (1) the change in r.m.s. misfit and (2) an asymptotic p-value obtained from the Kolmogorov–Smirnov (K–S) statistical test on the two sets of data residuals. We determined that the MT data can only detect sills that are unrealistically large (2.25 km3) with very low resistivities (0.1 or 0.3 Ωm). Smaller magma bodies (0.125 and 1 km3) were not detected; thus the MT data are not sensitive to small rhyolite magma bodies near the bottom of IDDP-1. Our tests gave similar results when evaluating the changes in r.m.s. misfit and the K–S test p-values, but the K–S test is a more objective method than appraising a relative change in r.m.s. misfit. Our resistivity model and resolution tests are consistent with the idea of rhyolite melt forming by re-melting of hydrothermally altered basalt on the edges of a deeper magma body.

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Matthew J. Comeau ◽  
Michael Becken ◽  
Alexey V. Kuvshinov ◽  
Sodnomsambuu Demberel

AbstractCrustal architecture strongly influences the development and emplacement of mineral zones. In this study, we image the crustal structure beneath a metallogenic belt and its surroundings in the Bayankhongor area of central Mongolia. In this region, an ophiolite belt marks the location of an ancient suture zone, which is presently associated with a reactivated fault system. Nearby, metamorphic and volcanic belts host important mineralization zones and constitute a significant metallogenic belt that includes sources of copper and gold. However, the crustal structure of these features, and their relationships, are poorly studied. We analyze magnetotelluric data acquired across this region and generate three-dimensional electrical resistivity models of the crustal structure, which is found to be locally highly heterogeneous. Because the upper crust (< 25 km) is found to be generally highly resistive (> 1000 Ωm), low-resistivity (< 50 Ωm) features are conspicuous. Anomalous low-resistivity zones are congruent with the suture zone, and ophiolite belt, which is revealed to be a major crustal-scale feature. Furthermore, broadening low-resistivity zones located down-dip from the suture zone suggest that the narrow deformation zone observed at the surface transforms to a wide area in the deeper crust. Other low-resistivity anomalies are spatially associated with the surface expressions of known mineralization zones; thus, their links to deeper crustal structures are imaged. Considering the available evidence, we determine that, in both cases, the low resistivity can be explained by hydrothermal alteration along fossil fluid pathways. This illustrates the pivotal role that crustal fluids play in diverse geological processes, and highlights their inherent link in a unified system, which has implications for models of mineral genesis and emplacement. The results demonstrate that the crustal architecture—including the major crustal boundary—acts as a first‐order control on the location of the metallogenic belt.


2019 ◽  
Vol 1 (2) ◽  
pp. 41
Author(s):  
Triana Triana ◽  
Tony Yulianto ◽  
Udi Harmoko ◽  
Iqbal Takodama

Magnetotelluric data has been carried out at the "WS" geothermal field to analyze the resistivity model resulting from 2D inversion of magnetotelluric data in TE, TM and TE-TM modes. Base on the three models produced, the mode is determined to produce the most representative model to assist in the interpretation of the "WS" geothermal system. There is a step of modes separation, namely TE (Tranverse Electric) and TM (Transverse Magnetic) modes in processing MT data. Each mode produces a 2D model with different conductivity properties. The analysis results of the three modes explain that TE mode is dominated by low resistivity with a range of values of 10-35 Ωm and medium resistivity with a value range of 35-250 Ωm and a vertical resistivity contrast. The TM mode describes the high resistivity in the Southwest and the center of the track with a value of more than 470 sehinggam resulting in lateral resistivity contrast. While the TE-TM mode produces a model that is not much different from TM mode, only the distribution of the resistivity value is a combination with TE mode. This mode describes the distribution of resistivity both vertically and laterally. Based on the analysis of the three modes, it can be concluded that the TE-TM mode is the mode that produces the most representative model. Interpretation model shows that from the TE-TM mode we have a low resistivity distribution (10-35 Ωm) represent a cap rock zone, reservoir rock with a medium resistivity distribution (35-380 Ωm), resistive zone with a high resistivity distribution (more than 380 Ωm), and the existence of the three of faults structures ro be a controller system of the "WS" geothermal.


2021 ◽  
Author(s):  
◽  
Francesco Civilini

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction.  Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions.   We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells.  Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>


Geophysics ◽  
1971 ◽  
Vol 36 (2) ◽  
pp. 363-381 ◽  
Author(s):  
B. J. Mitchell ◽  
M. Landisman

Data from 20 magnetotelluric stations were used to determine a resistivity‐depth distribution for the crust in western Texas. We computed one‐dimensional models, which adequately explained the apparent resistivities obtained from the elements of a rotated impedance tensor. Descriptions of the method of data processing and model computation are given. The derived model includes an anisotropic, low resistivity layer at depths slightly greater than 20 km. The resistivity model is compared with a seismic velocity interpretation for eastern New Mexico, and the low resistivity layer is inferred to coincide with a zone of low seismic velocity. The base of the low velocity‐low resistivity region occurs at the transition from intermediate to more basic rocks in the crust. The low resistivity values are interpreted as resulting from interstitial fluid. This fluid may cause low seismic velocities and may form a weak region in which continental earthquakes occur. Interstitial water at depth in the crust may be produced by the intersection of the crustal temperature‐depth curve and equilibrium boundaries of hydrothermal reactions.


2021 ◽  
Author(s):  
Matthew Joseph Comeau ◽  
Michael Becken ◽  
James A. D. Connolly ◽  
Alexander Grayver ◽  
Alexey V. Kuvshinov ◽  
...  

&lt;p&gt;We investigate how a conceptual hydrodynamic model consisting of fluid localization and stagnation by thermally activated compaction can explain low-resistivity anomalies observed in the lower crust (&gt;20 km depth). Electrical resistivity models, derived from magnetotelluric data collected across the intracontinental Bulnay region, a subset of a larger regional array across central Mongolia, are generated. They reveal low-resistivity (3 - 30 &amp;#937;m) domains with a width of ~25 km and a vertical extent of &lt;10 km in the lower crust, with their tops ~5 km below the brittle-ductile transition zone. In 3-D these features appear as laterally extended (tube-like) structures, 300 km long, rather than disconnected ellipsoids. The features are oriented parallel to the adjacent Bulnay fault zone segments and perpendicular to the far-field compressive tectonic stress (i.e., northward motion from China and Tibet). These low-resistivity domains are consistent with the presence of saline metamorphic fluids. Deeper features imaged with the data include a large upper mantle conductor that we attribute to an asthenospheric upwelling, and thin lithosphere, related to intraplate surface uplift and volcanism, in agreement with recent geodynamic modelling of lithospheric removal in this region.&lt;/p&gt;&lt;p&gt;Based on the observed thermal structure of the crust, and assuming the mean stress at the brittle-ductile transition is twice the vertical load, the hydrodynamic model predicts that fluids would collect in zones &lt;9 km below the brittle-ductile transition zone, and the zones would have a vertical extent of ~9 km, both in agreement with the resistivity models across the Bulnay region. The hydrodynamic model also gives plausible values for the activation energy for viscous creep (270 - 360 kJ/mol), suggesting that the mechanism is dislocation creep.&lt;/p&gt;&lt;p&gt;From the electrical resistivity models, the lower crustal viscous compaction-length is constrained to be ~25 km - in this region. Within the conceptual model, this length-scale is entirely consistent with independent estimates for the specific hydraulic and rheological properties of this region. In fact, this can be used to independently constrain acceptable ranges for the lower crustal effective viscosity, which is found to be low (on the order of 10^18 Pas). Accordingly, the results indicate that low-salinity fluids (likely 1 - 0.01 wt% NaCl), and correspondingly low porosities (likely 5 - 0.1 vol%), are the most plausible. These key findings suggest partial melts are not favoured to explain the anomalies. Overall, the results of this contribution imply that it is tectonic and compaction processes that control lower crustal fluid flow, rather than lithological or structural heterogeneity.&lt;/p&gt;


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. WA95-WA104 ◽  
Author(s):  
Dale F. Rucker ◽  
Meng H. Loke ◽  
Marc T. Levitt ◽  
Gillian E. Noonan

An electrical-resistivity survey was completed at the T tank farm at the Hanford nuclear site in Washington State, U.S.A. The purpose of the survey was to define the lateral extent of waste plumes in the vadose zone in and around the tank farm. The T tank farm consists of single-shell tanks that historically have leaked and many liquid-waste-disposal facilities that provide a good target for resistivity mapping. Given that the site is highly industrialized with near-surface metallic infrastructure that potentially could mask any interpretable waste plume, it was necessary to use the many wells around the site as long electrodes. To accommodate the long electrodes and to simulate the effects of a linear conductor, the resistivity inversion code was modified to assign low-resistivity values to the well’s location. The forward model within the resistivity code was benchmarked for accuracy against an analytic solution, and the inverse model was tested for its ability to recreate images of a hypothetical target. The results of the tank-farm field survey showed large, low-resistivity targets beneath the disposal areas that coincided with the conceptual hydrogeologic models developed regarding the releases. Additionally, in areas of minimal infrastructure, the long-electrode method matched the lateral footprint of a 3D surface-resistivity survey with reasonable fidelity. Based on these results, the long-electrode resistivity method may provide a new strategy for environmental characterization at highly industrialized sites, provided a sufficient number and density of wells exist.


2021 ◽  
Author(s):  
◽  
Francesco Civilini

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction.  Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions.   We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells.  Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>


2017 ◽  
Vol 846 ◽  
pp. 012016
Author(s):  
Nazli Ismail ◽  
Muhammad Yanis ◽  
Syafrizal Idris ◽  
Faisal Abdullah ◽  
Bukhari Hanafiah

2020 ◽  
Vol 8 (2) ◽  
pp. 146
Author(s):  
Olawale OlakunleOsinowo ◽  
Ahmed Kehinde Usman ◽  
Ayotunde Allen Omitoogun

This study applied Electrical Resistivity Tomography geophysical investigation technique to evaluate the gold mineralization potential of Iperindo in Ilesha Schist Belt, southwestern Nigeria, where commercial exploitation capable of generating revenue and employment for the inhabitants has been challenged by lack / inadequate subsurface geological/geophysical information. The filtered and inverted electrical resistivity data acquired through five (5) 336 m long E – W trending profiles, established 10 m apart from each other, delineate isolated near surface but thick (> 30 m) low resistivity zones, especially at the eastern and western ends of the study area. Some of the delineated low resistivity zones (3 – 200 𝛀m) present vertical sharp edges, likely created by vertical faults that flank the zones on both sides. The low resistivity of these zones could be attributed to the occurrence of conductive material such as gold and associated base metals which probably exist in pegmatitic veins within the zones.   


Minerals ◽  
2020 ◽  
Vol 10 (6) ◽  
pp. 559 ◽  
Author(s):  
María Gabarrón ◽  
Pedro Martínez-Pagán ◽  
Marcos Martínez-Segura ◽  
María Bueso ◽  
Silvia Martínez-Martínez ◽  
...  

The legacy of the mining industry has left a large number of tailing ponds in the Cartagena–La Unión mining district exposed to water and wind erosion, which causes serious environmental and health problems and requires remediation. Before applying any remediation technique, an intensive sampling of the materials infilling the pond is required to determine the geochemistry of the pond, which will condition the remediation process. However, sampling the large number of tailing ponds that compose the district could be expensive. Thus, the main objective of this study is to evaluate the usefulness of electrical resistivity tomography (ERT) as a non-invasive tool to provide an image of spatial subsurface resistivity distribution and its relation to the physicochemical composition of near-surface mine wastes. To achieve this objective, three short ERT profiles were conducted, and 12 samples in each profile were collected at different depths for its geochemical characterization. Several non-linear regression models were fitted to predict physicochemical properties and metal concentrations from electrical resistivity measures. As a result, a high resistivity area was depicted in the ERT profiles G2 and G3, while the low resistivity ERT profile G1 was also obtained in accordance with the site’s surficial characteristics. Relationships among low resistivity values and high salinity, clay content, high metal concentrations, and mobility were established. Specifically, calibrated models were obtained for electrical conductivity, particle sizes of 0.02–50 µm and 50–2000 µm, total Zn and Cd concentration, and bioavailable Ni, Cd, and Fe. The ERT technique was shown to be a useful tool for the approximation of the location and distribution of the highest ranges of fine particle sizes, moisture, and, to a lesser extent, metal accumulation in the near-surface waste materials.


Sign in / Sign up

Export Citation Format

Share Document