Simulated Annealing for volcano muography. (August 2021)

• Record Type: Journal Article
• Authors: Vesga-Ramírez, A.; Sanabria-Gómez, J.D.; Sierra-Porta, D.; Arana-Salinas, L.; Asorey, H.; Kudryavtsev, V.A.; Calderón-Ardila, R.; Núñez, L.A.
• Is Part Of: Journal of South American earth sciences. Volume 109(2021)
• Publication Date: 2021-08
• Subjects: 00-01 -- 86-08; Muography -- Muon tomography -- Muon radiography -- Geophysical inversion -- Simulated Annealing -- Volcanoes -- Cosmic ray techniques; Geology -- Latin America -- Periodicals; Earth sciences -- Latin America -- Periodicals; Geology -- Antarctica -- Periodicals; Earth sciences -- Antarctica -- Periodicals; Geology -- Caribbean Area -- Periodicals; Earth sciences -- Caribbean Area -- Periodicals; Géologie -- Amérique latine -- Périodiques; Sciences de la terre -- Amérique latine -- Périodiques; Géologie -- Antarctique -- Périodiques; Sciences de la terre -- Antarctique -- Périodiques; Géologie -- Caraïbes (Région) -- Périodiques; Sciences de la terre -- Caraïbes (Région) -- Périodiques; Earth sciences; Geology; Antarctica; Caribbean Area; Latin America; Periodicals; 550.5
• DOI: 10.1016/j.jsames.2021.103248
• Languages: English
• ISSNs: 0895-9811

Abstract

Abstract: Muography or muon radiography is a non-invasive emerging image technology relying on high energy atmospheric muons, which complements other standard geophysical tools to understand the Earth's subsurface. This work discusses a geophysical inversion methodology for volcanic muography, based on the Simulated Annealing algorithm, using a semi-empirical model of the muon flux to reach the volcano topography and a framework for the energy loss of muons in rock. The Metropolis-Simulated-Annealing algorithm starts from an "observed" muon flux and obtains the best associated inner density distribution function inside a syntetic model of the Cerro Machín Volcano (Tolima-Colombia). The estimated initial density model was obtained with GEOMODELER, adapted to the volcano topography. We improved this model by including rock densities from samples taken from the crater, the dome and the areas associated with fumaroles. In this paper we determined the minimum muon energy (a function of the arrival direction) needed to cross the volcanic building, the emerging integrated flux of muons, and the density profile inside a model of Cerro Machin. The present inversion correctly reconstructed the density differences inside the Machín, within a 1% error concerning our initial simulation model, giving a remarkable density contrast between the volcanic duct, the encasing rock and the fumaroles area. Highlights: Geophysical inversion methodology for volcanic muography based on the Simulated Annealing algorithm. The inversion algorithm correctly reconstructed the density differences inside the Machin, within a 1% error concerning our initial simulation model. We estimate the emerging integrated flux of muons and the density profile inside a model of Cerro Machin Volcano (Tolima, Colombia).