scholarly journals Structural and Hydrothermal Inferences from a Magnetotelluric Survey across Mt. Ruapehu, New Zealand

2021 ◽  
Author(s):  
◽  
Stacey Maree Dravitzki

<p>This thesis explains the electrical conductivity structure of Mt. Ruapehu. To identify hydrothermal or volcanic components of the volcano, data from 25 magnetotelluric sites are analyzed. Data collected are first analyzed in the time domain prior to conversion into the frequency domain. Here, data are remote referenced, and the impedance tensors, tippers, apparent resistivity and phase values are calculated. These components are then analyzed to identify major features within the data. The new phase tensor ellipse method is applied to identify influential features and determine the dimensionality of data. This analysis indicates where it is appropriate to apply 1 or 2 dimensional inversion schemes. Dimensionality analysis led to 1-D modelling of the determinant impedance at each site; and limited 2-D profiles across the Tongariro Volcanic Centre boundaries. These models are used to create a simple 3-D structural model of the volcano that is then forward modelled. The results of the 3-D forward modelling indicate that the dominating features of the volcano's electrical structure have been identified in the previous models. Crater Lake is the only possible hydrothermal system on Mt. Ruapehu identified in this study. It is also very unlikely that any large coherent bodies of magma exist in the near surface. However, a second thin conductor laying somewhere between 10 and 30 km deep beneath the eastern flank may contain 13% melt and is the probable driving heat force beneath the volcano. The structure of Mt. Ruapehu can be split into seven layers. A resistive surface layer (100 ohm m) of young volcanic debris within the Tongariro Volcanic Centre that is up to 500 m thick near the crater.  A conductive layer (10 - 30 ohm m) of wet, fractured and altered volcanic debris underlaying the younger debris throughout the Tongariro Volcanic Centre.  A layer of Tertiary sediment under the Tongariro Volcanic Centre that extends to the south and west. This layer is electrically indistinguishable from the previous layer and extends to approximately sea level. A resistive layer (400 ohm m), and consistent with greywacke basement covers the entire field area.  A second conductive layer (20 ohm m) is identified under the eastern flank of the volcano somewhere between the depths of 10 and 30 km. This layer is likely to be the heat and magma source driving the volcanic activity.  A surrounding resistive layer extends beyond and below the second conductive layer mentioned above. This surrounding layer is electrically similar to the greywacke above.  A very high resistivity layer (7000 ohm m) is identified below 80 km deep, and may be associated with the land/sea boundary or subduction zone to the east.</p>

2021 ◽  
Author(s):  
◽  
Stacey Maree Dravitzki

<p>This thesis explains the electrical conductivity structure of Mt. Ruapehu. To identify hydrothermal or volcanic components of the volcano, data from 25 magnetotelluric sites are analyzed. Data collected are first analyzed in the time domain prior to conversion into the frequency domain. Here, data are remote referenced, and the impedance tensors, tippers, apparent resistivity and phase values are calculated. These components are then analyzed to identify major features within the data. The new phase tensor ellipse method is applied to identify influential features and determine the dimensionality of data. This analysis indicates where it is appropriate to apply 1 or 2 dimensional inversion schemes. Dimensionality analysis led to 1-D modelling of the determinant impedance at each site; and limited 2-D profiles across the Tongariro Volcanic Centre boundaries. These models are used to create a simple 3-D structural model of the volcano that is then forward modelled. The results of the 3-D forward modelling indicate that the dominating features of the volcano's electrical structure have been identified in the previous models. Crater Lake is the only possible hydrothermal system on Mt. Ruapehu identified in this study. It is also very unlikely that any large coherent bodies of magma exist in the near surface. However, a second thin conductor laying somewhere between 10 and 30 km deep beneath the eastern flank may contain 13% melt and is the probable driving heat force beneath the volcano. The structure of Mt. Ruapehu can be split into seven layers. A resistive surface layer (100 ohm m) of young volcanic debris within the Tongariro Volcanic Centre that is up to 500 m thick near the crater.  A conductive layer (10 - 30 ohm m) of wet, fractured and altered volcanic debris underlaying the younger debris throughout the Tongariro Volcanic Centre.  A layer of Tertiary sediment under the Tongariro Volcanic Centre that extends to the south and west. This layer is electrically indistinguishable from the previous layer and extends to approximately sea level. A resistive layer (400 ohm m), and consistent with greywacke basement covers the entire field area.  A second conductive layer (20 ohm m) is identified under the eastern flank of the volcano somewhere between the depths of 10 and 30 km. This layer is likely to be the heat and magma source driving the volcanic activity.  A surrounding resistive layer extends beyond and below the second conductive layer mentioned above. This surrounding layer is electrically similar to the greywacke above.  A very high resistivity layer (7000 ohm m) is identified below 80 km deep, and may be associated with the land/sea boundary or subduction zone to the east.</p>


