Muography. Группа авторов
of data collection. Dynamics of geofluids such as magma, natural gas, underground water, and sea current will cause time‐dependent variations of the subterranean densimetric heterogeneity. These variations can be captured by time‐sequential muographic observations.
Figure 1.3 Photographs of airborne muography. Distant (a) and close‐up (b) views of the measurements are shown. The white inset in (a) indicates the region of the close‐up view in (b). The distance between the rotor blade and the cliff is a few meters. Exterior (c) and interior (d) views of the apparatus, an integrated form of an array of detectors and a helicopter, are shown. A highly qualified maneuvering skill is required for this operation.
1.3 PIONEERING WORKS
1.3.1 Early Works
In the early stage of cosmic ray studies, underground measurements were the most effective way to extend the energy range of the measured muon spectrum beyond 1 TeV. In these measurements, mine galleries located at various depths were utilized to measure the depth‐dependent muon flux since geological features of these mines were well studied. Inversely, if this depth‐dependent muon flux was used as a reference curve, the average density above the detector could be derived. The idea of using muons produced by cosmic rays as probes was first applied 75 years ago by E.P. George, who measured the thickness of the rock overburden above a tunnel of the hydroelectric plant in Snowy Mountain, Australia (George, 1955). George measured the reduction in the muon flux after passing through the rock. The apparatus consisted of Geiger counters but was unable to provide an image of any structure within the overlying rock.
The use of muography to reveal the internal structures of inaccessible objects in archeology was introduced in the late 1960s by a group led by Luis Alvarez to search for undiscovered chambers in the Pyramid of Chephern in Egypt (Alvarez et al., 1970). After the invention of a spark chamber with a digital readout, imaging became more realistic. The group recorded the trajectories of muons through the pyramid and studied their transmission to establish why the Pyramid of Chephren had only one burial chamber (the so‐called Belzoni chamber) while the pyramid made by his father Cheops had more complicated internal structures, including the King's and Queen's chambers and the Grand Gallery. This joint project between the United States and the United Arab Republic began in 1966. Calculations showed that if a hidden room was located above the Belzoni chamber, it could be observed in the same way that a void is highlighted by a darker area in an X‐ray image. However, no hidden rooms were found inside the pyramid. Instead, the group showed the potential of the technique by detecting the cap rock from the Belzoni chamber. This pioneering work paved the way for the application of muography in various fields.
Since TeV muons penetrate kilometric rock, the technique shown by Alvarez et al. (1970) was in principle applicable to mountains. This possibility was explored by focusing on detecting muons that traversed at angles almost parallel to the ground surface, which could be utilized to probe mountains by tracing the trajectories of muons emerging from the other side of the mountain (Nagamine et al., 1995). Muography cannot image the deep structure of a volcano such as the magma chamber; however, it can image shallow regions of the volcano, which can provide useful information for understanding how the eruption style might change. The first muographic image of the inside of a volcano suggested possible pathways for magma ejection by visualizing the shape and size of low‐density regions under a deposit of solidified magma (Tanaka et al., 2007). At the same time, the results showed visual evidence of the resolving power of muography, and its applicability to any targets smaller than volcanoes. The first time‐sequential muographic images that captured the motion of subterranean geofluid targeted the rainfall‐triggered permeation of water into the mechanical fracture zone of the active seismic fault (Tanaka et al., 2011). The results later motivated researchers to apply muography to monitoring underground water conditions (Tanaka & Sannomiya, 2012) and magmatic motion inside volcanoes (Oláh et al., 2019; Tanaka et al., 2014).
1.3.2 Magmatic Convection
By taking advantage of the resolving power of muography, we can address the following issues in volcanology: (i) conduit diameter (Tanaka et al., 2007a, 2007b, 2008; Tioukov et al., 2019) and three‐dimensional conduit location (Tanaka et al., 2010); (ii) depth of magma degassing (Tanaka et al., 2009a); (iii) magma convection and magma supply rate (Shinohara & Tanaka, 2011; Tanaka et al., 2009a); (iv) whether the magma pathway is plugged (Tanaka et al., 2007b; Tanaka & Yokoyama, 2008) or drained‐back (Kusagaya & Tanaka, 2015a; Tanaka et al., 2007a); (v) characterization of a high‐density spine inside a porous lava dome (Tanaka, 2016) or the magma intrusion underneath the volcano (Kusagaya & Tanaka, 2015b); (vi) the level of the magma head (Tanaka et al., 2014); (vii) volcanic plug explosion (Tanaka et al., 2009b) and formation (Oláh et al., 2019); (vii) lava dome explosion (Tanaka & Yokoyama, 2013); (ix) tephra deposition (Tanaka, 2020a); and (x) eruption forecast (Nomura et al., 2020).
Figure 1.4 Muographic image of Satsuma‐Iwojima volcano, Japan. The arrow indicates the location of the bubbly low‐density magma.
An example of how a muographic image can be used to understand the eruption dynamics is shown in Fig. 1.4. Statistical errors ranged from 0.02 to 0.2 g/cm3 above 250 m a.s.l. with the variation depending on the location in the volcano. Satsuma‐Iwojima, Japan, continuously emits large amounts of magmatic gases without a significant output of magma. The density gradually decreases up the conduit, and the top of the magma column at 400 m a.s.l. has the lowest density, indicating the presence of magma degassing, in agreement with the magmatic convection model (Shinohara & Tanaka 2011; Tanaka et al., 2009a). In this convection model, a magma conduit is connected to a deep magma chamber (Fig. 1.5), and in the upper part of the conduit, the gas escapes from the magma and exits the volcano. The degassed magma sinks, and, at the same time, new low‐density non‐degassed magma ascends from the bottom of the conduit and the cycle continues.
Degassed magma, which has a high proportion of bubbles, has been interpreted as being the low‐density region, and its dimensions (location and diameter) were compared with the following results of field measurements and laboratory and numerical modeling studies (Shinohara & Tanaka 2011). (i) The depth of the magma head observed 200 m below the crater floor was consistent with the degassing pressure of the magma and had a value of 0.5–3.0 MPa in Satsuma‐Iwojima (Kazahaya et al., 2002). (ii) High‐temperature (> 900 oC) volcanic gasses were continuously emitted from the vent (Shinohara & Tanaka, 2011). (iii) The oversaturation of volatiles in the melt was found, i.e., the degassing had been occurring under relatively low‐pressure conditions (Hedenquist et al., 1994). (iv) Gravity mapping of Satsuma‐Iwojima (with a residual profile derived from Bouguer anomalies assuming a density of 2.0 g/cm3) revealed that a low‐gravity region was located within the volcanic cone, which reached its lowest value in the crater region (Komazawa et al., 2005). The overall features of these results are consistent with the interpretation of this muographic result. Shinohara and Tanaka (2011) compared their muographic results with the results of gravity surveys of Masaya volcano, Nicaragua, and concluded that the low‐density region at the top of the magma column is a common feature of conduit magmatic convection. Conduit magmatic convection is often considered to be the degassing mechanism of basaltic volcanoes (Aiuppa et al., 2009; Burton et al., 2007; Oppenheimer et al., 2009; Shinohara & Witter,