Next Article in Journal
Synthesis and Biological Activity of 1,3,4-Oxadiazoles Used in Medicine and Agriculture
Next Article in Special Issue
The Variation of Counter-Electrojet Current at the Southeast Asian Sector during Different Solar Activity Levels
Previous Article in Journal
Prospects of Structural Similarity Index for Medical Image Analysis
Previous Article in Special Issue
Evaluating the Effect of Noise from Traffic on HYB Magnetic Observatory Data during COVID-19 Lockdown
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring of Geomagnetic and Telluric Field Disturbances in the Russian Arctic

by
Olga V. Kozyreva
1,2,
Vyacheslav A. Pilipenko
1,2,*,
Elena E. Marshalko
3,
Elena Yu. Sokolova
1 and
Mikhail N. Dobrovolsky
2
1
Institute of Physics of the Earth, 123242 Moscow, Russia
2
Geophysical Center, 119296 Moscow, Russia
3
Finnish Meteorological Institute, FI-00560 Helsinki, Finland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(8), 3755; https://doi.org/10.3390/app12083755
Submission received: 19 February 2022 / Revised: 25 March 2022 / Accepted: 6 April 2022 / Published: 8 April 2022

Abstract

:
The influence of space factors on technological systems in the Arctic (power transmission lines, oil/gas pipelines) has become critically important. To examine in depth these effects, an archive of digital 1 min data from Soviet/Russian magnetic stations deployed along the Arctic coast was created, starting from 1983 to the present. All data from various sources were converted to daily files in standard IAGA-2002 format and supplemented with quick-look magnetograms. Some of these data are included already in the existing world magnetic field databases, but not all. Examples of disturbances known to excite intense geomagnetically induced currents in power transmission lines were presented: irregular Pi3 pulsations and magnetic perturbation events. The database was augmented with the global 3D model of the Earth’s conductivity structure. The given example showed how the combined usage of the geomagnetic field database and the conductivity model enables one to synthesize the geoelectric field response to geomagnetic variations, and to assess the distortions of the pipeline-soil potential. To determine regions most susceptible to geomagnetic hazard, a map with normalized telluric fields was created for a uniform sinusoidally varying magnetic disturbance. This map showed that the largest electrotelluric potentials and field are induced in regions with a high resistivity (e.g., Kola Peninsula and Ural Mountains). This database can be also a useful support for space missions in the magnetosphere. The database is publicly available on the anonymous FTP site.

1. Possible Risks to Pipelines from Geomagnetic Variations

Space weather is a serious threat to the regular functioning of satellite and ground-based technological systems. One of the most significant space weather hazards is an electric geomagnetically induced current (GIC) in grounded low-resistance industrial structures—power transmission lines and pipelines. GICs are excited by abrupt changes in the geomagnetic field, i.e., large dB/dt. The most severe geomagnetic disturbances and GICs induced by them in technological systems are observed at auroral latitudes.
The actual driver of GIC is the variations in the geoelectric (telluric) field. Modern energy and pipeline systems are large networks with a complex topology that cover vast areas of the Earth’s surface, the local geoelectric properties of which can vary by up to five orders of magnitude. The electric conductivity of the Earth’s interior ranges from ~10−5–10−4 S/m in crystalline silicate rocks up to ~3 S/m in the ocean. However, it is not easy to obtain direct information about geoelectric fields. While variations in the geomagnetic field are monitored by a worldwide network of magnetometers, there are very few regular observations of telluric electric fields. Although the calculation of telluric fields and currents from geomagnetic data is feasible, it requires not only a sufficiently dense network of magnetometers, but also the availability of adequate information on the distribution of electrical conductivity in the Earth.
The information about the geoelectric structure of the Earth’s interior is scarce. For North America, this information was provided by a continental-scale magnetotelluric (MT) sounding project—EarthScope [1]. Within the framework of this project, an MT survey was carried out on a grid with a distance of 70–100 km between temporary sites, and the corresponding impedance tensors were estimated using these MT data. The impedances of the EarthScope project were further used in [2] for the estimation of the geoelectric field response to a uniform sinusoidally varying geomagnetic field. In a similar way, an array of MT impedance measurements in Alberta (Canada) was used to map the spatial structure of the geoelectric field resulting from a magnetic storm [3].
After the Quebec energy blackout of 1989, it was found that space weather not only creates a giant natural plasma laboratory, but also poses a serious threat to the uninterrupted functioning of space- and ground-based technological systems. Nowadays, research is being actively carried out worldwide on the influence of GIC on ground-based technological systems and possible measures to reduce the negative consequences. A negative impact of GICs on power transmission lines and high-voltage transformers is a thoroughly studied phenomenon (see review [4]). At the same time, there is very little scientific information on the impact of geomagnetic disturbances and associated telluric field variations on pipelines. The development of GIC in pipelines results in an increased corrosion at grounding points or at points with insulation damage, as well as disruption in the operation of cathodic protection. Commonly, the cathodic protection against electrocorrosion feeds a negative potential of ~1–2 V between the pipe and soil. Perturbations in the geomagnetic field induce fluctuations in the pipe-to-soil potential (PSP), which take the pipeline voltage beyond the safe range of electrocorrosion protection. During geomagnetic storm/substorms, changes in PSP can be large enough to expose pipeline sections for extended periods of time. For example, on a pipeline in Quebec (Canada), corrosion damage of a pipeline section developed after 5 years of operation only instead of the expected 20–30 years [5]. The extensive studies of the effect of space weather on pipelines were carried out on the Finnish gas pipeline [6,7]. The number of cases with significant (>10 A) currents in the pipeline in southern Finland correlated well with the number of sunspots, which confirmed the direct relationship of the GIC with the solar/geomagnetic activity [8]. In the trans-Alaska pipeline during strong magnetic disturbances, current surges up to 200 A occurred [9]. In Northern Alberta (Canada), failures in the cathodic protection associated with geomagnetic disturbances were identified [10]. In a pipeline in Northern Norway, it was found that the PSP fluctuates with an amplitude of ~5 V during magnetic disturbances [11].
Strong geomagnetic disturbances can cause quite significant effects in technological systems at mid-latitudes as well. In Australia and New Zealand, during severe storms, PSP fluctuations on pipelines for a long time were several times higher than the established limits [12,13]. In the pipeline in northern Bavaria, the peak current values upon a moderate geomagnetic activity reached ~12 A, and the PSP variations were ~3 V during the magnetic storm main phase [14]. The largest PSP variations associated with geomagnetic fluctuations were observed in the Bolivia-Brazil gas pipeline buried 0.5 m under the ground at the ends of the pipeline [15].
The Arctic zone of the Russian Federation is a promising testing ground for the studies of GIC influence on ground-based technological systems. In this territory, the GIC monitoring system along the power transmission line at the Kola Peninsula has permanently been operating since 2011 [16]. It is the only system that provides the GIC monitoring data for academic research. In this area, especially at the Yamal Peninsula and the Arctic shelf, the centers of the natural gas production and transportation were formed.
In addition to practical applications, ground-based magnetometers form the basis for the studies of electromagnetic processes in near-Earth space. The study of the mechanisms of solar–terrestrial relationships has also led to the creation of a new scientific discipline—the “space climate”. For the space climate studies, it is very important to have long-term time series of geomagnetic observations. The Arctic territory of the Russian Federation, which covers approximately 1/3 of the auroral oval, can be a viable support to satellite missions in the ionosphere and magnetosphere.
Under the harsh Arctic conditions, Russian scientific organizations support long-term observations (the list of institutes maintaining regular magnetometer observations at high latitudes is given at the end of the paper). IKFIA conducts complex observations at a network of stations in Eastern Siberia and the Far East. AARI has deployed magnetic stations along the Arctic coast up to Chukotka. PGI has a complex of observations on the Kola Peninsula and Svalbard. New magnetic stations were installed by GC [17]. IZMIRAN deployed several stations on the Yamal Peninsula. However, the data of these observations are scattered among various sources and stored in different formats. Although the monitoring of the geomagnetic field in the Russian Arctic zone is extremely important for space and applied science, the availability of data from Russian stations to the world research community is still limited.
This paper describes a new database, which includes an archive of magnetovariational observations in the Russian Arctic and gives examples of its use for scientific and applied problems. Particular attention is paid to the territory of Yamal and Western Siberia, where the active development of the infrastructure is in progress now. These experimental data can serve as a basis for studying the features of geomagnetically induced electric fields. The database is augmented with the global 3D model of the Earth’s conductivity structure. The information on the geoelectric structure makes it possible to obtain estimates of extreme GIC values in infrastructure facilities and assess the possible increase in electrocorrosion of pipelines and the reduction in their service life.

