Скачать книгу

Underscript i equals 1 Overscript 2 Endscripts sigma-summation Underscript j equals 1 Overscript 2 Endscripts sigma-summation Underscript k equals 1 Overscript 2 Endscripts s left-brace z Subscript k Baseline tangent Superscript negative 1 Baseline StartFraction x Subscript i Baseline y Subscript j Baseline Over z Subscript k Baseline upper R Subscript italic i j k Baseline EndFraction minus x Subscript i Baseline ln left-parenthesis upper R Subscript italic i j k Baseline plus y Subscript j Baseline right-parenthesis minus y Subscript j Baseline ln left-parenthesis upper R Subscript italic i j k Baseline plus x Subscript i Baseline right-parenthesis right-brace EndLayout comma"/>

      where upper R Subscript italic i j k Baseline equals StartRoot x Subscript i Superscript 2 Baseline plus y Subscript j Superscript 2 Baseline plus z Subscript k Superscript 2 Baseline EndRoot and s = s i s j s k with s 1 = −1 and s 2 = +1. This formula provides the gravity kernel written in equation 3.7. When the elevation model is provided in a rectangular grid, the gravitational attraction of the topography is calculated easily just by summing up all the gravitational contributions of prisms in each grid.

Schematic illustration of geometry of a rectangular prism and its gravitational attraction.

       László Oláh and Hiroyuki K. M. Tanaka

       Earthquake Research Institute, and International Muography Research Organization (MUOGRAPHIX), The University of Tokyo, Tokyo, Japan; and International Virtual Muography Institute, Global

      ABSTRACT

      The application of machine learning techniques on remote sensing data has opened new possibilities in automated classification and forecasting of geophysical phenomena, including those triggering natural hazards, such as landslides or volcano eruptions. Muographic imaging can reveal the structural and density changes in geological bodies in connection with subsurface geophysical phenomena. We present the recent progress in machine learning‐based muographic image processing and discuss how it may reveal subsurface volcanic phenomena and contribute to volcano eruption forecasting. Our analysis is based on data sets collected by the Multi‐Wire‐Proportional‐Chamber‐based Muography Observation System during the eruption episodes of Minamidake crater of Sakurajima volcano between October 2018 and July 2020. The area under the curve score of 0.761, provided by the receiver operating characteristic curve analysis of a convolutional neural network, shows correlation between the impending vulcanian eruptions of Minamidake crater and the muographic images captured on the previous days.

      4.1.1 The Concept of Machine Learning

      Machine learning (ML) is a process of automated construction of mathematical models based on training data, in order to accomplish decisions or predictions on test data without being explicitly programmed (Samuel, 1959; Géron, 2019). ML aims to perform the cognitive functions of the human brain for tasks that require a lot of fine‐tuning, or are too complex for conventional techniques and human sense organs, to explore hidden information. The first mathematical model of neurons was developed based on the “all‐or‐none” characteristics of nervous activity (McCulloch & Pitts, 1943). Some years later, the perceptron model was developed for storing and organizing information (Rosenblatt, 1958). Currently, support vector machines (SVM) (Vapnik, 1995) and deep neural networks (LeCun et al., 2015) are the most commonly used ML models.

      In practice, insufficient quality or quantity of training data or wrong algorithm can result in a suboptimal model. The training data are underfitted when the model is too simple to learn the underlying structure of data. The training data are overfitted when the model overgeneralizes those. The regularization of model can be achieved by means of tuning the parameters of learning algorithm, called hyperparameters. The following developments helped the regularization procedure and accelerated the running of complex ML models: (i) Pre‐training of a hidden layer using the unsupervised learning procedure for restricted Boltzmann machine reduces the dimensionality of complex data (G. E. Hinton & Salakhutdinov, 2006). (ii) Rectified Linear Units (ReLU) activation functions of hidden layers drastically speed up the training of deep neural networks (Nair & Hinton, 2010). (iii.) The dropout technique prevents the overtraining of deep neural networks (Srivastava et al., 2014).

      The latest improvements in ML techniques and the rise of graphics processing units revolutionized image and sound recognition (G. Hinton et al., 2012; Krizhevsky et al., 2012), medical imaging and radiation therapy (Esteva et al., 2017; Litjens et al., 2017; Sahiner et al., 2019), computational chemistry (Goh et al., 2017), biology and medicine (Ching et al., 2018), astrophysics (VanderPlas et al., 2012), high energy physics (A Living Review of Machine Learning for Particle Physics, 2020; Guest et al., 2018; Radovic et al., 2018), etc. The improvement of geophysical signal processing techniques is required for automated exploration and classification of various signals from Earth’s subsurface and forecasting of the geophysical phenomena including those triggering various natural hazards, such as volcano eruptions (Reath et al., 2016; Sparks, 2003), earthquakes (Geller, 1997), landslides (Korup & Stolle, 2014), or tsunamis (Titov et al., 2005). Although statistical algorithms have already achieved progress in the prediction of short‐term volcano eruptions (Brancato et al., 2016; Newhall & Hoblitt, 2002), it was natural to also consider using ML techniques. In the following sections, we focus on the prediction of volcano eruptions with ML techniques.

      4.1.2 Volcano Eruption Forecasting With Machine Learning

      Volcano eruption prediction (Sparks, 2003) has three main categories: (i) long‐term eruption prediction aims to forecast future catastrophic events based on the compilation of various scientific and historic data. Although occurrence of catastrophic eruptions is very rare and thus modeling of volcanic mechanism is not possible, this is a highly studied area that is under continuous development. (ii) Mid‐term eruption prediction aims to forecast eruptions that can influence daily human activities. These kinds of eruptions occur more frequently, thus the eruption mechanism can be explored and a model of eruption mechanism can be developed. (iii) Short‐term eruption prediction is based on observation of volcanic phenomena occurring anywhere from a few minutes to a few days before the eruptions. Short‐term eruptions occur frequently and they can have strong impacts on the surrounding areas, e.g. the 2014 eruption of Mount Ontake, Honsu, Japan, caused enormous damage to tourists and local population (Yamaoka et al., 2016).

      Automatic classification of seismic signals originating from active volcanoes (Falsaperla et

Скачать книгу