2021 ◽  
Vol 869 (1) ◽  
pp. 012058
Author(s):  
M Jannah ◽  
M Muhammad ◽  
Marwan ◽  
Z Jalil

Abstract The 2D subsurface identification work of iron sands in Benteng Lubuk, Krueng Raya was successfully studied using the very low frequency method based on resistivity mode (VLF-R). This study aims to identify iron sand deposits in coastal areas using electromagnetic inversion. The inversion process shows a conductivity zone of iron sand area, where the resistive layer is strongly covered by a conductive layer above it. High resistivity values were found at 80-100 m stations. This layer has a resistivity value between 20000 – 40000 m and the conductivity value tend to be low. It is estimated that at this point there will only be manifestations of iron sand or sea water intrusion, due to the location of the track close to the coastline.


1978 ◽  
Vol 15 (10) ◽  
pp. 1539-1546 ◽  
Author(s):  
A. Koziar ◽  
D. W. Strangway

The audiofrequency magnetotelluric (AMT) method has been used to study permafrost thickness near Tuktoyaktuk, N.W.T. in the Mackenzie Delta. In the frequency range of 10 Hz–10 kHz the permafrost behaves as a simple resistive layer over a conductive layer. This simple two-layer model can be inverted by asymptotic models to give a unique value for the thickness of the highly resistive frozen layer. In areas of simple layering, these results correlate well with drilling. In areas of sharp lateral variations in resistivity, depths tend to be underestimated. Unlike other electrical methods, AMT is not hampered by the presence of a surface melt layer in the summer if the conductivity–thickness product of this 'active layer' is less than about 0.03 mho (0.03 S).


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. WA179-WA188 ◽  
Author(s):  
Alan Yusen Ley-Cooper ◽  
James Macnae ◽  
Andrea Viezzoli

Most airborne electromagnetic (AEM) data are processed using successive 1D approximations to produce stitched conductivity-depth sections. Because the current induced in the near surface by an AEM system preferentially circulates at some radial distance from a horizontal loop transmitter (sometimes called the footprint), the section plotted directly below a concentric transmitter-receiver system actually arises from currents induced in the vicinity rather than directly underneath. Detection of paleochannels as conduits for groundwater flow is a common geophysical exploration goal, where locally 2D approximations may be valid for an extinct riverbed or filled valley. Separate from effects of salinity, these paleochannels may be conductive if clay filled or resistive if sand filled and incised into a clay host. Because of the wide system footprint, using stitched 1D approximations or inversions may lead to misleading conductivity-depth images or sections. Near abrupt edges of an extensive conductive layer, the lateral falloff in AEM amplitudes tends to produce a drooping tail in a conductivity section, sometimes coupled with alocal peak where the AEM system is maximally coupled to currents constrained to flow near the conductor edge. Once the width of a conductive ribbon model is less than the system footprint, small amplitudes result, and the source is imaged too deeply in the stitched 1D section. On the other hand, a narrow resistive gap in a conductive layer is incorrectly imaged as a drooping region within the layered conductor; below, the image falsely contains a blocklike poor conductor extending to depth. Additionally, edge-effect responses often are imaged as deep conductors with an inverted horseshoe shape. Incorporating lateral constraints in 1D AEM inversion (LCI) software, designed to improve resolution of continuous layers, more accurately recovers the depth to extensive conductors. The LCI, however, as with any AEM modeling methodology based on 1D forward responses, has limitations in detecting and imaging in the presence of strong 3D lateral discontinuities of dimensions smaller than the annulus of resolution. The isotropic, horizontally slowly varying layered-earth assumption devalues and limits AEM’s 3D detection capabilities. The need for smart, fast algorithms that account for 3D varying electrical properties remains.


2021 ◽  
Vol 28 (2) ◽  
pp. 247-256
Author(s):  
Siming He ◽  
Jian Guan ◽  
Xiu Ji ◽  
Hang Xu ◽  
Yi Wang