2. Description of the Russian Arctic Database

The presented database contains the results of observations of the geomagnetic field variations with a 1 min sampling period at stations in the Russian Arctic. The stations’ geographic and corrected geomagnetic coordinates are given in Table 1, and the map with the stations’ locations is shown in Figure 1.
The data of some stations are also available on the SuperMAG portal, some are part of the INTERMAGNET network, and some are included in the IMAGE and MAGDAS networks. Part of the data were taken from the Common Data Center “Analytical Center for Geomagnetic Data”. However, data from several stations, for example, on Yamal, are not included in these archives. Naturally, in each of the above archives, data are stored in their own format.
The measurements of the absolute values of the geomagnetic field are carried out only at stations included in the INTERMAGNET network. Other stations are variational, so only deviations from a reference level have a physical meaning in their records. The data quality and sensitivity are the same as in the data portals mentioned above. All the collected data have been converted into the standard ASCII format IAGA-2002. The names of the files have the format STAYYYYMMDDvmin.min, where STA is the station code from Table 1, YYYY—year, MM—month, DD—day. Such file names make it easy to select the required intervals for analysis. The data files are sorted by year and stored in a compressed format to speed up the download process. The data are converted into the geomagnetic coordinate system (H, E, Z): H axis is directed to the geomagnetic pole, E axis is directed to the geomagnetic east, and Z axis is directed vertically downward. The magnetograms of all stations were visually examined and evidently defective files were discarded. However, the remaining files may still contain level jumps, data gaps, and intermittent noise, because their correction requires a special analysis.
To quickly browse the data, check for the data availability for specific time intervals, and select events, the quick-look magnetograms are drawn. For a general characterization of conditions in interplanetary space from the OMNI database, the same graph shows variations in the solar wind plasma density Np, solar wind velocity V, and the vertical component of the interplanetary magnetic field Bz.
Examples of the constructed magnetogram plots are demonstrated below. The presented events are important for the problem of GIC excitation because they are known to be associated with intense GIC bursts, but their physical nature has not been completely clarified.

3. Examples of Geomagnetic Variations Causing Significant GIC Bursts

The largest magnetic disturbances on the Earth’s surface are caused by an extended large-scale auroral electrojet, intensified during substorms [18]. However, the most intense GIC bursts are produced by small-scale ionospheric current structures [19]. The energy of such impulsive or quasiperiodic disturbances is much lower than the energy of magnetospheric storms or substorms, but rapidly changing fields of such disturbances can cause GIC bursts of significant magnitude [20]. Examples of such impulsive and quasiperiodic disturbances of geomagnetic field are given below based on the data from the Russian Arctic database.
After the expansion phase of a substorm, large-scale magnetic variations are often superimposed by irregular quasiperiodic Pi3 pulsations with characteristic time scales of 3–15 min. Figure 2 shows quick-look graphs for the substorms on 9 February 2019 (left-hand panel) and 2 March 2019 (right-hand panel) from the Arctic magnetometer database. These Pi3 pulsations are not harmonic oscillations such as Pc5 waves, but a series of nonlinear magnetic pulses with steep fronts. Ps6 pulsations in the same frequency range are morphologically similar [21]. Due to steep fronts of the Pi3/Ps6 pulses, the associated field variability is high, up to dB/dt~20 nT/s. In many events with extremely large GICs, the currents reached their maximum value (up to hundreds of Amperes) during the subsequent series of Pi3/Ps6 pulses [22,23].
When considering the impacts of space weather, it is usually assumed that disturbed geomagnetic fields are spatially uniform throughout the area of a technological system. However, often against the background of a general field enhancement during substorms, isolated magnetic perturbation events (MPEs) are observed in the magnetic (>100 nT) and telluric (>1 V/km) fields [24]. The MPEs are spatially localized, so the peak at one station can be several times larger than the regional value of the field on average (at distances of about 500 km). Figure 3 shows examples of such isolated impulsive disturbances during moderate magnetic activity on 9 September 2019 (at ~16 UT) and 17 December 2019 (at ~05 UT). These magnetic pulses are usually associated with a local intensification of the auroral luminosity. Such pulses may induce sporadic intense bursts of GIC. The physical processes that determine the generation of the MPEs are poorly understood.
There are other geomagnetic disturbances that are hazardous to power transmission and pipeline networks. Due to the high variability of the geomagnetic field of Pc5 pulsations (3–10 min), the magnitude of associated periodic GICs in power lines may increase and exceed 10 A [25]. The long-lasting Pc5 wave activity can be more dangerous than short-term bursts of GICs during substorm onsets due to cumulative effects—corrosion and premature aging of high-voltage transformers.

