logo
    The Thermal Plumbing System of Stromboli Volcano, Aeolian Islands (Italy) Inferred From Electrical Conductivity and Induced Polarization Tomography
    3
    Citation
    69
    Reference
    10
    Related Paper
    Citation Trend
    Abstract:
    Abstract We performed the first 3D island‐scale tomography of the electrical conductivity of Stromboli volcano (Aeolian Islands, Italy) using 2D acquisition lines (37.2 km) and a total of 18,880 measurements and 2,402 unique electrode locations. This 3D data set was inverted using a Gauss‐Newton algorithm, parallel‐processing on an unstructured tetrahedral mesh containing 678,420 finite‐element nodes and 3,580,145 elements to account for the topography of the volcanic island. The tomogram exhibits a conductive body (10 −2 –1.0 S m −1 ) consistent with the location of CO 2 and temperature anomalies observed at the ground surface. It corresponds to the hydrothermal system with high electrical conductivity associated with alteration. In order to confirm this interpretation, a 2.5D large‐scale induced polarization tomography was performed crossing the volcano. The joint interpretation of the conductivity and normalized chargeability is done with a petrophysical model previously tested and verified at both shield‐ and strato‐volcanoes. This model implies that alteration (through the effect of the cation exchange capacity associated with clay minerals and zeolites) plays a strong role in both controlling the electrical conductivity and normalized chargeability at Stromboli volcano. A temperature tomogram, derived from the geoelectrical measurements, is consistent with surface temperature anomalies and the Very Long Period (VLP) seismicity related to the mild‐explosive activity. This survey displays at 600 m a.s.l. a lateral shift in the highest temperature location, also corresponding to the source of VLP seismicity. Structural boundaries have a major role in the hottest hydrothermal fluids rising below the active crater terrace of Stromboli volcano.
    Keywords:
    Electrical Resistivity Tomography
    Petrophysics
    Electrical Resistivity Tomography
    Electrical Impedance Tomography
    Sandbox (software development)
    Electrical capacitance tomography
    Abstract Determining soil hydraulic properties is complex, posing ongoing challenges in managing subsurface and agricultural practices. Electrical resistivity tomography (ERT) is an appealing geophysical method to monitor the subsurface due to its non‐invasive, easy‐to‐apply and cost‐effective nature. However, obtaining geoelectrical tomograms from raw measurements requires the inversion of an ill‐posed problem, which causes smoothing of the actual structure. Furthermore, the spatial resolution is determined from the distances in the electrode placement, thus inherently upscaling the obtained structure. This study explores the applicability of physics‐informed neural networks (PINNs) for upscaling permeability and petrophysical relations and monitoring water dynamics at heterogeneous soils using time‐lapse geoelectrical data. High‐resolution numerical simulations mimicking water infiltration were used as benchmarks. Synthetic ERT surveys with electrode spacing 10 times larger than the numerical model resolution were conducted to provide 2D electrical tomograms. The tomograms were fed to a PINNs system to obtain the permeability, petrophysical relations, and water content maps. An additional PINNs system incorporating water content measurements was trained to examine measurement sensitivity. Results have shown that the PINNs system could produce reliable results regarding the upscaled permeability and petrophysical relations fields. Water dynamics at the subsurface was accurately predicted with an average error of ∼3%. Adding water content measurements to PINNs training improved the system outcomes, mainly at the ERT low sensitivity zones. The PINNs system reduced water saturation errors by more than 30% compared to the common practice of directly translating the geoelectrical tomograms to water saturations using known, homogeneous petrophysical relations.
    Petrophysics
    Electrical Resistivity Tomography
    Water saturation
    Saturation (graph theory)
    Relative permeability
    Citations (0)
    In most cases of reservoir-induced seismicity, seismicity follows the impoundment, large lake-level changes, or filling at a later time above the highest water level achieved until then. We classify this as initial seismicity. This "initial seismicity" is ascribable to the coupled poroelastic response of the reservoir to initial filling or water level changes. It is characterized by an increase in seismicity above preimpoundment levels, large event(s), general stabilization and (usually) a lack of seismicity beneath the deepest part of the reservoir, widespread seismicity on the periphery, migrating outwards in one or more directions. With time, there is a decrease in both the number and magnitudes of earthquakes, with the seismicity returning to preimpoundment levels. However, after several years some reservoirs continue to be active; whereas, there is no seismicity at others. Preliminary results of two-dimensional (similar to those by Roeloffs, 1988) calculations suggest that, this "protracted seismicity" depends on the frequency and amplitude of lake-level changes, reservoir dimensions and hydromechanical properties of the substratum. Strength changes show delays with respect to lake-level changes. Longer period water level changes (~ 1 year) are more likely to cause deeper and larger earthquakes than short period water level changes. Earthquakes occur at reservoirs where the lake-level changes are comparable or a large fraction of the least depth of water. The seismicity is likely to be more widespread and deeper for a larger reservoir than for a smaller one. The induced seismicity is observed both beneath the deepest part of the reservoir and in the surrounding areas. The location of the seismicity is governed by the nature of faulting below and near the reservoir.
    Citations (43)