Abstract. In spread spectrum induced polarization (SSIP) data processing, attenuation of background noise from the observed data is the essential step that improves the signal-to-noise ratio (SNR) of SSIP data. The time-domain spectral induced polarization based on pseudorandom sequence (TSIP) algorithm has been proposed to improve the SNR of these data. However, signal processing in background noise is still a challenging problem. We propose an enhanced correlation identification (ECI) algorithm to attenuate the background noise. In this algorithm, the cross-correlation matching method is helpful for the extraction of useful components of the raw SSIP data and suppression of background noise. Then the frequency-domain IP (FDIP) method is used for extracting the frequency response of the observation system. Experiments on both synthetic and real SSIP data show that the ECI algorithm will not only suppress the background noise but also better preserve the valid information of the raw SSIP data to display the actual location and shape of adjacent high-resistivity anomalies, which can improve subsequent steps in SSIP data processing and imaging.


2013 ◽  
Vol 135 (5) ◽  
Author(s):  
S. A. Sina ◽  
T. Farsadi ◽  
H. Haddadpour

In this study, the aeroelastic stability and response of an aircraft swept composite wing in subsonic compressible flow are investigated. The composite wing was modeled as an anisotropic thin-walled composite beam with the circumferentially asymmetric stiffness structural configuration to establish proper coupling between bending and torsion. Also, the structural model consists of a number of nonclassical effects, such as transverse shear, material anisotropy, warping inhibition, nonuniform torsional model, and rotary inertia. The finite state form of the unsteady aerodynamic loads have been modeled based on the indicial aerodynamic theory and strip theory in the subsonic compressible flow. Novel Mach dependent exponential approximations of the indicial aerodynamic functions have been developed. The extended Galerkin’s method was used to construct the mass, stiffness, and damping matrices of the nonconservative aeroelastic system. Eigen analysis of the system was performed to obtain the aeroelastic instability (divergence and flutter) boundaries. Also, solving the equations of motion in the time domain leads to the aeroelastic response of wing in different flight speeds. The obtained results are compared with the available results in the literature, which reveals an excellent agreement. The numerical results obtained in this article seek to clarify the effects of geometrical and material couplings and flight Mach number on the aeroelastic instability and response of composite wings in subsonic compressible flow.


Author(s):  
Sang Woo Kim ◽  
Svein Sævik ◽  
Jie Wu

Abstract This paper addresses the performance evaluation of an empirical time domain Vortex Induced Vibrations (VIV) model which has been developed for several years at NTNU. Unlike the frequency domain which is the existing VIV analysis method, the time domain model introduces new vortex shedding force terms to the well known Morison equation. The extra load terms are based on the relative velocity, a synchronization model and additional empirical coefficients that describe the hydrodynamic forces due to cross-flow (CF) and In-line (IL) vortex shedding. These hydrodynamic coefficients have been tuned to fit experimental data and by considering the results from the one of existing frequency domain VIV programs, VIVANA, which is widely used for industrial design. The feature of the time domain model is that it enables to include the structural non-linearity, such as variable tension, and time-varying flow. The robustness of the new model’s features has been validated by comparing the test results in previous researches. However, the riser used in experiments has a relatively small length/diameter (L/D) ratio. It implies that there is a need for more validation to make it applicable to real riser design. In this study, the time domain VIV model is applied to perform correlation studies against the Hanøytangen experiment data for the case of linear sheared current at a large L/D ratio. The main comparison has been made with respect to the maximum fatigue damage and dominating frequency for each test condition. The results show the time domain model showed reasonable accuracy with respect to the experimental and VIVANA. The discrepancy with regard to experiment results needs to be further studied with a non-linear structural model.


2017 ◽  
Vol 2017 ◽  
pp. 1-14 ◽  
Author(s):  
Djamel Boubaya

The Maghnia plain in western Algeria is filled by Plio-Quaternary and Miocene sediments that rest unconformably on a basement of Jurassic rocks. Electrical sounding (VES), magnetic data, well information, and hydrogeological data have been used to explore for groundwater potential in the Maghnia plain. The interpretation of Schlumberger sounding data was first calibrated with the lithology of available nearby wells. Four geoelectrical layers were identified within the study area. They are a thin near surface topsoil layer with variable resistivities, a moderate resistive aquifer (15–30 ohm-m), a resistive aquifer (40–70 ohm-m), and a conductive clay layer (1–10 ohm-m). Near Sidi Mbarek, the geoelectric section is reduced to three layers: a topsoil layer, a conductive layer corresponding to the Miocene marls, and a deep resistive layer that correlates with the Oxfordian sandstones. The interpretation of VES data and the enhancement techniques of magnetic data enabled the identification of a number of unmapped faults that occur near recharge zones close to adjacent mountains. This study enabled us to study the extension of the known Plio-Quaternary aquifer of the Maghnia plain and to explore the possible existence of a second deep groundwater aquifer in Oxfordian sandstones.