4. Modeling Disturbances of the Geoelectric Field

Using information on the geoelectric structure of the underlaying crust, one may reproduce the dynamics of the telluric response for real events from geomagnetic observations, even in the absence of synchronous electric field records. The construction of a detailed global model of the Earth’s electrical conductivity is a task for the future; therefore, nowadays, calculations have to use some approximate scheme [26]. A comparison of various methods showed that with a good accuracy for calculation of the telluric field, one can use the impedance relation in the plane wave approximation. This approximation, also named the strong skin-effect condition, is valid when the horizontal scale of the disturbance is much larger than the skin-depth d = 2 / μ 0 ω σ of the electromagnetic field with cyclic frequency ω in the conducting Earth with electrical conductivity σ (Wait–Price criterion) [27,28]. Under this condition, the actual spatial structure of geomagnetic disturbance does not influence the impedance relationship, so the primary disturbance should not necessarily be a “plane vertically incident” wave. For practically important calculations of GIC, the situation is greatly simplified by the fact that integral estimates of the potential difference throughout an extended system (at least several hundred kilometers) are needed, and therefore the necessary estimates can be performed with a sufficient accuracy with a relatively rare network of magnetometers and a rough model of conductivity.
The geoelectric field in the plane wave approximation is coupled with geomagnetic variations through the surface MT impedance Z(ω). In media with a low conductivity, the negative effects of strong magnetic disturbances are most probable, because the excited telluric fields and potential differences turn out to be larger, so the induced currents mainly flow through conductive elements of industrial networks. On the Earth’s surface, the impedance relation holds between the spectral amplitudes of the horizontal electric E = {Ex, Ey} and magnetic B = {X, Y} fields as follows: E = μ o 1 Z B , where Z is the impedance matrix. In a simplified approach, the use of 1D impedance Z(ω) in the impedance boundary condition may be sufficient, namely
( E x ( ω ) E y ( ω ) ) = μ 0 1 ( 0 Z ( ω ) Z ( ω ) 0 ) ( X ( ω ) Y ( ω ) )
For uniform conductivity of the Earth, the impedance Z ( ω ) = i ω μ 0 ρ (where ρ = 1/σ is the electric resistivity) depends on the period of the incident wave as a power law Z ( T ) T 1 / 2 . In the vicinity of coastal areas with a sharp gradient of conductivity, the use of 3D impedances may be necessary to account for horizontal inhomogeneity [29]. The telluric field enhancement in such regions may be as high as several times. The voltage on a virtual transmission line modeled using a piecewise 1-D conductivity profile was estimated to differ by 10–20% (depending on orientation relative to the geoelectric field polarization) compared with modeling that accounted for 3-D conductivity structure [3].
Here, using the recorded time series of B(t) = {H, E} in geomagnetic coordinates and frequency domain MT impedance Z(ω), we synthesized the temporal variations in the electric field E(t) = {EH, EE} by a two-step procedure. First, a spectral transform of the recorded time series was used in the impedance condition (1), and then, the inverse Fourier transform was applied to the spectrum of the telluric field. The selected example corresponds to the substorm with Pi3 activity recorded at a local latitudinal profile of stations at Yamal. The magnetograms from stations along the latitudinal profile 150° during the event on 6 January 2019 are shown in Figure 4. The peak magnetic disturbance in the H component is reached at BEY station ΔH~350 nT.
To calculate the geoelectric fields that could be excited by this geomagnetic disturbance, we used information on the impedances of the Earth’s crust calculated on the basis of the model from [30]. This global 3D model was compiled using the results of many regional MT sounding projects, as well as current ideas about the structure of the outer solid shells of the Earth and their electrical properties. The model describes the structure of the electrical conductivity of the Earth in the depth range of 0–100 km, including seawater, crust, and partly the upper mantle. This model was augmented by the data of the upper 10 km of the SMAP conductivity model of Fennoscandia [31], which was not originally included in [30]. To describe the conductivity distribution at depths larger than 100 km, the global 1D conductivity profile from [32] was used. Figure 5 (left-hand panel) demonstrates the vertical conductivity profiles σ(z) at magnetometer stations’ locations at Yamal, whereas Figure 5 (right-hand panel) shows the corresponding frequency-dependent impedance Z(f). At station SKD, the dependence is close to a power law, which corresponds to the homogeneous conductivity. At other stations, the impedances at f > 0.5 mHz deviates noticeably from a power law dependence.
Figure 6 demonstrates the results of the synthesis of electric field time series during the substorm event on 6 January 2019 for station BEY. The superposed Pi3 pulsations on a magnetic bay are most evident in the eastward magnetic component E. These pulsations provide a high variability of the geomagnetic field, dH/dt ~ 5 nT/s, dE/dt ~ 6.5 nT/s. The synthetic electric field is also highly variable, though its spectral content is evidently different from dB/dt variations. The highest peak-to-peak variations in the telluric field are EH~EE ~ 0.4 V/km. The occurrence of electric field peaks does not coincide with the moment of maximum magnetic field disturbance (cf. Figure 4).
The spectra of magnetic and telluric field variations at station BEY during the 6 January 2019 event (shown in Figure 6) are presented in Figure 7. The spectral amplitude of magnetic field variations decays with frequency and has a spectral peak of both components in the band 2–3 mHz caused by Pi3 pulsations. The electric field has similar behavior with spectral peaks in the same band. In contrast, the spectra of magnetic field variability dH/dt and dE/dt have a weakly decaying dependence on frequency with only a slightly evident spectral peak at the Pi3 frequency. Thus, this comparison shows that in the telluric field variations, the high-frequency part of the spectrum is substantially suppressed as compared with dB/dt fluctuations [33]. Indeed, the spectral power of horizontal magnetic components drops ~12.5 times from 1 mHz to 6 mHz, whereas that of the electric field drops ~8 times, and that of dB/dt drops ~3.5 times only.
The considered event and other observations [22,23] showed that the most intense GICs in power transmission lines are excited by geomagnetic fluctuations in the Pi3/Ps6 band. The formation of the highlighted frequency band in the GIC response most vulnerable to geomagnetic variations can be understood from the following simple consideration, illustrated in Figure 8. GICs can be imagined as currents induced in the electric contour formed by a power transmission line, grounding, and conductive crust. Very approximately, the induced current in this contour according to Faraday’s law is proportional to temporal variations in the geomagnetic field, J ~ dB/dt. That is why the magnetic field variability dB/dt is a commonly used measure of possible GIC, because information on telluric field E is not easily available. On the other hand, GIC can be imagined as a “discharge current” via a grounded technological system between distant points with different telluric field potentials. These approaches do not contradict each other, because in both cases, J ( f ) f α B ( f ) , where α = 1 for the dB/dt scenario and α = 1/2 for the telluric field in homogeneous Earth. The spectral power of magnetic variations B(f) typically decays with frequency f, whereas the power of both dB/dt variability and E-field fluctuations grows with frequency. The resulting GIC is a convolution of those two factors and must have a maximum at some intermediate frequencies. Dedicated studies have shown that this frequency band falls into the range of 3–15 min, that is, the Pi3/Ps6 band.

