Phreatic and phreatomagmatic eruptions represent some of the greatest hazards occurring on volcanoes. They result from complex interactions at a depth between rock, water, and magmatic fluids. Understanding and assessing such processes remain a challenging task, notably because a large‐scale characterization of volcanic edifices is often lacking. Here we focused on Miyakejima Island, an inhabited 8‐km‐wide stratovolcano with regular phreatomagmatic activity. We imaged its plumbing system through a combination of four geophysical techniques: magnetotellurics, seismicity, self‐potential, and thermal image. We thus propose the first comprehensive interpretation of the volcanic island in terms of rock properties, temperature, fluid content, and fluid flow. We identify a shallow aquifer lying above a clay cap (<1 km depth) and reveal its relation with magmatic‐tectonic features and past eruptive activity. At greater depths (2–4.5 km), we infer a seismogenic resistive region interpreted as a magmatic gas‐rich reservoir (≥370°C). From this reservoir, gases rise through a fractured conduit before being released in the fumarolic area at ∼180°C. During their ascent, these hot fluids cross a ∼1.2‐km‐long liquid‐dominated zone causing local steam explosions. Such magmatic‐hydrothermal interaction elucidates (i) the origin of the long‐period seismic events and (ii) the mixing mechanism between magmatic and hydrothermal fluids, which was previously observed in the geochemical signature of fumaroles. Our results demonstrate that combining multidisciplinary large‐scale methods is a relevant approach to better understand volcanic systems, with implications for monitoring strategies.
We analyzed ground deformation prior to the eruptions in 2006 and 2010 at Merapi volcano, Central Java, Indonesia. Ground deformation was monitored by electronic distance measurement (EDM) by measuring the slope distance toward 12 reflectors installed near the summit from five benchmarks on flanks every day. A large change of slope distance (CSD) was detected on the southeast and south baselines and a minor CSD was detected on the north and northwest baselines during the pre-eruptive stages of both the 2006 and 2010 eruptions. We applied a block movement model to the south and southeast baselines and a spherical pressure source model to the CSDs on the north and northwest baselines using the finite element method (FEM). The rates of block movement southward and the volume change of the pressure source increased on April 7, 2006 and continued at constant rates until the appearance of a new lava dome on April 26. Prior to the eruption in 2010, the block movement southeastward and the volume increase of the pressure source accelerated in the middle of October, and acceleration continued until the first outburst on October 26, 2010. Temporal patterns of the block movement and the increase in the volume of the pressure source correlate well with the increase in seismicity of VT and MP earthquakes. The pressure sources were obtained at a depth of 2 +/- 0.5 km below the summit, and this position corresponds to the aseismic zone of VT earthquakes. Magma injection at the shallow part of this region causes an increase in the volume of the pressure source, and inflation of the ground of the summit triggered gravitational slip southeastward or southward of the ground surface. The volumes increases of the pressure sources were 9.7 +/- 1 M m(3) and 17.6 +/- 0.8 M m(3) in 2006 and 2010, respectively. The volume increase is related to the scale and type of the eruption. The effusive eruption in 2006 had a volcanic explosivity index (VEI) of 2 and the explosive eruption in 2010 had a VEI of 4. The directions and amounts of the block movement are strongly related to topography, hydrothermally weak zone and underground gap near the summit between West and East Domes. (C) 2018 Published by Elsevier B.V.
Moment tensor representation is useful for interpreting source processes from seismic and/or geodetic observations. However, there remains difficulty to determine the source processes because some models cannot be distinguished from moment tensors and the range of moment tensor represented by particular process is generally limited. We examine magma movement between two ellipsoids with consideration of mass conservation and pressure balance. The resultant moment tensor component ratio (MTCR) range is clearly distinguished from that generated by expanding ellipsoids. However, there is a range of the MTCR that cannot be explained by either fluid expansion or movement; in such cases, partial pressure recovery after fluid movement provides a plausible explanation. The MTCR for fluid movement between two ellipsoids with randomly varying geometry and relative orientation is shown to exhibit a concentrated trend representing a linear combination of a double‐couple source and an opening crack. Although such an MTCR is generally interpreted as a crack with both tensile and shear dislocations, we demonstrate that it can be generated by fluid movement between two thin ellipsoids. Our results provide a physically reasonable mechanism for interpreting the MTCR obtained from observed data with a comprehensive view of possible source processes.