2009 ◽  
Vol 1165 ◽  
Author(s):  
Sebastian Lehmann ◽  
David Fuertes Marrón ◽  
José Manuel Merino Álvarez ◽  
Maximo Léon ◽  
Michael Tovar ◽  
...  

AbstractChalcopyrite-based devices show highest conversion efficiencies among present thin film architectures with values of 20% at laboratory scale. This outstanding performance has been achieved for quaternary Cu(Inx,Ga1-x)Se2 (x˜0.7) compound material. However, a strong correlation between the performance and the gallium content or, in other words, low versus high bandgap materials has been recognized. One critical issue in this discussion is the formation of a copper-depleted near-surface phase with 1:3:5 and 1:5:8 stoichiometries. In earlier reports, surface phases with corresponding compositions have been found on CuInSe2, CuGaSe2 and Cu(Inx,Ga1-x)Se2 thin films. These near-surface phases show a positive influence on the performance of cells based on low bandgap Cu(Inx,Ga1-x)Se2 material due to n-type inversion and band gap widening compared to bulk properties. A tendency towards a neutral or even a negative impact of the near-surface phase on wide band gap material (high gallium content) has recently been reported [1]. Nevertheless, the structural models of copper-poor chalcopyrite-related compounds have been controversially discussed in literature but a stannite-type structural model is most suitable as will be presented. In any case, the relation of the structural properties between chalcopyrite and 1:3:5 phases is crucial for the performance of related devices.In this contribution we will report about the structural analysis of the Cu(Inx,Ga1-x)3Se5 solid solution series by means of anomalous x-ray scattering using synchrotron radiation, powder and single crystal neutron diffraction. Contributions of the isoelectronic species Cu+ and Ga3+ could be separated by these experiments. Bulk samples synthesized from the elements and heat treated at 650°C after the main reaction step - the latter in order to allow equilibrium structure formation - were investigated. Structural data like lattice parameters, tetragonal distortion and cation distribution were obtained for the complete Cu(Inx,Ga1-x)3Se5 solid solution series. The stannite-type structural model was assigned to all members of the investigated 1:3:5s which will be strengthened by simulations. We observed that the tetragonal distortion vanishes for compositions close to a gallium content as used for highest efficiency Cu(Inx,Ga1-x)Se2 devices. However, the tetragonal distortion depends critically on the cation distribution which is in turn controlled by the thermal history of the sample, as we have recently reported for pure CuGaSe2 [1]. This means that we can plot a direct correlation for the misfit between chalcopyrite and 1:3:5 phases depending on the gallium content and the thermal treatment of the considered thin films. These results will widen the understanding of the chalcopyrite-based thin film photovoltaic devices.[1] S. Lehmann et al., Phys. Stat. Sol. A (in press)


Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. WB19-WB35 ◽  
Author(s):  
Cyril Schamper ◽  
Fayçal Rejiba ◽  
Roger Guérin

Electromagnetic induction (EMI) methods are widely used to determine the distribution of the electrical conductivity and are well adapted to the delimitation of aquifers and clayey layers because the electromagnetic field is strongly perturbed by conductive media. The multicomponent EMI device that was used allowed the three components of the secondary magnetic field (the radial [Formula: see text], the tangential [Formula: see text], and the vertical [Formula: see text]) to be measured at 10 frequencies ranging from 110 to 56 kHz in one single sounding with offsets ranging from 20 to 400 m. In a continuing endeavor to improve the reliability with which the thickness and conductivity are inverted, we focused our research on the use of components other than the vertical magnetic field Hz. Because a separate sensitivity analysis of [Formula: see text] and [Formula: see text] suggests that [Formula: see text] is more sensitive to variations in the thickness of a near-surface conductive layer, we developed an inversion tool able to make single-sounding and laterally constrained 1D interpretation of both components jointly, associated with an adapted random search algorithm for single-sounding processing for which almost no a priori information is available. Considering the complementarity of [Formula: see text] and [Formula: see text] components, inversion tests of clean and noisy synthetic data showed an improvement in the definition of the thickness of a near-surface conductive layer. This inversion code was applied to the karst site of the basin of Fontaine-Sous-Préaux, near Rouen (northwest of France). Comparison with an electrical resistivity tomography tends to confirm the reliability of the interpretation from the EMI data with the developed inversion tool.


Sign in / Sign up

Export Citation Format

Share Document