5. Map of Telluric Field Response to a Geomagnetic Disturbance over Northern Russia

For qualitative operative estimates of expected vulnerability of a technological system to geomagnetic variations, it is desirable to have a map with the telluric response amplitude to a given disturbance of the magnetic field. A map, similar to [2], but for the territory of the Russian Federation, was constructed here. The map is based on the geoelectric model [30]. The spatial distribution of the average conductivity in the upper 10 km layer of the Earth’s crust derived from this model for the territory of the Russian Federation is shown in Figure 9. The lowest conductivities, ~10−5–10−4 S/m, correspond to the Precambrian basement outcrops (Fennoscandian Shield) and to the exposed nuclei of the Paleozoic and Mesozoic orogenic structures: Kola Peninsula, Karelia, Urals, Novaya Zemlya, Taimyr, and Chukotka. In the Yamal region covered by thick Mesozoic sedimentary cover, the conductivity is generally quite high.
The spatial distribution of the electric field response to geomagnetic disturbance was calculated from the corresponding impedances using harmonic approximations of geomagnetic disturbances. The simplest approximation of the driving magnetic field was taken in the form of a plane wave with an amplitude of 1 nT and oscillating with a fixed frequency. The magnetic variations were assumed to be in the longitudinal direction, though for the considered 1-D algorithm of the telluric electric field calculation, the actual polarization is not important. The calculation of the electric field according to the available MT impedance was performed on a 0.25° × 0.25° grid using local vertical profiles of conductivity in each node from the model [30]. The calculated spectral amplitudes were multiplied by the spectral width to obtain the amplitude (in mV/km) of the excited telluric field. The examples given above indicate that the most intense GICs are provided by geomagnetic pulsations with time scales of about 100–1000 s, so these periods were chosen for modeling.
Figure 10 shows a map of the magnitude of synthetic geoelectric fields that could be excited by spatially homogeneous magnetic field variations with T = 100 s throughout the territory of the Russian Federation. According to this “telluric hazard” map, the highest values of geomagnetically induced fields up to EH ~ 5 mV/km (normalized to the exciting magnetic field) are typical for the northwest of the Russian Federation. In this region, the magnitudes of the telluric field for different points can differ 3–4 times under the same amplitudes of the magnetic disturbance. In addition, increased values should occur in the Urals. In the Yamal region, the telluric fields are weak, ~0.5 mV/km.
A similar map, but for the driving magnetic disturbance with T = 1000 s, is shown in Figure 11. It demonstrates that a lower-frequency variation excites a lower amplitude of telluric currents. The maps in Figure 9, Figure 10 and Figure 11 prove that the geological structure significantly affects the amplitude of geoelectric fields generated by magnetic disturbance. The constructed map of normalized amplitudes of geoelectric responses to external plane-wave geomagnetic excitation provides helpful information for system operators in regard to space weather disturbances.
The horizontal scale of even localized magnetospheric disturbances on the ground (such as magnetic impulses composing Pi3/Ps6 pulsations) is about 200–300 km. For the average resistivity of the crust in the Arctic region of ρ~3–50 Ohm·m, the strong skin-effect condition remains valid for periods under consideration. The impedance approach may cause doubt only for long-period disturbances with time scales 5–20 min above those of high-resistivity regions, such as the Kola Peninsula with ρ~104–105 Ohm·m. In this specific region, the full wave solution of the electromagnetic disturbance structure in the crust may be needed.
Should the PSP changes during strong substorms be taken into account when maintaining the cathodic potential to protect against electrocorrosion? This question may be answered using the constructed map. Different analytical and numerical methods have been elaborated to calculate the induction in a grounded conductor [34,35,36,37,38]. Those methods demand prior information on the telluric field intensity, polarization, and spatial structure. According to calculations in [35], currents in the south of Finland for an electric field of 1.0 V/km should cause a GIC of ~50 A in the pipeline and up to 25 A in its grounding, which, in the case of insulation failure, may well lead to a noticeable enhancement of corrosion. A very crude estimate using the results of the measured PSP in oil pipelines was carried out in [39]. They showed that the telluric field calculated in the approximation of uniform conductivity is linearly proportional to the measured PSP. The relation found shows that the potential created by telluric fields can exceed the maintained PSP~2 V at E > 0.5 V/km. According to the constructed map (Figure 10), applying this result to the Yamal region (the area of oil/gas productive Mesozoic sediments), a geomagnetic disturbance of ~500 nT with a time scale of ~100 s is required to cause a disturbance of PSP menacing the cathodic protection.

