ANNALS OF GEOPHYSICS, 60, 2, 2017, G0219, doi: 10.4401/ag-7062 The italian magnetic repeat station network: results from the 2012.5 ‘Reduced network’ completion Guido Dominici1,*, Antonio Meloni1 1 Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy Article history Received October 26, 2016; accepted December 7, 2016. Subject classification: Earth’s magnetic field, Geomagnetic map, Magnetic surveys, Magnetic repeat stations, Secular variation, Global and regional models ABSTRACT In Italy INGV (Istituto Nazionale di Geofisica e Vulcanologia) has systematically undertaken the task of making periodic measure- ments of the Earth’s magnetic field on a network of more than 110 points (repeat stations), with an average spacing around 55-60 km. Measurements are repeated regularly every 5 years and the last publi- shed data reports and magnetic maps, refer to 2010.0, are available in Dominici et al, 2012. At that time the report referred to a survey of 131 repeat stations (including 2 observatories, 11 stations in Alba- nia, 3 stations in Corsica and 1 in Malta) carried out between 2009 and 2010, also with the purpose of updating the national magnetic cartography. At the epoch 2012.5 a selection of stations, from the first order Italian Magnetic Network, selected on the basis of the lowest amplitude values of anomaly with respect to a ‘normal’ field, was repeated. The number of the selected stations amounts to 25, distri- buted according to an almost uniform national geographical covera- ge, with an average spacing around 100 km. Secular variation and analytical expressions, such as second order polynomials, in latitude and longitude for all field elements, were determined and coefficients were obtained for the spatial field variation and secular variation. We describe here the characteristics of this reduced network with the data elaboration procedure, normal fields and maps, and compare the results with other magnetic field models 1. Introduction Magnetic surveys are normally undertaken for magnetic maps compilation, and associated data ela- borations: in particular for main field and magnetic anomalies representation, for observation of geoma- gnetic field elements secular variation and for the com- putation of regional magnetic models. All surveys are based on two typical magnetic networks, usually cal- led first order (repeat stations) and second order. Re- peat stations are devoted to main field representation; they are constituted of accurately located points at the Earth’s surface, where the three components of the ge- omagnetic field are regularly measured to the highest possible accuracy [see for example: Newitt et al. 1996, Lanza and Meloni 2006]. Regional magnetic models, also called normal fields, [see Haines et al. 1990] are computed on the basis of repeat stations to represent the main geomagnetic field elements pattern over a certain area. Second order networks are denser mea- surement point networks that are generally finalized to improve crustal magnetic anomalies representation. The magnetic characteristics of the Italian terri- tory have been the subject of several investigations. A brief history of Italian magnetic surveys dates back to the seventeenth century and includes the first survey made in 1640, when Declination (D) was measured over 21 stations by Fathers Borri and Martini. In more recent times the magnetic survey of 1881-1892 by Chi- stoni and Palazzo can be considered the first ‘modern’ 3-component national magnetic survey in the unified kingdom of Italy. At that time D, I (Inclination) and H (Horizontal Intensity) were measured on 284 sta- tions by the Ufficio Centrale di Meteorologia e Geo- dinamica (UCMG). More detailed information on this and older magnetic measurement compilations, can be found in Cafarella et al. [1992a], [1992b]. For what concerns recent magnetic networks in Italy a fundamental starting point was the 1979.0 survey. In that occasion a repeat station network of about 120 data points was assessed [Molina et al. 1985] and moreover a full second order magnetic data network, with a higher spatial density of mea- surements, was deployed. This last one consists of horizontal intensity (H), vertical intensity (Z) and to- tal field (F), with measurements taken in the period 1977-1981 on 2552 stations within the framework of the Finalized Project Geodinamica (PFG-CNR), see Molina et al. [1985]. The 1979.0 network (all above mentioned measu- rements were reduced to this central epoch) has ge- G0219 DOMINICI ET AL. nerated a full representation of magnetic elements, a magnetic maps compilation and a set of coefficients for normal field regional models to represent the main magnetic field elements pattern over Italy. For what concerns Declination, this was measu- red over 1529 points and measurements were, under- taken in the 30s of nineteen hundred by the Military Geographical Institute (IGM), for military demands. This set of declination measurements were successi- vely updated for a publication released in 1973.0. Infor- mation on this network can be found in Talamo [1975]. In later years national magnetic surveys were all dedicated to the repetition of the first order network of about 120 repeat stations including 2 Observatories. From 1994 onwards surveys also included 11 stations in Albania, from 2005 3 stations in Corsica and from 2010 1 in Malta, as shown in Figure 1 that reports situation at 2010.0. The main goal of these studies was the determi- nation of secular variation and consequently the main field generated in the Earth’s core. These most recent surveys have been described in Meloni et al. [1988, 1994], Coticchia et al. [2001] and for the years after 2005.0 Do- minici et al. [2007, 2012], that can be considered the la- test surveys at the national level. All magnetic stations, constituting the so-called second order network, have been used since then, only updated on the basis of the first order network (repeat stations) results in order to bring up to date the magnetic cartography. This was done for all the latest years at 2000.0, 2005.0, 2010.0. [De Santis et al. 2003, Dominici et al. 2012]. If we compare Italian repeat station density (an average of 1 station over 3000 km2) with similar Eu- ropean magnetic networks, we find that other nations usually rely on a sensibly smaller density of repeat stations. For example France relies on about 1 station over 20,000 km2 [Mandea 2004], Spain on 1 station over 12,500 km2. Several workshops in the European ma- gnetic operators community were convened to discuss problems regarding magnetic measurements, in order to exchange experiences and coordinate the individual efforts, and also to discuss what measurements density would be considered optimal for a national underta- king. One successful initiative at the European level is known as MagNetE (Magnetic Network of Europe) and was supported by a resolution of IAGA (Interna- tional Association of Geomagnetism and Aeronomy). MagNetE workshops were organised every two years, in order to ensure continued cooperation between the participants from about 20 European countries. The latest workshops were held in May 2011 in Rome, in Prague in 2013 and in Budapest in 2015. In this paper we present the results of a recent ma- gnetic survey undertaken in 2011-2012 on a selection of the above mentioned repeat stations of the first order Italian Magnetic Network. The number of these se- lected stations amounts to 25, compared to the about 120 stations typically measured in the last 40 years or so. In the next paragraphs we will describe this ‘redu- ced’ network, how points were selected, results of the measurements, and we will discuss advantages and di- sadvantages of this reduced density magnetic network. 2. Magnetic repeat station distribution on a ‘reduced’ network Geomagnetic repeat station surveys are mainly de- voted to secular variation knowledge (from the short time scale 2-5 years onwards) and the understanding of its pattern over the interested area, usually starting with the national area. With this information available, local magnetic models and magnetic maps can be updated. The average distance of magnetic stations, the time interval for their reoccupation and all the pro- cedures needed for a full measurement of magnetic field elements, are the basis of the Newitt et al report [Newitt et al. 1996] where all details for magnetic re- peat stations surveys are described. This report and subsequent discussions and agreements, have brought to the suggestion that an average distance between stations of 125 km and a 2 year time interval for reoc- cupation, are considered satisfactory in almost all ca- Figure 1. The original location of 131 repeat stations in Italy including 2 Observatories. The latest full survey undertaken at 2010.0 also inclu- ded 11 stations in Albania, 3 stations in Corsica and 1 in Malta. Circles indicate Observatory and variometric station areas of influence. 2 ITALIAN REDUCED MAGNETIC NETWORK ses; these parameters are also considered to be a good compromise between costs related to the survey ope- rations and value of the results. The same would also be true if we consider the survey results accuracy. In fact a good survey should not be completed in a too long time window, since secular variation keeps going on during the whole effort. It must be considered that the time window should be kept as short as possible in order to make minimum correction for secular varia- tion; furthermore scientific institutions generally can rely only on one or two teams in the field, so one year total measurement time is considered acceptable. In the ‘reduced’ network described here we star- ted taking in consideration that first order magnetic surveys must rely on stations whose magnetic elemen- ts, once measured, should report the magnetic value representative of the main field for that determined re- gion. This is because magnetic measurements provide not only a measurement of the earth’s core field (at lar- ge the main field) but also an insight into the magnetic structure of the Earth’s crust, that is closely connected to its geological structure. For this reason the station selection should be done in order to include only sta- tions representative of the main field. The need to co- ver almost uniformly all the national area was also a constrain. Following above said criteria at the end the number of the selected stations amounts to 25, distri- buted according to a national geographical coverage, with an average spacing around 100 km (Figure 2). For what concerns field instrumentations, also in 2012.5 survey a DIM (Declination Inclination Magnetometer) a portable theodolite fitted with a fluxgate magneto- meter, was used for the angular measurements of the Earth’s magnetic field (therefore declination D and in- clination I, with 0.1 minute of arc accuracy) and a pro- ton precession magnetometer for the measurement of the field intensity F (0.01 nT accuracy). All other ma- gnetic elements can be derived. Survey was repeated on known repeat stations and in the general case azi- muth were already determined and known. However in a few cases, stations could not be reoccupied (i.e. a station mark was lost) and a new mark was manu- factured onsite. In these cases azimuth determination was made anew; this was also done on those stations in which the visibility of the old azimuth marks was limited or absent. A gyroscopic theodolite was used in all these cases. This instrument was preferred to Global Positioning System (GPS) instruments because it can also be used for measurements in tunnels or roof cove- red places. Moreover even where, in presence of bushes of tall stem, GPS could show problems for the satellites limited visibility. In addition the two instruments show comparable precisions of tenth of minutes. For each repeat station a monograph is realized. In the monograph the site is described, instructions are given to easily reach it, what land marks are useful for the approach, what azimuths are visible and all infor- mation that could be useful for identification. 3. Magnetic measurement reduction procedures, normal fields and maps Magnetic elements observed at the repeat stations are reduced firstly to 02 UT using diurnal variation cor- rection from digital data obtained at magnetic observa- tories working at the time of measurements in the field. For what concerns the 2012.5 survey these were Duronia (Dur) and Castello Tesino Observatory (Cts). Portable variometric stations were also used for rapid magnetic variation corrections, when stations were considered far away from the two indicated Observatories; in this case a portable time variation fluxgate station was installed (all variometric stations are reported in Figure 2). The value of a generic element E (i.e., D,I or F) Figure 2. The geographical distribution of the 25 reduced network stations as described in this paper for the 2012.5 reduced network. Colored circles indicate Observatory and variometric station areas of influence. Pink refers to CTS Observatory area of influence: Yellow refers to DUR Observatory area of influence. Green refers to Variometric OTT station area of influence; Light blue refers to Variometric GIB station area of influence. In case of overlapping regions the two variometric records were considered. 3 DOMINICI ET AL. at station s, i.e. reduced, at 02UT was calculated fol- lowing the: Es(02UT)d = Evar (02UT)d+ [Es(t) d - Evar(t)d] Where: Es(02UT)d = Value of element E at station s reduced at time 02UT of day of measure d Evar (02UT)d = Value of element E at variometer station s at time 02UT of day of measure d Es(t) d = Value of element E observed at station s at time t at that day d Evar(t)d = Value of element E at variometer station at time t at that day d Secondly data are reduced to Castello Tesino Obs. for the fixed epoch, i.e. 2012.5 following the: Es (2012.5) = Eobs (2012.5)+ [Es (02UT) d - Eobs (02UT)d] Where: Es (2012.5) = Value of element E at station s reduced at epoch 2012.5 Eobs (2012.5) = Mean Value of element E at Observatory at epoch 2012.5 Eobs (02UT)d = Value of element E measured at Obser- vatory at time 02UT of day of measurement d. The classical technique of SHA (Spherical Har- monic Analysis) represents mathematically the con- figuration of the Earth’s magnetic field in terms of spherical harmonics. This can be applied only to the entire Earth’s globe and the characteristic minimum wavelength associated with this technique, used for example in IGRF (International Geomagnetic Refe- rence Field), when maximum degree is 12, is around 4000 km [IAGA 2010]. This is the reason why, on a li- mited part of the globe, regional magnetic models are normally used. In the Italian area different techniques were used during the years. In many cases simple poly- nomial second order, together with other methods, for example SCHA (Spherical Cap Harmonic Analysis) see De Santis et al. 1997, and other related techniques, were applied at the local scale. In this paper a regional magnetic model was con- structed using a mathematical second order polyno- mial, that we call ‘normal field’. At the scale of the case here discussed we expect that the normal field inclu- des almost all the contribution from the Earth’s core (equivalent to a SHA with maximum degree 12-14). However it cannot be excluded that a deep lithosphe- ric contribution could be involved and represented in the normal field coefficients [see also Korte and Lusur 2012]. The normal field was computed independently for each element and results, for a generic element E, from the expression: E(ϕ,λ) = ao + a1ϕ + a2λ + a3ϕ 2 + a5ϕλ ϕ and λ are latitude and longitude respectively and ai are the coefficients. These coefficients were compu- ted by means of a least squares fit over the observa- tional values. When the total number of repeat sta- tions was used, as for all latest surveys (for example 2010.0, 2005.0, etc…) in order to avoid contamination from stations located in anomalous areas, the Chau- venet rejection criterion was generally used for remo- ving anomalous values from the computation of the least squares fit coefficients. Essentially the mean re- sidual of each individual station was calculated using the expression σ =± υυ[ ] / (n−6) where [υυ] is the sum of the square differences between observed and computed values, n is the number of stations used in the inversion and 6 is the number of coefficients. Stations with residuals greater than 2σ were neglected and the calculation repeated for the remainder. This Table 1. Coefficients of the normal geomagnetic field in Italy for 2012.5; ϕ = (Lat. – 42°) in minutes λ = (Long. – 12°) in minutes D, I in minutes F, H, Z in nT D (‘) = 145.28 + 0.03227 ϕ + 0.19034 λ - 0.00007 ϕ² - 0.00005 λ² + 0.00014 λϕ F (nT) = 46273.7 + 5.70875 ϕ + 1.17744 λ - 0.00191 ϕ² + 0.00055 λ² - 0.00020 λϕ H (nT) = 24390.3 - 9.46341 ϕ - 0.26936 λ - 0.00017 ϕ² + 0.00070 λ² + 0.00061 λϕ Z (nT) = 39319.2 + 12.65370 ϕ + 1.56750 λ - 0.00480 ϕ² + 0.00027 λ² - 0.00090 λϕ I (‘) = 3490.91 + 1.09820 ϕ + 0.07938 λ - 0.00031 ϕ² - 0.00003 λ² - 0.00010 λϕ 4 ITALIAN REDUCED MAGNETIC NETWORK procedure was repeated until no residuals greater than 2σ were left. For the 2010.0 survey the stations ne- glected and the mean residual σ for all the remaining stations were 39 and σ = 2.535’ for D, 31 stations and σ = 25.1nT for F, 25 stations and σ = 22.18nT for H, 28 stations and σ = 27.1nT for Z. In the case of the reduced network (at 2012.5), reported in this paper, this criterion was not applied. This point was taken in consideration since the Chauvenet criterion is a means of assessing whether one experimental data (an out- lier) from a set of observations, is likely to be spurious; we considered that this should not be applied on an objectively more limited data set. It was considered that deletion of outlier data should not be done espe- cially in small sets where a normal distribution cannot be assumed. Coefficients of the normal geomagnetic field in Italy for 2012.5 are reported in Table 1. In order to test the validity of a reduced network approach, at least in the case of the Italian area, we made a comparison of the normal polynomial field in- terpolation, called here INGRF with the IGRF for total intensity F. These comparisons are reported in Figure 3. For this test we had first compared the normal field for the epoch 2010 (INGRF), when all repeat stations were available (panel a) and afterwards using only the 25 stations (INGRF_red), that were successively used for the 2012.5 reduced network (panel b), both models were compared to IGRF. For the same epoch (2010.0) we found that the use of a limited number of stations (that were afterwards used for the reduced network) had not introduced sensibly larger differen- ces with respect to the use of the full set of available repeat stations data. In panel c the only available 25 reduced network station were used for a new polyno- mial normal field (INGRF_red) at 2012.5. This last pa- nel shows almost same patterns with respect to IGRF 2012.5 similarly to the 2010.0 case. All panels show larger differences at the edges and better agreements at the center of the considered areas. This effect can be interpreted as a larger ability of polynomial models to represent the earth’s core field at the map’s center than at the edges. We conclude that a comparison of the two 2010.0 models, obtained from the original and reduced networks, shows a moderate difference across all the Italian area. 4. Conclusions In this paper we presented the results of a ma- gnetic survey undertaken in 2012 on a selection of the Italian network of magnetic repeat stations (first or- der network). The original number of stations for the previous full 2010 survey was 131 including 2 Obser- vatories including 11 stations in Albania, 3 stations in Corsica and 1 in Malta, as shown in Figure 1. From this network some stations were selected for the 2012.5 that was designed as a ‘reduced’ survey. We have de- scribed how the reduced network stations were chosen and we stressed that all reduced network points were selected on the basis of their minimum deviation from the expected earth’s core magnetic field, and conside- ring the geographical distribution in order to cover uniformly almost all the Italian area. At the end the number of selected stations amounts to 25, distributed over the Italian area with an average spacing around 100 km (Figure 2). Regional modeling of magnetic data from a repe- Figure 3. At 2010.0 Difference represented with a colour scale between F normal field polynomial interpolation from INGRF on all stan- dard repeat stations, and 2010.0 F from IGRF (panel a). At 2010.0 F normal field polynomial interpolation INGRF_red computed on only 25 stations that were selected later for the 2012.5 survey and IGRF (panel b). Finally In panel c the comparison made at 2012.5 between IGRF and INGRF_red from 2012.5 measurements. 5 DOMINICI ET AL. at station network was also the mathematical techni- que used in the case of this reduced network at 2012.5, in order to represent the main field (essentially the core field) as accurately as possible at the local scale. In this study we used only simple polynomials that were considered to give an adequate amount of in- formation to represent the main field. An empirical rule given by Bullard [1967] relates the number of coefficien- ts of a global SH model with the minimum wavelength that a global model can represent and gives the details represented by the regional model. According to Bul- lard’s rule such relation brings to a minimum wavelen- gth of about 800 km for the second order degree poly- nomial representation for the area here considered. In Figure 3 all comparisons have shown that the use of a li- mited number of stations (that were afterwards used for the reduced network) had not introduced sensibly larger differences with respect to the use of the full set of avai- lable repeat stations data when compared to IGRF. References Bullard, E.C. (1967): The removal of trend from ma- gnetic surveys, Earth Planet. Sci. Lett., 2, 293-300. Coticchia, A., A. De Santis, A. Di Ponzio, G.Domi- nici, A.Meloni, M. Pierozzi, M. Sperti, (2001) La Rete Magnetica italiana e la Carta Magnetica d’I- talia al 2000.0, Bollettino di Geodesia e Scienze affini, 4, 261-291. Cafarella, L., A. De Santis and A. Meloni (1992a): Secu- lar variation from historical geomagnetic field me- asurements, Phys. Earth Planet. Inter., 73, 206-221. Cafarella, L., A. De Santis and A. Meloni (1992b): Il Ca- talogo geomagnetico storico italiano, Pubblicazio- ni dell’Istituto Nazionale di Geofisica vol. 7, 160 pp. De Santis, A., M. Chiappini, G. Dominici e A. Meloni (1997): Regional geomagnetic field modelling: the contribution of the Istituto Nazionale di Geofisica, Annali di Geofisica, XL, 5, 1161-1169. De Santis, A., L. Gaya-Piqué, G. Dominici, A. Meloni, J.M. Torta, and R. Tozzi (2003): ITalian Geomagne- tic Reference Field ITGRF: update for 2000 and se- cular variation model up to 2005 by autoregressive forecasting, Annals of Geophysics, 46, 3, 491-500. Dominici, G., A. Meloni, M. Miconi, M. Pierozzi and M. Sperti (2007): La Rete Magnetica Italiana e la Carta Magnetica d’Italia al 2005.0, Bollettino di Ge- odesia e Scienze affini,1,1-47. Dominici, G., A. Meloni, A. Di Ponzio, M. Miconi (2012): Italian Magnetic Network and magnetic re- ference fields at 2010.0, Annals of Geophysics, 55, 6, doi: 10.4401/ag-5411. Haines, G.V. (1990): Regional magnetic field model- ling: a review, J. Geomag. Geoel.,42, 1001-1018. IAGA (International Association of Geomagnetism and Aeronomy) Working Group V-MOD (2010): International Geomagnetic Reference Field: the eleventh generation, Geophysical Journal Interna- tional, 183, 1216-1230. Korte, M. and V. Lesur (2012): Repeat station data compared to a global geomagnetic field model, An- nals of Geophysics, 55, 1101-1111. Lanza, R., and A. Meloni (2006): The Earth’s Magnetism: An Introduction for Geologists, Springer, 290 pp. Mandea, M. (2004): Reseau magnétique de répétition de la France 2002, Bureau Central de Magnetisme terrestre, Paris, Bullettin n. 20. Meloni, A., O. Battelli, G. Dominici, S. Arca and A. Mar- chetta (1988): Italian magnetic network at 1985.0, Bollettino di Geodesia e Scienze Affini, 4, 339-350. Meloni, A., Battelli, O., De Santis, A. and Dominici, G. (1994): The 1990.0 magnetic repeat station survey and normal reference fields for Italy, Annali di Geo- fisica, Vol. XXXVII, N. 5, 949-967. Molina, F., E. Armando, R. Balia, O. Battelli, E. Bozzo, G. Budetta, G. Caneva, M. Ciminale, N. De Floren- tiis, A. De Santis, G. Dominici, M. Donnaloia, A. Elena, V. Iliceto, R. Lanza, M. Loddo, A. Meloni, E. Pinna, G. Santarato, R. Zambrano (1985): Geoma- gnetic Survey of Italy. Repeat Stations Network and Magnetic Maps: a Short Report, Annales Geophysi- cae, 3 (3), 365-368. Newitt, L.R., Barton C.E. and J. Bitterly (1996): Guide for magnetic repeat station surveys, IAGA Working Group, 115. Shepard, D. (1968): A two-dimensional interpolation function for irregularly-spaced data, In: Proceedin- gs 23rd National Conference ACM, ACM, 517-524. Talamo R (1975): Le carte magnetiche d’Italia delle isodinamiche della H e delle isogone dell’Istituto Militare e loro aggiornamento al 1973.0, Bollettino di Geodesia e Scienze Affini, Vol. XXXIV, n.1,p.34. ___________ *Corresponding author: Guido Dominici Istituto Nazionale di Geofisica e Vulcanologia, Roma, Italy; tel. +39-6-51860312 - fax +39-6-51860397 email: guido.dominici@ingv.it 2017 by Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved 6