Accelerations in the Local Magnetic Field on the Adriatic Tectonic Microplate ()
1. Introduction
In 2014 the observatory was built on the edge of the village of Sv. Peter above River Dragonja in the Slovenian part of Istria (45.459˚N, 13.685˚E, H = 196 m). After the international organization IAGA (International Association of Geomagnetism and Aeronomy), it has international code PIA (Piran, Slovenia) (hereinafter: the Observatory). We have been systematically looking for a suitable place to set up such an observatory since 2007 [1]. The method of its construction was mainly influenced by the presence of atmospheric electrical discharges [2], as well as the choice of its measuring instruments [3], which was specially adapted to the conditions in this area. The Observatory is located in the southwestern part of Slovenia, in the area with the highest annual density of lightning in Europe [4]. In addition to atmospheric discharges, it is also affected by vandalism, forest fires and floods due to surface waters. It is situated in a Natura 2000 area [5] and is not protected by a fence. This direct contact with nature allows us to monitor changes therein.
The setting up of the Observatory was most influenced by the general state of affairs in the society for which such an observatory is intended. International recommendations [6] and our experiences with the first geomagnetic observatory, SNV Observatory at Sinji vrh above Ajdovščina, were of great help in finding another additional site in Slovenia. The site of the SNV Observatory was 52.6 km from the PIA Observatory with the azimuth of 21.5 angular degrees.
On 20 May 2012, an unexpected change in the local magnetic field was measured at Sinji vrh. At that time we were just testing the magnetometer, so the measurement results have not been converted into real units of measurement and with a time delay of two hours. A detailed analysis of this event showed that the registration took place during a strong earthquake with magnitude ML = 5.1 (Richter earthquake scale) and with the epicenter near the city of Bologna in Italy. The change in all three components of the magnetic field (Figure 1) was very similar to the record of the seismograph [7]. The epicenter was 260 km from the measurement site with the azimuth 216.7 angular degrees. The registration itself was done with a three-axis fluxgate magnetometer that was not made for seismic zones because its sensor was suspended. That is why we chose a different type of magnetometer for the Observatory in Sv. Peter. The Slovenian part of Istria is built of flysch. The layer of hard limestone translates into seismic waves very well. For this reason, the sensor of the magnetometer was placed on a measuring pillar in a soft layer of sandstone, which is not a standard solution [8]. For the expected maximum amplitudes of seismic waves in the Observatory area, such a solution proved to be the right one. Earthquakes do not move the magnetometer sensor and therefore the measuring instrument for measuring changes of the magnetic field does not act as a seismograph (Figure 1). Therefore, even during an earthquake, it actually measures only changes in the magnetic field [9].
2. Greater Level of Noise in the Local Magnetic Field
The geomagnetic field on the territory of Slovenia contains a larger amount of noise than on the territory of neighboring countries. In July 2015 the SNV Observatory was still in function. At that time the nearest geomagnetic observatory with a long tradition was the TYH Observatory (Tihany, Hungary) [10]. A comparison of the measured values in last three hours 3 July 2015 shows that the noise level at the PIA Observatory is 38% higher than the level at the TYH Observatory. This level is increasing from the Adriatic coast towards the interior of Slovenia, so that at the SNV Observatory it is already 20% higher than at the PIA Observatory (Figure 2).
Figure 1. Measuring of the change in the local magnetic field at Sinji vrh above Ajdovščina during earthquake in Bologna (Italy) on 20 May 2012. Graf presents the variations of all three orthogonal components X (north), Y (east) and Z (nadir) of the geomagnetic field.
Figure 2. Magnetograms of the geomagnetic field density F [nT], drawn on the basis of measurements on 3 July 2015, at the SNV (Sinji vrh, Slovenia), PIA (Piran, Slovenia) and TYH (Tihany, Hungary) Observatories.
The geomagnetic field on the territory of Slovenia contains noise in component Z, which in the northern hemisphere acts vertically downward and in component X, which acts in the direction of the geographical pole (south-north) [11]. As this noise is a special feature of the Slovenian territory, we first had to verify the low noise level of the initial measuring instrument—an Oberhauser proton magnetometer [12] and then make comparative measurements at five locations with approximately the same distance along the same longitude from the location of the Observatory near the Adriatic coast to the middle of the country. In doing so, we found that the calmest area in Slovenia is on the edge of the Slovenian part of Istria. This was not shown by the GIS (geographic information system) analysis [6] or by the preliminary measurements [13].
The initial study of atmospheric discharges [3] was followed by the study of the influence of transitions of weather fronts on the local magnetic field [14] and artificial sources of noise in the magnetic field [15] [16]. The influence of the Moon could be studied in more detail only during the solar minimum at the transition from the 24th to the 25th solar cycle [17]. Knowing those influences allows us to exclude them from further processing of measurement data of the change of the local magnetic field.
Any relative changes in the energy density in the Earth’s magnetic field due to changes in tectonic stresses are very small—they are part of the noise. They can be evaluated on the basis of changes in the energy of the Earth’s magnetic field between individual days [18]. The energy density in the magnetic field is determined by Equation (1), which also applies to the energy density in the geomagnetic field wGM.
(1)
In Equation (1) are:
- W: energy in the magnetic field;
- V: volume;
- μ0: magnetic permeability of free space;
- B: magnetic flux density.
Two to three days before an earthquake occurs, the energy density in the local magnetic field increases greatly, which could be used for a short-term earthquake prediction. But an increase in this density is not yet a serious indicator of an earthquake. When an earthquake does occur, the density of the earth’s magnetic field is previously increased, but with an increased density, it is not necessary for an earthquake to occur.
3. Geomagnetic Jerks
Mathematical models of the change in the Earth’s magnetic field on its surface rely on secular changes in the magnetic field, which are changes in its components over long periods of time. In shorter intervals, however, these changes are uneven and therefore cannot be predicted. The sources of all changes in the Earth’s magnetic field are either inside of the Earth or in space. For geomagnetic jerks, the interior of the Earth is recognized as the place of their origin [19] [20]. Recent theoretical research suggests that they result from irregularities in the flow of the Earth’s liquid core [21] [22]. Hydrodynamic waves in the form of torsional oscillations are formed in its outer core. They are caused by bubbles that travel slowly towards the outer edge of the area of convection currents, interfering with their normal flow [23]. The consequences of this uneven flow in the outer core are various geophysical phenomena: small deviations of the Earth’s rotational axis with respect to land surface (Chandler wobble) and change in its speed with associated change in its climate [19]. There is also the possibility that geomagnetic jerks are the cause of stronger earthquakes and vice versa [24] [25].
Periods of secular changes in which the changes are constant are interconnected by geomagnetic jerks. In the graph of the secular changes of the geomagnetic field geomagnetic jerks are noticed as sudden changes in the shape of the letter V. They were first described in the 1970s [26]. Some of them occur only in a limited area of the Earth’s crust, such as individual continents or only parts of them, and occur every few months or years. But there are also those that are registered everywhere on the globe and happen unannounced every ten years on average.
In general, impulses are a sudden change in the form of an increase or decrease in a more or less uniform change in the observed phenomenon. Geomagnetic impulses or jerks are sudden changes that defined in the value of at least the first derivation of the time in the secular change of the Earth’s magnetic field [27]. Each of the components of the Earth’s magnetic field can be described in a second-degree polynomial Equation (2). In this polynomial equation geomagnetic jerks are recognized as a change in its coefficients. The closer these changes approach the shape of the letter V, the more the first derivation of the time is equal to the step and the second derivation of the time is equal to the pulse.
(2)
Geomagnetic jerks of a global nature have been investigated using both frequency analysis [28] and wavelet analysis [29]. For the individual observatories the initial research was conducted on the bases of annual mean values of measurements, from which the temperature drift of electronic measuring instruments is eliminated. Research of the geomagnetic jerks has improved the understanding of the mechanisms that take place inside the Earth. Measurements of the geomagnetic field with the help of state-of-the-art purpose build satellites have improved this understanding [30].
4. Accelerations in the Local Magnetic Field
The graph of mean daily values of the magnetic field, determined directly from measurements of magnetometers on the Observatory, shows that the strength of the local magnetic field F(t) [nT] in 2020 was constantly increasing (Figure 3). A larger number of sudden changes are also visible. The acceleration of the change in magnetic field strength is represented by the second derivative d2F(t) [nT/day2]. The value of this derivative in the whole of 2020 exceeded the value of 6.0 nT/day2 24 times, of which twice 12.0 nT/day2: the 34th day (March 22, 2020) and the 340th day (December 28, 2020).
The MSA (Mean Square Amplitude) algorithm was used for spectral analysis of FFT (Fast Fourier Transform) above the daily mean values of the local magnetic
Figure 3. Graph of daily mean values of the local magnetic field F(t) and its second derivative d2F(t) in Slovenia in 2020.
field in 2020 measured on the Observatory. The properties of this selected algorithm for conversion from time space to frequency space must be taken into account when analyzing the results obtained. After contributing to the spectral power density, the periods are divided into two consecutive groups (Figure 4). They are separated by a period of 14.07692 days, which is about half of the lunar month, which lasts an average of 29.53059 days. Equation (1) for the energy density in a magnetic field also proves that these two groups exist if the polynomial Equation (2) is replaced in it instead of the magnetic field density.
In 2020, seismographs of the national network of seismic observatories of the Republic of Slovenia registered 555 local earthquakes with a magnitude of ML ≥ 1 or an average of 1.5 earthquakes/day [31]. Local earthquakes are all earthquakes that have their epicenter on the territory of Slovenia (20,271 km2) or in its immediate vicinity. An earthquake occurs in its immediate vicinity if its epicenter is not more than 50 km from the nearest populated area in Slovenia.
The year 2020 was marked by two major earthquakes. The first occurred on the 82nd day (March 22, 2020) in Zagreb with the largest magnitude ML = 5.5 [32], and the second one on the 363rd day (December 28, 2020) near Petrinja with the largest magnitude ML = 6.2 [33]. Both earthquakes were on the territory of the neighboring Republic of Croatia. On December 28, 2020 there were no local earthquakes in Slovenia. For this reason the earthquake in Petrinja is only marked with m = 70 earthquakes/day (Figure 5).
5. Discussion
For a period of a few days, the power contribution of all the periods that make up the year-round change in the mean daily values of the local magnetic field is negligibly small. The transition to the part with larger amplitudes is at a period
Figure 4. Sequence of periods of less than 100 days or 0.27 years obtained by spectral analysis of FFT above the mean daily values of the local magnetic field F(t) [nT] in 2020. Due to a better presentation of the obtained results of conversion into frequency space, these are fitted together with the Basis spline (B-spline).
Figure 5. Graph of the second derivative d2F(t) of daily mean values of the local magnetic field F(t) and the number of local earthquakes in Slovenia in 2020.
of the half of the lunar month (Figure 4). The period of 7.4 days or the mean length of each lunar phase is very poorly presented in the power spectrum. Further research would be needed to confirm the link between the gravitational influence of the Moon and geophysical phenomena associated with accelerations in the local magnetic field. The sequence of period pairs indicates the pulsation generated by frequencies longer than one year. Therefore, eventual further research for additional confirmation the findings presented in this article, should cover a period of at least a few years.
In 2020, no major earthquake occurred on the territory of Slovenia, but there were two strong earthquakes on the territory of neighboring Croatia. The first one occurred 48 days after the appearance of the first large impulses in the geomagnetic field with acceleration greater than 12.0 nT/day2 was earthquake with local characteristic. The second one, 23 days after the second such large impulse in 2020, was earthquake with regional characteristic. Therefore, we find that all the geophysical phenomena related to accelerations in local magnetic field should be observed on a wider territory than that of Slovenia.
6. Conclusion
Reliable earthquake prediction would be desirable, but it will only be possible if we first develop appropriate algorithm and test it on a limited territory. The algorithm should be based on different measurements, which arise from different physical principles. Among them seismological measurements would be basic. The said territory should be a seismological and tectonically limited area without volcanic activity such as the Adriatic tectonic microplate. Only the algorithm with great reliability could eventually be applied to outer areas of the world.
Acknowledgements
The authors would like to thank all those who, as volunteers, helped set up and maintain the PIA Observatory. We would also like to thank all those who help write this article on behalf of the professional journal OJER.