6. Conclusions

This paper describes the new freely available database of geomagnetic field monitoring in the Russian Arctic, which may be helpful for fundamental research of the near-Earth space electrodynamics and for practical applications. To assess possible risks associated with space weather phenomena for technological facilities, a thorough analysis of the geoelectric structure of the region is necessary, as inhomogeneities in the distribution of terrestrial electrical conductivity can significantly affect the pattern of GICs and PSPs, modifying their amplitudes and modulating the spectral composition. The constructed map of the distribution of the synthesized telluric fields indicates domains of the Russian Arctic with the most intense responses. The estimates obtained show that the horizontal geoelectric field, induced by harmonic oscillations B(t) with an amplitude of 1 nT and period of 100 s, varies greatly in magnitude in the range from <1 mV/km to ~120 mV/km. For the considered moderate substorm, the variation in the telluric field in the Yamal Peninsula exceeded 2 V/km. The influence of geomagnetic variations caused by magnetospheric disturbances should be taken into account when maintaining the cathodic protection system of pipelines, as the effect of GIC and PSP can manifest itself both directly during the development of disturbance as failure of the equipment and have a long-term cumulative nature as electrocorrosion or aging of industrial transformers [40].
The obtained estimates of the possible amplitudes of the geoelectric field can be used in calculating the extreme values of the induced GIC to optimize the cathodic protection of pipelines from the destructive effects of strong geomagnetic disturbances. The constructed map in combination with the Arctic archive of geomagnetic data allows one to quickly assess possible risks to regional infrastructure under various disturbances of the geomagnetic field. Using the proposed magnetic field database and the conductivity model, any researcher can calculate the disturbance of telluric fields in the Russian Arctic for any time period.
The entire database, graphs of magnetograms for quick-look viewing, and maps of telluric fields are freely available on the anonymous FTP site. The database is constantly being supplemented and improved, so the authors will appreciate any suggestions for its improvement and development.

7. Data Provision

Digital geomagnetic data were deployed and maintained by the following institutes:
  • AANII—Arctic and Antarctic Research Institute, Rosgidromet, Saint-Petersburg;
  • GC—Geophysical Center, Russian Academy of Sciences, Moscow;
  • IKIR—Institute of Cosmophysics and Radiowaves, Russian Academy of Sciences, Kamchatka;
  • IKFIA—Institute of Cosmophysical Research and Aeronomy, Siberian Branch of Russian Academy of Sciences, Yakutsk;
  • ISZF—Institute of Solar-Terrestrial Physics, Siberian Branch of Russian Academy of Sciences, Irkutsk;
  • IZMIRAN—Institute of Terrestrial Magnetism and Radiowaves, Russian Academy of Sciences, Moscow;
  • PGI—Polar Geophysical Institute, Kola Scientific Center of Russian Academy of Science; Apatity.
