adg vol5 n02 Iuli 431_438.pdf ANNALS OF GEOPHYSICS, VOL. 45, N. 2, April 2002 431 Advanced magnetic visualization of Mt. Vesuvius shallow plumbing system by probability tomography Teresa Iuliano (1), Paolo Mauriello (2) and Domenico Patella (1) (1) Dipartimento di Scienze Fisiche, Università «Federico II», Napoli, Italy (2) Istituto per le Tecnologie Applicate ai Beni Culturali, CNR, Roma, Italy Abstract This paper reports the results of the application of a new magnetic probability tomography imaging method to the aeromagnetic data set collected by AGIP on Mt. Vesuvius in 1978. A magnetic dipole occurrence probability function is defined, in order to contain all information about the whole class of equivalent sources compatible with the original survey data set. An advanced 3D visual system is also applied for the first time to visualise the spacial pattern of the magnetic dipole occurrence probability function. It is shown that in the case of Vesuvius, the shallow plumbing system is characterised by a high probability of occurrence of a dipole at 2 km of depth b.s.l., nearly directed along the present main earth’s magnetic field. This leads to hypothesize that the top terminal portion of the volcanic apparatus is completely filled with magnetised volcanic material. 1. Introduction Probability tomography was first introduced to image self-potential anomaly sources (Patella, 1997a,b) and then extended to analyse geoelectrical (Mauriello et al., 1998; Mau- riello and Patella, 1999a), natural electro- magnetic induction (Mauriello and Patella, 1999b, 2000), gravity (Mauriello and Patella, 2001a,b) and magnetic data (Mauriello and Patella, 2001c). It is an interpretation tool, which accounts for the inherently uncertain nature of the geophysical inversion problem, due to the many different source configurations that can be modelled compatibly with a given data set. Probability tomography has been already applied to the analysis of self-potential, geo- electrical, gravity and magnetic data sets collec- ted in the volcanic area of Mt. Vesuvius (Naples, Italy) (Di Maio et al., 1998; Patella and Mauriello, 1999; Iuliano et al., 2001). The purpose of this paper is to repropose the to-mographic imaging of the magnetic sources beneath Mt. Vesuvius, using a new formulation of the magnetic probability tomography theory (Mauriello and Patella, 2001c) coupled with a 3D object-oriented Advanced Visual System (AVS). In the following sections, the new mag- netic probability tomography imaging method and the most important aspects of a 3D AVS visual development tool are described. Then Mailing address: Prof. Domenico Patella, Dipartimento di Scienze Fisiche, Università «Federico II», Complesso Universitario Monte S. Angelo, Via Cintia, 80126 Napoli, Italy; e-mail: patella@na.infn.it Key words probability tomography – magnetic anomaly – Vesuvius volcano 432 Teresa Iuliano, Paolo Mauriello and Domenico Patella the results of the application to the Vesuvius aeromagnetic data set collected by Cassano and La Torre (1987) are presented and discussed. 2. Probability tomography theory 2.1. The 3D magnetic occurrence probability function We assume a coordinate system with the (x, y)- plane at sea level and the z-axis positive up- wards. Let B(r) be a static magnetic induction field, evaluated at a grid of points r S, where S is a portion of the earth’s surface characterized by a topography function z (x, y). In terms of total magnetization vector M(r), B(r) can be expressed as (Jackson, 1998; Mauriello and Patella, 2001c) where V is a volume containing all the magnetic sources at r V and n is the unit vector in the direction of r-r . Assuming the magnetic field on the earth’s surface as due to a finite distribution of Q elementary magnetic sources, eq. (2.1) can be discretized as follows: where n q is the unit vector in the direction r r q and the q-th generic element is a small volume V q centred at r q with magnetization M(r q ) and magnetic moment d q given by Using the principles of the probability tomo- graphy (Patella, 1997a,b), at first the total power u , associated with the component B u (r) of B(r) along a generic direction identified by the unit vector u, is defined as (2.4) where the explicit expressions of the three uv (r – r q ) functions ( = x, y, z) are given by (2.5a) (2.5b) (2.5c) with n qv ( = x, y, z) being the components of the unit vector nq. Then, a three-component Magnetic Occur- B r n n M r' M r' r r' ( ) = ( )[ ] ( )µ0 3 4 3 dV V B r n n d d r r q q q q q ( ) = ( ) = 3 1 3 q Q d M rq= ( ) µ0 4 dVq q . u u S B dS= ( ) = ( ) 2 r S d B dSqv u uv Sv x y zq Q = ( ) ( ) ( )==1 r r rq , , uy( ) =r rq = 1 r rq 3 3n n n nqy qx qy qzi u j u k u j u+ +( )[ ] uz r rq( ) = n n n nqz qx qy qz r r i u j u k u k u q = + +( )[ ] 1 3 3 ux( ) =r rq n n n nqx qx qy qz= + +( )[ ]3 1 3 r r i u j u k u i u q (2.1) (2.2) (2.3) 433 Advanced magnetic visualization of Mt. Vesuvius shallow plumbing system by probability tomography rence Probability (MOP) function uv (M)(r q ) is defined as (Mauriello and Patella, 2001c) ( = x, y, z), (2.6) where ( = x, y, z), (2.7) and g (z), called the topographic surface regularization factor, is (2.8) Each uv (M) function satisfies the bounding condition 1 uv (M)(r q ) +1 and three uv (M) values ( = x, y, z) can be computed at each r q . Each value is interpreted as the probability, with which the homologous component M ( = x, y, z) of a M-field can at r q be considered responsible for the measured component for the magnetic field. In geophysical exploration, either the z- component or the modulus of the earth’s magnetic field is usually measured and a scalar secondary field (the anomalous field) is evaluated in order to identify local sources of magnetic anomaly. In the case of the so-called total field survey, the scalar anomalous field is obtained by subtracting from the measured modulus of the earth’s magnetic field the uv M( )( ) =rq uv M u Y Y X X uvC B g z dxdy ( ) = ( ) ( ) ( )r r rq , modulus of the known primary magnetic field. Since the secondary field is always a very small fraction of the primary field, it can be readily shown that in the case of a total field survey the scalar secondary field is the projection of the secondary vector field along the direction of the primary vector field (Blakely, 1996). This direction can generally be assumed uniform within the areas normally considered in geophysical prospecting (Parasnis, 1997). Thus, from now on we will consider B(r) as the secondary vector field, and the scalar com- ponent of B(r) along any fixed direction as the object of analysis in the new tomographic approach. 3. 3D probability tomography procedure The 3D tomography procedure for imaging the sources of a magnetic field, measured on a generally non-flat topography, consists in a reiterated computation code, involving the scanner functions u and the B u (r) field data set. In practice, since we do not know the position of the real sources generating the anomalous magnetic field, we use a synthetic source of unitary strength to scan the x, y, z half-space below the survey area (tomospace), in order to search where the real sources can be located in probabilistic sense. The scanning operation is performed by computing the cross- correlation integrals in eq. (2.6) for each point (x q ,y q ,z q ) of a regular grid within the tomospace. At each point, the value of each integral is interpreted as the probability of occurrence of the relative magnetic source component, whose positive or negative orientation depends on whether it is > 0 or < 0. By scanning the tomospace, using, e.g., a sequence of horizontal slices spaced from each other by a constant depth interval, we can finally obtain a 3D image of the equivalent magnetic source distribution underground in a probabilistic sense. In order to improve the filtering property of the scanning procedure, for each r q of the tomospace it is advisable to use varying sizes of the integration surface in eq. (2.6). The smallest size is the [– X, X ] × [– Y, Y ] domain fully uv MC( ) = u Y Y X X uv Y Y X X B g z dxdy g z dxdy= ( ) ( ) ( ) ( )r r rq / 2 2 1 2 ,, g z z x z y .( ) = +( ) +( ) 1 2 2 434 Teresa Iuliano, Paolo Mauriello and Domenico Patella containing the surface trace of the magnetic response of the scanning source element placed at r q . The greatest surface is, of course, the largest rectangle fitting to the whole survey area. The highest (r q ) is taken with its sign as the most appropriate source occurrence probability at r q . 4. Advanced data visualization In the most general term, data visualization is the process of representing data as a graphic display on a computer screen. In practice, data visualization gives an insight into data allowing the interpreter to graphically view structure and changes in the data. Generally speaking, there are two things the interpreter wants to do with data: analyse it to obtain information and present it to others to share that information. In either case, an effective visualization is an important tool. An Advanced Visual System (AVS) extends the obvious advantages of standard visualization by sophisticated 3D representations of data to be developed. Modern 3D AVS softwares are considered object-oriented, visual development tools, enabling one to build reusable application components and complex visualization ap- plications. AVS softwares become very useful when they can support different object-oriented techniques, such as encapsulation of data and methods, class inheritance, templates and instances, object hierarchies and polymorphism. Very shortly, a 3D AVS tool can be considered a visual development environment that can be used to connect, define, assemble and manipulate objects and related applications through mouse- driven operations. The objects and applications that are connected and assembled control how data is processed and how it is displayed. The interpreter can even add a user interface to create a complete application that can be delivered as a stand-alone application. We have applied a 3D AVS software (Advanced Visual Systems Inc., 1998) to vis- ualize the results of the application of the magnetic probability tomography to the volcanic area of Mt. Vesuvius. Of course, as we cannot reproduce all the dynamical capabilities of the 3D AVS, we show in the following section only a particular perspective of the tomo- graphic images from the most effective point of view. 5. Mt. Vesuvius magnetic tomography and AVS imaging Vesuvius is considered one of the most risky active volcanoes in the world, because of its volcanic history and closeness to the city of Naples, Italy. The results achieved so far (Iuliano et al., 2001, and references therein) indicate that the shallow part of the volcano is made of a unique central plumbing system, entirely filled of altered volcanics in the summit portion. This peculiarity makes the volcano extremely hazardous, since explosion becomes a highly probable event in case of renewal of extruding activity. The magnetic study can notably help to ascertain whether such a condition really occurs. To this purpose, we show the results of the application of the probability tomography to the aeromagnetic total field data set collected in 1978 by the Italian Oil Company AGIP to study the volcano-geothermal structure of Mt. Vesuvius. A cesium optical pumping magnetometer was used at the constant altitude of 1460 m a.s.l. Figure 1 shows the original total magnetic field map (after Cassano and La Torre, 1987). Subtracting a constant regional field in the survey area with intensity 30 100 nT, inclination 56.30° and declination 1°, the residual effect is a large and intense positive anomaly with peak value of approximately 1800 nT, located in the central part of the volcanic system. The 3D AVS pictures in fig. 2a-c show the pattern of the MOP functions ux (M), uy (M) and uz (M), respectively, drawn by plotting the MOP isosurface of the lower limit of the coloured palette corresponding to the decimal step where the highest absolute MOP values occur. Only absolute MOP values equal or greater than 0.6 have been considered in fig. 2a-c. In particular, fig. 2c shows a well defined nucleus centred at about 2 km b.s.l., including the lowest computed values of uz (M), falling in the range from – 0.9 to – 1. Accordingly, fig. 2b shows a smaller nucleus, centred again at about 2 km b.s.l., including the highest computed values of uy (M), falling in the 435 Advanced magnetic visualization of Mt. Vesuvius shallow plumbing system by probability tomography 0 5 10 15 20 25 x (km) 0 5 10 15 20 25 30 35 y ( k m ) 29900 30500 31100 31700 Total magnetic field (nT) Castellamare S.Anastasia Mt Somma Vesuvius Torre Annunziata Pomigliano Naples Trecase well Fig. 1. The total field aeromagnetic map in the volcanic area of Mt. Vesuvius (after Cassano and La Torre, 1987). 436 Teresa Iuliano, Paolo Mauriello and Domenico Patella z (km) -1 -3 -5 -7 z (km) -1 -3 -5 -7 a b ux (M) 0.6 0.6 uy (M) 0.7 1 uz (M) 0.9c Fig. 2a-c. 3D AVS imaging of the magnetization oc- currence probability function along the x-axis (a); y-axis (b) and z-axis (c). The top slice is the exper- imental map of the residual total magnetic field. z (km) -1 -3 -5 -7 range from 0.6 to 0.7. No signal above the threshold of 0.6 is present in the ux (M) tomography of fig. 2a. Combining all these features, the relevant information is the high occurence probability of an inclined dipole centred at x = 12 km, y = 17 km and z = 2 km, with dominant negative vertical component and positive hor- izontal component directed along the y-axis. Figure 3a,b shows the comparison between the experimental map of the residual total magnetic field (fig. 3a) and the synthetic map (fig. 3b) obtained assuming a magnetic dipole centred at x = 12 km, y = 17 km and z = 2 km, with magnetic moment of 110 nT, inclination 56.30° and declination – 1°. The misfit between the two maps is in modulus less than 2%. Therefore, the synthetic dipole can represent a simple model of the observed residual magnetic anomaly field in the Vesuvius volcanic area. Concluding, this leads to hypothesize the existence of a uniformly magnetised material filling the top terminal portion of the volcano, exactly where gravity data indicate a low density (Iuliano et al., 2001). 437 Advanced magnetic visualization of Mt. Vesuvius shallow plumbing system by probability tomography REFERENCES ADVANCED VISUAL SYSTEMS INC. (1998): AVS/Express manual. Waltham, MA, USA. BLAKELY, R.J. (1996): Potential Theory in Gravity and Magnetic Applications (Cambridge University Press), pp. 461. CASSANO, E. and P. LA TORRE (1987): Geophysics, in: Somma Vesuvius, edited by R. SANTACROCE, Quad. Ric. Sci., 114/8, 175-196, Consiglio Nazionale delle Ricerche, Rome. DI MAIO, R., P. MAURIELLO, D. PATELLA, Z. PETRILLO, S. PISCITELLI and A. SINISCALCHI (1998): Electric and electromagnetic outline of the Mount Somma-Vesuvius structural setting, J. Volcanol. Geotherm. Res., 82, 219-238. Fig. 3a,b. Comparison between the observed map of the residual total magnetic field (a) and the synthetic map (b) obtained assuming a magnetic dipole centred at x = 12 km, y = 17 km and z = 2 km, with magnetic moment 110 nT, inclination 56.30° and declination – 1°. 0 5 10 15 20 25 x (km) 0 5 10 15 20 25 30 35 y ( k m ) 29900 30500 31100 31700 Total magnetic field (nT) Castellamare S.Anastasia Mt Somma Vesuvius Torre Annunziata Pomigliano Naples Trecase well y ( k m ) 0 5 10 15 20 25 x (km) 0 5 10 15 20 25 30 35 Castellamare S.Anastasia Mt Somma Vesuvius Torre Annunziata Pomigliano Naples Trecase well a b IULIANO, T., P. MAURIELLO and D. PATELLA (2001): Looking inside Mount Vesuvius by potential fields integrated probability tomographies, J. Volcanol. Geotherm. Res. (in press). JACKSON, J.D. (1998): Classical Electrodynamics (J. Wiley and Sons, New York), pp. 832. MAURIELLO, P. and D. PATELLA (1999a): Resistivity anomaly imaging by probability tomography, Geophys. Prospect., 47, 411-429. MAURIELLO, P. and D. PATELLA (1999b): Principles of probability tomography for natural-source electro - magnetic induction fields, Geophysics, 64, 1403-1417. MAURIELLO, P. and D. PATELLA (2000): A physical pattern recognition approach for 2D electromagnetic induction studies. Ann. Geofis., 43 (2), 343-359. MAURIELLO, P. and D. PATELLA (2001a): Gravity probability 438 Teresa Iuliano, Paolo Mauriello and Domenico Patella tomography: a new tool for buried mass distribution imaging, Geophys. Prospect., 49, 1-12. MAURIELLO, P. and D. PATELLA (2001b): Localization of maximum-depth gravity anomaly sources by a distribution of equivalent point masses, Geophysics, 66, 1431-1437. MAURIELLO, P. and D. PATELLA (2001c): Magnetic buried sources discrimination by probability tomography, Geophys. Prospect. (submitted). MAURIELLO, P., D. MONNA and D. PATELLA (1998): 3D geoelectric tomography and archaeological applications, Geophys. Prospect., 46, 543-570. PARASNIS, D.S. (1997): Principles of Applied Geophysics (Chapman & Hall, London), pp. 429. PATELLA, D. (1997a): Introduction to ground surface self- potential tomography, Geophys. Prospect., 45, 653-681. PATELLA, D. (1997b): Self-potential global tomography including topographic effects, Geophys. Prospect., 45, 843-863. PATELLA, D. and P. MAURIELLO (1999): The geophysical contribution to the safeguard of historical sites in active volcanic areas. The Vesuvius case-history, J. Appl. Geophys., 41, 241-258.