The data of the GIC monitoring system along the power transmission line at the Kola Peninsula are at (http://www.eurisgic.org, accessed on 18 February 2022).
The data of some stations from the database are available on the SuperMAG portal (https://supermag.jhuapl.edu, accessed on 18 February 2022), the INTERMAGNET network (https://www.intermagnet.org, accessed on 18 February 2022), and the IMAGE (https://space.fmi.fi/image) and MAGDAS (http://www.serc.kyushu-u.ac.jp/magdas, accessed on 18 February 2022) arrays.
Part of the data were taken from the Common Data Center “Analytical Center for Geomagnetic Data” (http://ckp.gcras.ru, accessed on 18 February 2022).
The parameters of interplanetary space are taken from the OMNI database (https://omniweb.gsfc.nasa.gov, accessed on 18 February 2022).
The description of the IAGA-2002 format is at (https://www.ngdc.noaa.gov/IAGA/vdat/IAGA2002/, accessed on 18 February 2022).

Author Contributions

Conceptualization and writing—V.A.P., data processing and analysis—O.V.K.; telluric field modeling and mapping—E.E.M.; review and editing—E.Y.S.; database maintenance—M.N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Russian Science Foundation grant 21-77-30010 (OK, VP, MD), and by project 314670 of the Academy of Finland (EM).

Informed Consent Statement

Not applicable.

Data Availability Statement

The entire database is freely available on the anonymous FTP site ftp://door.gcras.ru/ftp_anonymous/ARCTICA_Rus accessed on 18 February 2022. The usage of the database must be referenced https://doi.org/10.2205/Rus-Arctic-1-min-DB.

Acknowledgments

We appreciate useful discussions with A. Zaitsev and help from A. Grudnev.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schultz, A. EMScope: A continental scale magnetotelluric observatory and data discovery resource. Data Sci. J. 2009, 8, IGY6–IGY20. [Google Scholar] [CrossRef]
  2. Bedrosian, P.A.; Love, J.J. Mapping geoelectric fields during magnetic storms: Synthetic analysis of empirical United States impedances. Geophys. Res. Lett. 2015, 42, 10160–10170. [Google Scholar] [CrossRef] [Green Version]
  3. Cordell, D.; Unsworth, M.J.; Lee, B.; Hanneson, C.; Milling, D.K.; Mann, I.R. Estimating the geoelectric field and electric power transmission line voltage during a geomagnetic storm in Alberta, Canada using measured magnetotelluric impedance data: The influence of three-dimensional electrical structures in the lithosphere. Space Weather 2021, 19, e2021SW002803. [Google Scholar] [CrossRef]
  4. Pilipenko, V.A. Space weather impact on ground-based technological systems. Sol. Terr. Phys. 2021, 7, 68–104. [Google Scholar] [CrossRef]
  5. Boteler, D.H.; Trichtchenko, L. Telluric influence on pipelines. In Oil and Gas Pipelines: Integrity and Safety Handbook; Revie, R.W., Ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015; Chapter 21; pp. 275–285. [Google Scholar] [CrossRef]
  6. Pulkkinen, A.; Viljanen, A.; Pajunpaa, K.; Pirjola, R. Recordings and occurrence of geomagnetically induced currents in the Finnish natural gas pipeline network. J. Appl. Geophys. 2001, 48, 219–231. [Google Scholar] [CrossRef]
  7. Viljanen, A.; Pulkkinen, A.; Pirjola, R.; Pajunpaa, K.; Posio, P.; Koistinenet, A. Recordings of geomagnetically induced currents and a nowcasting service of the Finnish natural gas pipeline. Space Weather 2006, 4, S10004. [Google Scholar] [CrossRef]
  8. Viljanen, A. Geomagneticalliy induced currents in the Finnish natural gas pipelines. Geophysica 1989, 25, 135–159. [Google Scholar]
  9. Campbell, W.H. Observation of electric currents in the Alaska oil pipeline resulting from auroral electrojet current sources. Geophys. J. R. Astron. Soc. 1980, 61, 437–449. [Google Scholar] [CrossRef]
  10. Shapka, R. Geomagnetic effects on modern pipeline systems. In Proceedings of the Solar-Terrestrial Predictions Workshop, Ottawa, ON, Canada, 18–22 May 1992; pp. 163–170. [Google Scholar]
  11. Henriksen, J.F.; Elvik, R.; Gransen, L. Telluric currents corrosion on buried pipelines. In Proceedings of the 8th Scandinavien Corrosion Congress, Helsinki, Finland, 15–23 August 1978; Volume II, pp. 167–176. [Google Scholar]
  12. Marshall, R.A.; Waters, C.L.; Sciffer, M.D. Spectral analysis of pipe-to-soil potentials with variations of the Earth’s magnetic field in the Australian region. Space Weather 2010, 8, S05002. [Google Scholar] [CrossRef]
  13. Ingham, M.; Rodger, C.J. Telluric field variations as drivers of variations in cathodic protection potential on a natural gas pipeline in New Zealand. Space Weather 2018, 16, 1396–1409. [Google Scholar] [CrossRef]
  14. Brasse, H.; Junge, A. The influence of geomagnetic variations on pipelines and an application for large-scale magnetotelluric depth sounding. J. Geophys. 1984, 55, 31–36. [Google Scholar]
  15. De Morales, J.F.; Paulina, I.; Alves, L.R.; Denardini, C.M. Evaluation of possible corrosion enhancement due to telluric currents: Case study of Bolivia-Brazil pipeline. Ann. Geophys. 2020, 38, 881–888. [Google Scholar] [CrossRef]
  16. Sakharov, Y.A.; Danilin, A.N.; Ostafiychuk, R.M. Registration of GIC in power systems of the Kola Peninsula. In Proceedings of the 7th Symposium on Electromagnetic Compatibility and Electromagnetic Ecology, St. Petersburg, Russia, 26–29 June 2007; pp. 291–293, ISBN 9781424412693. [Google Scholar]
  17. Gvishiani, A.D.; Soloviev, A.A.; Sidorov, R.V.; Krasnoperov, R.I.; Grudnev, A.A.; Kudin, D.V.; Karapetyan, J.K.; Simonyan, A.O. Successes of the organization of geomagnetic monitoring in Russia and the near abroad. Vestn. Otd. Nauk Zemle 2018, 10, NZ4001. [Google Scholar] [CrossRef]
  18. McPherron, R.L. Physical processes producing magnetospheric substorms and magnetic storms. In Geomagnetism; Jacobs, J., Ed.; Academic Press: London, UK, 1991; pp. 593–739. [Google Scholar]
  19. Belakhovsky, V.; Pilipenko, V.; Engebretson, M.; Sakharov, Y.A.; Selivanov, V. Impulsive disturbances of the geomagnetic field as a cause of induced currents of electric power lines. J. Space Weather Space Clim. 2019, 9, A18. [Google Scholar] [CrossRef] [Green Version]
  20. Kozyreva, O.; Pilipenko, V.; Krasnoperov, R.; Baddeley, L.; Sakharov, Y.A.; Dobrovolsky, M. Fine structure of substorm and geomagnetically induced currents. Ann. Geophys. 2020, 63, GM219. [Google Scholar] [CrossRef]
  21. Connors, M.; Rostoker, G.; Sofko, G.; McPherron, R.L.; Henderson, M.G. Ps 6 disturbances: Relation to substorms and the auroral oval. Ann. Geophys. 2003, 21, 493–508. [Google Scholar] [CrossRef] [Green Version]
  22. Apatenkov, S.V.; Pilipenko, V.A.; Gordeev, E.I.; Viljanen, A.; Juusola, L.; Belakhovsky, V.B.; Sakharov, Y.A.; Selivanov, V.N. Auroral omega bands are a significant cause of large geomagnetically induced currents. Geophys. Res. Lett. 2020, 47, e2019GL086677. [Google Scholar] [CrossRef]
  23. Chinkin, V.E.; Soloviev, A.A.; Pilipenko, V.A.; Engebretson, M.J.; Sakharov, Y.A. Determination of vortex current structure in the high-latitude ionosphere with associated GIC bursts from ground magnetic data. J. Atmos. Sol. Terr. Phys. 2021, 212, 105514. [Google Scholar] [CrossRef]
  24. Engebretson, M.J.; Steinmetz, E.S.; Posch, J.L.; Pilipenko, V.A.; Moldwin, M.B.; Connors, M.G.; Boteler, D.H.; Mann, I.R.; Hartinger, M.D.; Weygand, J.M.; et al. Nighttime magnetic perturbation events observed in Arctic Canada: 2. Multiple-instrument observations. J. Geophys. Res. Space Phys. 2019, 124, 7459–7476. [Google Scholar] [CrossRef]
  25. Yagova, N.V.; Pilipenko, V.A.; Sakharov, Y.A.; Selivanov, V.A. Spatial scale of geomagnetic Pc5/Pi3 pulsations as a factor of their efficiency in generation of geomagnetically induced currents. Earth Planets Space 2021, 73, 88. [Google Scholar] [CrossRef]
  26. Puthe, C.; Kuvshinov, A. Towards quantitative assessment of the hazard from space weather: Global 3D-modellings of the electric field induced by a realistic geomagnetic storm. Earth Planets Space 2013, 65, 1017–1025. [Google Scholar] [CrossRef]
  27. Kaufman, A.A.; Keller, G.V. The Magnetotelluric Sounding Method; Elsevier Science: New York, NY, USA, 1981; 595p. [Google Scholar]
  28. Berdichevsky, M.N.; Dmitriev, V.I. The Magnetotelluric Response Functions. In Models and Methods of Magnetotellurics; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  29. Marshalko, E.; Kruglyakov, M.; Kuvshinov, A.; Murphy, B.S.; Rastätter, L.; Ngwira, C.; Pulkkinen, A. Exploring the influence of lateral conductivity contrasts on the storm time behavior of the ground electric field in the eastern United States. Space Weather 2020, 18, e2019SW002216. [Google Scholar] [CrossRef] [Green Version]
  30. Alekseev, D.; Kuvshinov, A.; Palshin, N. Compilation of 3D global conductivity model of the Earth for space weather applications. Earth Planets Space 2015, 67, 108. [Google Scholar] [CrossRef] [Green Version]
  31. Korja, T.; Engels, M.; Zhamaletdinov, A.A.; Kovtun, A.A.; Palshin, N.N.; Smirnov, M.Y.; Tokarev, A.D.; Asming, V.E.; Vanyan, L.L.; Vardaniants, I.L.; et al. Crustal conductivity in Fennoscandia—A compilation of a database on crustal conductance in the Fennoscandian Shield. Earth Planets Space 2002, 54, 535–558. [Google Scholar] [CrossRef] [Green Version]
  32. Kuvshinov, A.; Grayver, A.; Tøffner-Clausen, L.; Olsen, N. Probing 3-D electrical conductivity of the mantle using 6 years of Swarm, CryoSat-2 and observatory magnetic data and exploiting matrix Q-responses approach. Earth Planets Space 2021, 73, 67. [Google Scholar] [CrossRef]
  33. Kozyreva, O.; Pilipenko, V.; Sokolova, E.; Sakharov, Y.; Epishkin, D. Geomagnetic and telluric field variability as a driver of geomagnetically induced currents. In Problems of Geocosmos—2018, Springer Proceedings in Earth and Environmental Sciences; Yanovskaya, T.B., Kosterov, A., Bobrov, N.Y., Divin, A.V., Saraev, A.K., Zolotova, N.V., Eds.; Springer Nature: Cham, Switzerland, 2020; pp. 297–307. [Google Scholar] [CrossRef]
  34. Lehtinen, M.; Pirjola, R. Currents produced in earthed conductor networks by geomagnetically induced electric fields. Ann. Geophys. 1985, 3, 479–484. [Google Scholar]
  35. Pulkkinen, A.; Pirjola, R.; Boteler, D.; Viljanen, A.; Yegorov, I. Modelling of space weather effects on pipelines. J. Appl. Geophys. 2001, 48, 233–256. [Google Scholar] [CrossRef]
  36. Boteler, D.H. A new versatile method for modelling geomagnetic induction in pipelines. Geophys. J. Int. 2013, 193, 98–109. [Google Scholar] [CrossRef] [Green Version]
  37. Trichtchenko, L.; Boteler, D.H. Modelling of geomagnetic induction in pipelines. Ann. Geophys. 2020, 20, 1063–1072. [Google Scholar] [CrossRef]
  38. Trichtchenko, L.; Boteler, D.H. Specification of geomagnetically induced electric fields and currents in pipelines. J. Geophys. Res. 2001, 106, 21039–21048. [Google Scholar] [CrossRef]
  39. Hejda, P.; Bochnicek, J. Geomagnetically induced pipe-to-soil voltages in the Czech oil pipelines during October–November 2003. Ann. Geophys. 2005, 23, 3089–3093. [Google Scholar] [CrossRef] [Green Version]
  40. Gummow, R.; Eng, P. GIC effects on pipeline corrosion and corrosion control systems. J. Atmos. Sol. Terr. Phys. 2002, 64, 1755–1764. [Google Scholar] [CrossRef]
Figure 1. Stations in the Arctic zone of the Russian Federation monitoring geomagnetic field variations.
Figure 1. Stations in the Arctic zone of the Russian Federation monitoring geomagnetic field variations.
Applsci 12 03755 g001
Figure 2. An example of irregular quasiperiodic Pi3 pulsations against the background of the expansion phase of a substorm on 9 February 2019 (left) and 2 March 2019 (right). The top panel in both plots shows variations in the solar wind plasma density Np, solar wind velocity V, and vertical component of the interplanetary magnetic field Bz. Stack-plot magnetograms of the H-component from stations around Yamal Peninsula are shown in the bottom panel. Geomagnetic coordinates of the stations are indicated on the lower panels next to the station code. On the magnetograms, empty triangles show local noon at each station, and black diamonds show local midnight.
Figure 2. An example of irregular quasiperiodic Pi3 pulsations against the background of the expansion phase of a substorm on 9 February 2019 (left) and 2 March 2019 (right). The top panel in both plots shows variations in the solar wind plasma density Np, solar wind velocity V, and vertical component of the interplanetary magnetic field Bz. Stack-plot magnetograms of the H-component from stations around Yamal Peninsula are shown in the bottom panel. Geomagnetic coordinates of the stations are indicated on the lower panels next to the station code. On the magnetograms, empty triangles show local noon at each station, and black diamonds show local midnight.
Applsci 12 03755 g002
Figure 3. An example of localized magnetic disturbances (MPE) recorded on 9 September 2019 (left-hand) and 17 December 2019 (right-hand). The format of the figure is the same as in Figure 2.
Figure 3. An example of localized magnetic disturbances (MPE) recorded on 9 September 2019 (left-hand) and 17 December 2019 (right-hand). The format of the figure is the same as in Figure 2.
Applsci 12 03755 g003
Figure 4. Magnetograms (H component) from stations along the latitudinal profile 150°: VIZ-BEY-AMD-NAD during 14–22 UT on 6 January 2019. Black diamonds denote local midnight.
Figure 4. Magnetograms (H component) from stations along the latitudinal profile 150°: VIZ-BEY-AMD-NAD during 14–22 UT on 6 January 2019. Black diamonds denote local midnight.
Applsci 12 03755 g004
Figure 5. The conductivity profiles σ(z) and impedance dependence on frequency Z(f) for stations at Yamal derived from the global conductivity model. Below 100 km, profiles for each stations merge into one.
Figure 5. The conductivity profiles σ(z) and impedance dependence on frequency Z(f) for stations at Yamal derived from the global conductivity model. Below 100 km, profiles for each stations merge into one.
Applsci 12 03755 g005
Figure 6. Combined magnetic and telluric field variations at station BEY on 6 January 2019 in the following order (from top to bottom): magnetic H and E components, magnetic field variability dH/dt and dE/dt, and electrotelluric EH and EE components.
Figure 6. Combined magnetic and telluric field variations at station BEY on 6 January 2019 in the following order (from top to bottom): magnetic H and E components, magnetic field variability dH/dt and dE/dt, and electrotelluric EH and EE components.
Applsci 12 03755 g006
Figure 7. The amplitude spectra of horizontal magnetic (H and E components) and telluric (EH and EE components) field variations at station BEY during the 6 January 2019 event. Bottom panel shows the spectra of the magnetic field variability dH/dt and dE/dt.
Figure 7. The amplitude spectra of horizontal magnetic (H and E components) and telluric (EH and EE components) field variations at station BEY during the 6 January 2019 event. Bottom panel shows the spectra of the magnetic field variability dH/dt and dE/dt.
Applsci 12 03755 g007
Figure 8. Schematic illustration of the formation of GIC-effective frequency interval in the spectral content of geomagnetic variations.
Figure 8. Schematic illustration of the formation of GIC-effective frequency interval in the spectral content of geomagnetic variations.
Applsci 12 03755 g008
Figure 9. Map of the spatial distribution of the average conductivityσ [S/m] in the upper 10 km of the conductivity model [30] in the territory of the Russian Federation.
Figure 9. Map of the spatial distribution of the average conductivityσ [S/m] in the upper 10 km of the conductivity model [30] in the territory of the Russian Federation.
Applsci 12 03755 g009
Figure 10. Map of the horizontal geoelectric field magnitude |EH|, excited by uniform magnetic field variations with period T = 100 s and magnitude 1 nT.
Figure 10. Map of the horizontal geoelectric field magnitude |EH|, excited by uniform magnetic field variations with period T = 100 s and magnitude 1 nT.
Applsci 12 03755 g010
Figure 11. Map of the horizontal geoelectric field magnitude |EH|, excited by a uniform magnetic field with period T = 1000 s and magnitude of 1 nT.
Figure 11. Map of the horizontal geoelectric field magnitude |EH|, excited by a uniform magnetic field with period T = 1000 s and magnitude of 1 nT.
Applsci 12 03755 g011
Table 1. Geographic and corrected geomagnetic coordinates of stations in the Russian Arctic for the epoch 2015.
Table 1. Geographic and corrected geomagnetic coordinates of stations in the Russian Arctic for the epoch 2015.
StationCodeGeographic
Coordinates
Geomagnetic CoordinatesSupporting Institution
LatitudeLongitudeLatitudeLongitude
AmdermaAMD69.4761.4265.4137.7AARI
BaranovBRN 79.27101.7574.2174.3AARI
BarentsburgBBG78.0714.2575.6108.6AARI
BelyBEY73.3070.0069.5146.8IZMIRAN
CheluskinCCS77.72104.2871.3175.0AARI
ChokurdakhCHD70.62147.8964.8212.4ИКФИА
DiksonDIK73.5480.5667.9155.7AARI
HeissHIS80.6258.0574.7145.1AARI
KamennyCKA68.5073.6063.3148.2AARI
KharasaveyKHS71.1166.8667.5143.1AARI
KotelnyKTN75.94137.7169.9201.2MAGDAS
LovozeroLOZ 67.9735.0864.1114.9AARI
LoparskayaLOP68.2533.0864.1114.1IZMIRAN
NadymNAD65.5372.5061.8145.6IZMIRAN
NorilskNOK69.4088.4064.3162.2AARI
PevekPBK 70.08170.9065.0228.8AARI
SabettaSAB71.4272.1365.2147.1AARI
SalekhardSKD66.5266.6763.0141.7IZMIRAN
SeyahaSEY70.1072.5064.8147.7IZMIRAN
SchmidtCPS 68.88180.6364.2235.1IKIR
TiksiTIK71.59128.7865.7196.9ИКФИА
VizeVIZ79.2976.5874.9154.9AARI
White SeaWSE66.5533.1063.2111.5GC/MSU
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kozyreva, O.V.; Pilipenko, V.A.; Marshalko, E.E.; Sokolova, E.Y.; Dobrovolsky, M.N. Monitoring of Geomagnetic and Telluric Field Disturbances in the Russian Arctic. Appl. Sci. 2022, 12, 3755. https://doi.org/10.3390/app12083755

AMA Style

Kozyreva OV, Pilipenko VA, Marshalko EE, Sokolova EY, Dobrovolsky MN. Monitoring of Geomagnetic and Telluric Field Disturbances in the Russian Arctic. Applied Sciences. 2022; 12(8):3755. https://doi.org/10.3390/app12083755

Chicago/Turabian Style

Kozyreva, Olga V., Vyacheslav A. Pilipenko, Elena E. Marshalko, Elena Yu. Sokolova, and Mikhail N. Dobrovolsky. 2022. "Monitoring of Geomagnetic and Telluric Field Disturbances in the Russian Arctic" Applied Sciences 12, no. 8: 3755. https://doi.org/10.3390/app12083755

APA Style

Kozyreva, O. V., Pilipenko, V. A., Marshalko, E. E., Sokolova, E. Y., & Dobrovolsky, M. N. (2022). Monitoring of Geomagnetic and Telluric Field Disturbances in the Russian Arctic. Applied Sciences, 12(8), 3755. https://doi.org/10.3390/app12083755

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop