Search

Extreme space weather events can produce geomagnetically induced currents (GICs) that flow along long conductor systems such as power grids and pipelines. GICs can interrupt the operation of these systems by damaging transformers and increasing the rate of pipeline corrosion. GICs depend on the enhancement of the geoelectric field, where the magnitude, geographic location, and occurrence depend on the solar wind-magnetospheric-ionospheric coupling processes and the ground conductivity. Thus, countries at high latitudes are potentially more vulnerable to GICs due to higher geomagnetic activity.Magnetic field measurements (e.g. dB/dt) have been used for a long time as a proxy for GICs with some success. However, studying the behaviour of the geoelectric field is required to fully understand the GIC response during extreme space weather events. In this study, the geoelectric field was computed using data from the IMAGE network and a 1D conductivity model (SMAP) across Sweden. A 19-year statistical analysis of the daily maximum magnitude of the geoelectric field (E) has been performed for solar cycles 23 and 24. We examine the temporal and spatial distribution of the geoelectric field in Sweden to determine the importance of including the ground conductivity when assessing the strongest events and their implications to GICs. We found that the daily maximum E is more frequently registered in the dusk sector related to the eastward convection electrojet and a second relative maximum is observed in the dawn sector related to westward convection electrojet. The stronger E values are related to well-known extreme space weather events. However, these events produce different responses at different latitudes due to the changes in ground conductivity and ionospheric response. Therefore, the strongest geoelectric fields at different geographical locations are not all driven by the same events.

In this study, we perform rigorous three-dimensional (3-D) ground electric field (GEF) modeling in Fennoscandia for three days of the Halloween geomagnetic storm (29-31 October 2003) using magnetic field data from the IMAGE magnetometer network and a 3-D conductivity model of the region. To explore the influence of the inducing source model on the 3-D GEF simulations, we considered three different approaches to approximate the source. Within the first two approaches, the source varies laterally, whereas in the third method, the GEF is calculated by implementing the time-domain realization of the magnetotelluric intersite impedance method. We then compare GEF-based geomagnetically induced current (GIC) with observations at the Mäntsälä natural gas pipeline recording point. We conclude that a high correlation between modeled and recorded GIC is observed for all considered approaches. The highest correlation is achieved when performing a 3-D GEF simulation using a “conductivity-based” laterally nonuniform inducing source. Our results also highlight the strong dependence of the GEF on the earth’s conductivity distribution.

We have examined the most intense external (magnetospheric and ionospheric) and internal (induced) |dH/dt| (amplitude of the 10 s time derivative of the horizontal geomagnetic field) events observed by the high-latitude International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometers between 1994 and 2018. While the most intense external |dH/dt| events at adjacent stations typically occurred simultaneously, the most intense internal (and total) |dH/dt| events were more scattered in time, most likely due to the complexity of induction in the conducting ground. The most intense external |dH/dt| events occurred during geomagnetic storms, among which the Halloween storm in Oct 2003 featured prominently, and drove intense geomagnetically induced currents (GIC). Events in the prenoon local time sector were associated with sudden commencements (SC) and pulsations, and the most intense |dH/dt| values were driven by abrupt changes in the eastward electrojet due to solar wind dynamic pressure increase or decrease. Events in the premidnight and dawn local time sectors were associated with substorm activity, and the most intense |dH/dt| values were driven by abrupt changes in the westward electrojet, such as weakening and poleward retreat (premidnight) or undulation (dawn). Despite being associated with various event types and occurring at different local time sectors, there were common features among the drivers of most intense external |dH/dt| values: pre-existing intense ionospheric currents (SC events were an exception) that were abruptly modified by sudden changes in the magnetospheric magnetic field configuration. While proper description of the fast changes during SC events appears to require 1 s data, pulsations and substorms may be sufficiently described by 10 s |dH/dt|. 1 min data, however, significantly underestimates the |dH/dt| peaks. Our results contribute towards the ultimate goal of reliable forecasts of dH/dt and GIC.

Space weather, like solar eruptions, can be hazardous to Earth’s electric grids via geomagnetically induced currents (GIC). In worst cases they can even cause city-wide power outages. GIC is a complicated phenomenon, closely related to the time derivative of the geomagnetic field. However, behavior the time derivative is chaotic and has proven to be challenging to predict. In this study we look at the geomagnetic field orientations at different magnetometer stations in the Fennoscandian region during active space weather conditions. We aim to characterize the magnetic field behavior, to better understand the drivers behind strong GIC events. One of our main findings is that the direction of time derivative of the geomagnetic field has a very short “reset time“, only a few minutes. We conclude that this result gives insight on the time scale of the ionospheric current systems, which are the primary driver behind the time derivative’s behavior.

During auroral substorms, the electric currents flowing in the ionosphere change rapidly, and a large amount of energy is dissipated in the auroral ionosphere. An important part of the auroral current system is the auroral electrojets whose profiles can be estimated from magnetic field measurements from low-earth orbit satellites. In this paper, we combine electrojet data derived from the Swarm satellite mission of the European Space Agency with the substorm database derived from the SuperMAG ground magnetometer network data. We organize the electrojet data in relation to the location and time of the onset and obtain statistics for the development of the integrated current and latitudinal location for the auroral electrojets relative to the onset. The major features of the behaviour of the westward electrojet are found to be in accordance with earlier studies of field-aligned currents and ground magnetometer observations of substorm temporal statistics. In addition, we show that, after the onset, the latitudinal location of the maximum of the westward electrojet determined from Swarm satellite data is mostly located close to the SuperMAG onset latitude in the local time sector of the onset regardless of where the onset happens. We also show that the SuperMAG onset corresponds to a strengthening of the order of 100 kA in the amplitude of the median of the westward integrated current in the Swarm data from 15 min before to 15 min after the onset.

Solar eruptions and other types of space weather effects can pose a hazard to the high voltage power grids via geomagnetically induced currents (GICs). In worst cases, they can even cause large-scale power outages. GICs are a complex phenomenon, closely related to the time derivative of the geomagnetic field. However, the behavior of the time derivative is chaotic and has proven to be tricky to predict. In our study, we look at the dynamics of the geomagnetic field during active space weather. We try to characterize the magnetic field behavior, to better understand the drivers behind strong GIC events. We use geomagnetic data from the IMAGE (International Monitor for Auroral Geomagnetic Effect) magnetometer network between 1996 and 2018. The measured geomagnetic field is primarily produced by currents in the ionosphere and magnetosphere, and secondarily by currents in the conducting ground. We use the separated magnetic field in our analysis. The separation of the field means that the measured magnetic field is computationally divided into external and internal parts corresponding to the ionospheric and telluric origin, respectively. We study the yearly directional distributions of the baseline subtracted, separated horizontal geomagnetic field, Delta H, and its time derivative, d Delta H/dt. The yearly distributions do not have a clear solar cycle dependency. The internal field distributions are more scattered than the external field. There are also clear, station-specific differences in the distributions related to sharp conductivity contrasts between continental and ocean regions or to inland conductivity anomalies. One of our main findings is that the direction of d Delta H/dt has a very short "reset time", around 2 min, but Delta H does not have this kind of behavior. These results hold true even with less active space weather conditions. We conclude that this result gives insight into the time scale of ionospheric current systems, which are the primary driver behind the time derivative's behavior. It also emphasizes a very short persistence of d Delta H/dt compared to Delta H, and highlights the challenges in forecasting d Delta H/dt (and GIC).

During auroral substorms the electric currents flowing in the ionosphere change rapidly and a large amount of energy is dissipated in the auroral ionosphere. An important part of the auroral current system are the auroral electrojets whose profiles can be estimated from magnetic field measurements from Low Earth Orbit satellites. In this paper we combine electrojet data derived from the Swarm satellite mission of ESA with the substorm database derived from the SuperMAG ground magnetometer network data. We organize the electrojet data in relation to the location and time of the onset and obtain statistics for the development of the integrated current and latitudinal location for the auroral electrojets relative to the onset. The major features of the behaviour of the westward electrojet are found to be in accordance with earlier studies of field aligned currents and ground magnetometer observations of substorm time statistics. In addition we show that after the onset the latitudinal location of the maximum of the westward electrojet determined from Swarm satellite data is mostly located close to the SuperMAG onset latitude in the local time sector of the onset regardless of where the onset happens. We also show that the SuperMAG onset corresponds to a strengthening of the order of 100 kA in the amplitude of the median of the westward integrated current in the Swarm data from 15 minutes before to 15 minutes after the onset.

We must be able to predict and mitigate against geomagnetically induced current (GIC) effects to minimize socio-economic impacts. This study employs the space weather modeling framework (SWMF) to model the geomagnetic response over Fennoscandia to the September 7-8, 2017 event. Of key importance to this study is the effects of spatial resolution in terms of regional forecasts and improved GIC modeling results. Therefore, we ran the model at comparatively low, medium, and high spatial resolutions. The virtual magnetometers from each model run are compared with observations from the IMAGE magnetometer network across various latitudes and over regional-scales. The virtual magnetometer data from the SWMF are coupled with a local ground conductivity model which is used to calculate the geoelectric field and estimate GICs in a Finnish natural gas pipeline. This investigation has lead to several important results in which higher resolution yielded: (1) more realistic amplitudes and timings of GICs, (2) higher amplitude geomagnetic disturbances across latitudes, and (3) increased regional variations in terms of differences between stations. Despite this, substorms remain a significant challenge to surface magnetic field prediction from global magnetohydrodynamic modeling. For example, in the presence of multiple large substorms, the associated large-amplitude depressions were not captured, which caused the largest model-data deviations. The results from this work are of key importance to both modelers and space weather operators. Particularly when the goal is to obtain improved regional forecasts of geomagnetic disturbances and/or more realistic estimates of the geoelectric field.

Geomagnetically Induced Currents (GICs) are a space weather hazard that can negatively impact large ground-based infrastructures such as power lines, pipelines, and railways. They are driven by the dynamic spatiotemporal behaviour of currents flowing in geospace, which drive rapid geomagnetic disturbances on the ground. In some cases, geomagnetic disturbances are highly localised and spatially structured due to the dynamical behaviour of geospace currents and magnetosphere-ionosphere (M-I) coupling dynamics, which are complex and often unclear.In this work, we investigate and quantify the spatial structure of large geomagnetic depressions exceeding several hundred nT according to the 10 strongest events measured over Fennoscandia by IMAGE. Using ground magnetometer measurements we connect these spatially structured geomagnetic disturbances to possible M-I coupling processes and identify their likely magnetospheric origin. In addition, the ability for these disturbances to drive large GICs is assessed by calculating their respective geoelectric fields in Sweden using the SMAP ground conductivity model. To compliment the observations, we also utilise high resolution runs (>7 million cells) of the Space Weather Modeling Framework (SWMF) to determine to what extent global MHD models can capture this behaviour.

Geomagnetic variations are mainly produced by external currents in the ionosphere and magnetosphere, and secondarily by induced (internal/telluric) currents in the conducting Earth. Large geomagnetically induced currents (GIC) are associated with large time derivatives of the horizontal magnetic field. Recent results show that the time derivative is typically dominated by the contribution from the telluric currents. Our study aims to find measures to quantify the behaviour of external and internal currents and their time derivatives during large GIC events. Results of this study show that strong external currents have quite narrow directional distributions. Angular variation is larger for internal currents, and especially for their time derivatives. For external currents angular variation is larger at higher latitudes. Similar behaviour is not seen with internal currents.

Non

Faraday's law of induction is responsible for setting up a geoelectric field due to the variations in the geomagnetic field caused by ionospheric currents. This drives geomagnetically induced currents (GICs) which flow in large ground‐based technological infrastructure such as high‐voltage power lines. The geoelectric field is often a localized phenomenon exhibiting significant variations over spatial scales of only hundreds of kilometers. This is due to the complex spatiotemporal behavior of electrical currents flowing in the ionosphere and/or large gradients in the ground conductivity due to highly structured local geological properties. Over some regions, and during large storms, both of these effects become significant. In this study, we quantify the regional variability of dB/dt using closely placed IMAGE stations in northern Fennoscandia. The dependency between regional variability, solar wind conditions, and geomagnetic indices are also investigated. Finally, we assess the significance of spatial geomagnetic variations to modeling GICs across a transmission line. Key results from this study are as follows: (1) Regional geomagnetic disturbances are important in modeling GIC during strong storms; (2) dB/dt can vary by several times up to a factor of three compared to the spatial average; (3) dB/dt and its regional variation is coupled to the energy deposited into the magnetosphere; and (4) regional variability can be more accurately captured and predicted from a local index as opposed to a global one. These results demonstrate the need for denser magnetometer networks at high latitudes where transmission lines extending hundreds of kilometers are present.

Geomagnetically induced currents (GIC) are directly described by ground electric fields, but estimating them is time-consuming and requires knowledge of the ionospheric currents as well as the three-dimensional distribution of the electrical conductivity of the Earth. The time derivative of the horizontal component of the ground magnetic field (dH/dt) is closely related to the electric field via Faraday's law, and provides a convenient proxy for the GIC risk. However, forecasting dH/dt still remains a challenge. We use 25 years of 10 s data from the North European International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network to show that part of this problem stems from the fact that instead of the primary ionospheric currents, the measured dH/dt is dominated by the signature from the secondary induced telluric currents nearly at all IMAGE stations. The largest effects due to telluric currents occur at coastal sites close to highly-conducting ocean water and close to near-surface conductivity anomalies. The secondary magnetic field contribution to the total field is a few tens of percent, in accordance with earlier studies. Our results have been derived using IMAGE data and are thus only valid for the involved stations. However, it is likely that the main principle also applies to other areas. Consequently, it is recommended that the field separation into internal (telluric) and external (ionospheric and magnetospheric) parts is performed whenever feasible, i.e., a dense observation network is available.

High speed streams (HSS) and associated co-rotating interaction regions (CIR) in the solar wind are one of the major drivers of geomagnetic activity, especially during declining phases of sunspot cycles and near sunspot minima. We have identified 51 HSS/CIR driven geomagnetic storms that coincide with a Dst drop to less than -50nT during the period 2009-2018 and we investigate their impact on ionospheric current systems. Our approach is to study the evolution of the global scale current systems, i.e. the auroral electrojets and Region-1/2 field-aligned currents (FAC), with the SuperMAG magnetometers and AMPERE satellite data, respectively. The events are studied with a superposed epoch analysis centered at the storm onset to see the general behavior of the current system globally and in four different MLT sectors: noon, dusk, midnight and dawn. A minor enhancement of the integrated FAC was observed in the midnight, dawn and dusk sector 3 hours before the storm onset. The largest FAC and variability was observed in the dusk sector, and the integrated FAC maximum occurred in the middle of the storm main phase, 4 hours before the Dst minimum. This result will be compared to the evolution and behavior of the electrojet currents from superMAG. In the future a similar study will be conducted for ICME geomagnetic storms and compared to the HSS/CIR-related storms.

Fast geomagnetic variations of periods from seconds to hours and days are primarily produced by currents in the ionosphere and magnetosphere. There is always an associated secondary (internal, telluric) current system induced in the conducting ground and contributing to the total variation field measured by ground magnetometers. Mathematically, it is possible to fully explain the variation field by two equivalent current systems, one at the ionospheric altitude and another just below the ground. In practice, this separation is feasible using dense magnetometer networks.A common way in space physics has been to implicitly neglect the internal part and interpret the ground field only in terms of ionospheric currents. As known from previous studies, this is often a reasonable assumption, since a typical internal contribution is about 30%. However, the situation is much different when the time derivative of the magnetic field (dB/dt) is considered. For the north European IMAGE magnetometer network, the internal part exceeds the external one nearly at all stations. The largest effects due to telluric currents occur at coastal sites close to highly-conducting ocean water and at inland locations close to highly-conducting near-surface anomalies.This finding gives a new viewpoint for studies of geomagnetically induced currents (GIC), which are closely related to dB/dt. One key question is to understand which are the ionospheric drivers of big GIC events. We will demonstrate how the telluric currents can strongly modify field variations and especially dB/dt, and how this is correspondingly seen in equivalent current patterns. Consequently, it is recommended that the field separation is performed whenever it is feasible, i.e. a dense observation network is available.

The strongest event of geomagnetically induced currents (GIC) detected by the North-West Russian GIC network occurred during the main phase of the magnetic storm on June 28-29, 2013. Extremely high values, 120 A, were recorded in the 330 kV transformers on Kola Peninsula in the 04--07 magnetic local time (MLT) sector. The Defense Meteorological Satellite Program (DMSP) spacecraft took a sequence of ultraviolet (UV) auroral images in the southern hemisphere and observed multiple omega bands. The ionospheric equivalent electric currents based on the International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network reveal a sequence of current vortex pairs moving eastward with the speed of 0.5-2.5 km/s, that fits to the electrodynamics scheme of omega bands. Although the temporal variations of the associated current system are slow, the omega bands can be responsible for strong magnetic variations and GIC due to fast propagations of currents in the azimuthal direction. The first steps towards the statistica study of the highest GIC recorded at Vykhodnoy transformer show that about 50% of events have properties similar to the comprehensively studied 29 June 2013 case.

Geomagnetically induced currents (GICs) are directly described by ground electric fields, but estimating them is time-consuming and requires knowledge of the ionospheric currents and the three-dimensional (3D) distribution of the electrical conductivity of the Earth. The time derivative of the horizontal component of the ground magnetic field (dH∕dt) is closely related to the electric field via Faraday's law and provides a convenient proxy for the GIC risk. However, forecasting dH∕dt still remains a challenge. We use 25 years of 10 s data from the northern European International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network to show that part of this problem stems from the fact that, instead of the primary ionospheric currents, the measured dH∕dt is dominated by the signature from the secondary induced telluric currents at nearly all IMAGE stations. The largest effects due to telluric currents occur at coastal sites close to high-conducting ocean water and close to near-surface conductivity anomalies. The secondary magnetic field contribution to the total field is a few tens of percent, in accordance with earlier studies. Our results have been derived using IMAGE data and are thus only valid for the stations involved. However, it is likely that the main principle also applies to other areas. Consequently, it is recommended that the field separation into internal (telluric) and external (ionospheric and magnetospheric) parts is performed whenever feasible (i.e., a dense observation network is available).

Geomagnetically induced currents (GICs) have been long recognized as a ground effect arising from a chain of space weather events. GICs have been measured and modeled in many countries, resulting in a considerable amount of data. Previous statistical analyses have proposed various types of distribution functions to fit long-term GICs data sets. However, these extensive statistical approaches have been only partially successful in fitting the data sets. Here we use modeled GICs data sets calculated in four countries (Brazil, South Africa, United Kingdom, and Finland) using data from solar cycle 23 to show a plausible function based on a nonextensive statistical model of the q-exponential Tsallis function. The fitted q-exponential parameter is approximately the same for all locations, and the Lilliefors test shows good agreement with the q-exponential fits. From this fit, we compute that the likely numbers of extreme GICs events over the next ten solar cycles are 1–2 for both Finland and United Kingdom, at least one for Brazil and less than one event for South Africa. Our results indicate that the nonextensive statistics are a general characteristic of GICs, suggesting that the ground current intensity has a strong temporal correlation and long-range interaction.

This paper is the primary deliverable of the very first NASA Living With a Star Institute Working Group, Geomagnetically Induced Currents (GIC) Working Group. The paper provides a broad overview of the current status and future challenges pertaining to the science, engineering, and applications of the GIC problem. Science is understood here as the basic space and Earth sciences research that allows improved understanding and physics-based modeling of the physical processes behind GIC. Engineering, in turn, is understood here as the “impact” aspect of GIC. Applications are understood as the models, tools, and activities that can provide actionable information to entities such as power systems operators for mitigating the effects of GIC and government agencies for managing any potential consequences from GIC impact to critical infrastructure. Applications can be considered the ultimate goal of our GIC work. In assessing the status of the field, we quantify the readiness of various applications in the mitigation context. We use the Applications Readiness Level (ARL) concept to carry out the quantification.

The measurement and collection of digital magnetic field data in Europe extends back to the 1970s, providing over 30 years of data for the analysis of severe space weather. Although paper records can potentially extend these datasets back by over a century, few digitized records are currently available for use in extreme studies. Therefore we rely on theoretical arguments and modelling to elucidate the largest likely variations of the magnetic field. We assess the relationship, during the three largest storms in the digital era, between variations in the horizontal magnetic field and the largest measured Dst index to estimate likely magnetic variations for more extreme storms in northern and mid-latitude Europe. We examine how Geomagnetically Induced Currents (GIC) flow in the UK and French networks during recent severe storms and analyse the sensitivity of these flows to changes in grid parameters. The maximum GIC computed at any one node in the French and UK grids are 44 A and 208 A respectively. Sensitivity tests show that while gross changes of the whole network structure, such as disconnecting parts of the network, reduces the mean GIC per node, changes in GIC at individual nodes have distinct behaviors implying that local effects are network-dependent and require detailed modelling to sufficiently characterize GIC. In addition, the scale factors we have derived allow GIC results from recent storms to be up-scaled to estimate the potential risk to the system from more extreme events, such as the Carrington storm in 1859.

Between 7 and 8 September 2017, Earth experienced extreme space weather events. We have combined measurements made by the IMAGE magnetometer array, ionospheric equivalent currents, geomagnetically ...

High-frequency ( ≈ minutes) variability of ground magnetic fields is caused by ionospheric and magnetospheric processes driven by the changing solar wind. The varying magnetic fields induce electrical fields that cause currents to flow in man-made conductors like power grids and pipelines. Under extreme conditions the geomagnetically induced currents (GIC) may be harmful to the power grids. Increasing our understanding of the extreme events is thus important for solar-terrestrial science and space weather. In this work 1-min resolution of the time derivative of measured local magnetic fields (|dBh∕dt|) and computed electrical fields (Eh), for locations in Europe, have been analysed with extreme value analysis (EVA). The EVA results in an estimate of the generalized extreme value probability distribution that is described by three parameters: location, width, and shape. The shape parameter controls the extreme behaviour. The stations cover geomagnetic latitudes from 40 to 70° N. All stations included in the study have contiguous coverage of 18 years or more with 1-min resolution data. As expected, the EVA shows that the higher latitude stations have higher probability of large |dBh∕dt| and |Eh| compared to stations further south. However, the EVA also shows that the shape of the distribution changes with magnetic latitude. The high latitudes have distributions that fall off faster to zero than the low latitudes, and upward bounded distributions can not be ruled out. The transition occurs around 59–61° N magnetic latitudes. Thus, the EVA shows that the observed series north of ≈ 60° N have already measured values that are close to the expected maxima values, while stations south of ≈ ° N will measure larger values in the future.

To learn about the geometry and sources of the ionospheric current systems which generate strong geomagnetically induced currents, we categorize differential equivalent current systems (DEC) for events with strong dB/dt by decomposing them into the contributions of electrojet-type and vortex-type elementary systems. By solving the inverse problem we obtain amplitudes and locations of these elementary current systems. One-minute differences of the geomagnetic field values at the IMAGE magnetometer network in 1996–2000 are analysed to study the spatial distributions of large dB/dt events. The relative contributions of the two components are evaluated. In particular, we found that the majority of the strongest dB/dt events (100–1000nT/min) appear to be produced by the vortex-type current structures and most of them occur in the morning LT hours, probably caused by the Ps6 pulsation events associated with auroral omega structures. For strong dB/dt events the solar wind parameters are shifted toward strong (tens nT) southward IMF, enhanced velocity and dynamic pressure, in order for the main phase of the magnetic storms to occur. Although these events appear mostly during magnetic storms when the auroral oval greatly expands, the area of large dB/dt stays in the middle part of the auroral zone; therefore, it is connected to the processes taking part in the middle of the magnetosphere rather than in its innermost region populated by the ring current. Key words. Geomagnetism and paleomagnetism (rapid time variations) – Ionosphere (auroral ionosphere; ionospheric disturbances)

We study the role of ionospheric induction in different commonly observed ionospheric situations. These include an intensifying electrojet, westward travelling surge (WTS) and Ω-band. We use data based, realistic models for these phenomena and calculate the inductive electric fields that are created due to the temporal variations of ionospheric currents. The ionospheric induction problem is solved using a new calculation technique that can handle non-uniform, time-dependent conductances and electric fields of any geometry. We find that in some situations inductive effects are not negligible and the ionospheric electric field is not a pure potential field, but has a significant induced rotational part. In the WTS and Ω-band models the induced electric field is concentrated in a small area, where the time derivatives are largest. In the electrojet model the induced field is significant over a large part of the jet area. In these examples the induced electric field has typical values of few mV/m, which amounts to several tens of percents of the potential electric field present at the same locations. The induced electric field is associated with ionospheric and field aligned currents (FAC), that modify the overall structure of the current systems. Especially the induced FAC are often comparable to the non-inductive FAC, and may thus modify the coupling between the ionosphere and magnetosphere in the most dynamical situations. We also present some examples with very simple ionospheric current systems, where the effect of different ionospheric parameters on the induction process is studied.

One important contribution to the magnetic field measured at satellite altitude and at ground level comes from the external currents. We used the total field data sampled by the Overhauser Magnetometer on CHAMP and the horizontal magnetic field measurements of the IMAGE ground-based magnetometer network to study the ionospheric Hall current system in the auroral regions. For the CHAMP data a current model consisting of a series of lines and placed at a height of 110km is fitted to the magnetic field signature sampled on the passage across the polar region. The derived current distributions depend, among others, on season and on the local time of the satellite track. At dawn/dusk the auroral electrojets can be detected most clearly in the auroral regions. Their intensity and location are evidently correlated with the A E activity index. For a period of almost two years the results obtained from space and the currents determined from ground-based observations are studied. For the full IMAGE station array a newly-developed method of spherical elementary current systems (SECS) is employed to compute the 2-D equivalent current distribution, which gives a detailed picture of an area covering latitudes 60° – 80° N and 10° – 30° E in the auroral region. Generally, the current estimates from satellite and ground are in good agreement. The results of this survey clearly show the average dependence of the auroral electrojet on season and local time. This is particularly true during periods of increased auroral activity. The correlation coefficient of the results is close to one in the region of sizeable ionospheric current densities. Also the ratio of the current densities, as determined from above and below the ionosphere, is close to unity. It is the first time that the method of Hall current estimate from a satellite has been validated quantitatively by ground-based observations. Among others, this result is of interest for magnetic main field modelling, since it demonstrates that ground-based observations can be used to predict electrojet signatures in satellite magnetic field scalar data. Key words. Ionosphere (auroral Ionosphere; electric fields and currents; ionosphere-magnetosphere interactions)

In addition to field observations and numerical models, geomagnetic induction effects can be studied by scaled analogue model experiments. We present here results of analogue model studies of the auroral electrojet with an Earth model simulating the Arctic Ocean and inland conductivity structures in northern Fennoscandia. The main elements of the analogue model used were salt water simulating the host rock, an aluminium plate corresponding to the ocean and graphite pieces producing the inland highly conducting anomalies. The electrojet was a time-harmonic line current flowing at a (simulated) height of 100 km above northern Fennoscandia. The period simulated was 9 min. The analogue model results confirmed the well-known rapid increase of the vertical field when the coast is approached from the continent. The increase of the horizontal field due to induced ocean currents was demonstrated above the ocean, as well as the essentially negligible effect of these currents on the horizontal field on the continent. The behaviour of the magnetic field is explained with a simple two-dimensional thin-sheet model. The range, or the adjustment distance, of the ocean effect inland was found to be some hundreds of kilometers, which also agrees with earlier results of the Siebert-Kertz separation of IMAGE magnetometer data. The modelled inland anomalies evidently had too large conductivities, but on the other hand, their influence decayed on scales of only some tens of kilometers. Analogue model results, thin-sheet calculations, and field observations show that the induction effect on the horizontal magnetic field Bx near the electrojet is negligible. On the other hand, the vertical component Bz is clearly affected by induced currents in the ocean. Evidence of this is the shift of the zero point of Bz 0-1° southwards from the maximum of Bx. The importance of these results are discussed, emphasizing the determination of ionospheric currents.

Geomagnetically induced currents (GICs) have been long recognized as a ground effect arising from a chain of space weather events. GICs have been measured and modeled in many countries, resulting in a considerable amount of data. Previous statistical analyses have proposed various types of distribution functions to fit long-term GICs data sets. However, these extensive statistical approaches have been only partially successful in fitting the data sets. Here we use modeled GICs data sets calculated in four countries (Brazil, South Africa, United Kingdom, and Finland) using data from solar cycle 23 to show a plausible function based on a nonextensive statistical model of the q-exponential Tsallis function. The fitted q-exponential parameter is approximately the same for all locations, and the Lilliefors test shows good agreement with the q-exponential fits. From this fit, we compute that the likely numbers of extreme GICs events over the next ten solar cycles are 1–2 for both Finland and United Kingdom, at least one for Brazil and less than one event for South Africa. Our results indicate that the nonextensive statistics are a general characteristic of GICs, suggesting that the ground current intensity has a strong temporal correlation and long-range interaction.

One of the major challenges pertaining to extreme geomagnetic storms is to understand the basic processes associated with the development of dynamic magnetosphere-ionosphere currents, which generate large induced surface geoelectric fields. Previous studies point out the existence of localized peak geoelectric field enhancements during extreme storms. We examined induced global geoelectric fields derived from ground-based magnetometer recordings for 12 extreme geomagnetic storms between the years 1982 and 2005. For the present study two important extreme storms, 29 October 2003 and 13 March 1989, are shown. The primary purpose of this paper is to provide further evidence on the existence of localized peak geoelectric field enhancements and to show that the structure of the geoelectric field during these localized extremes at single sites can differ greatly from globally and regionally averaged fields. Although the physical processes that govern the development of these localized extremes are still not clear, we discuss some possible causes.

A solid understanding of the solar wind control of ground magnetic field disturbances is essential for utilizing the existing long time series of ground data to obtain information on solar wind-magnetosphere-ionosphere coupling. We have used 20 years of International Monitor for Auroral Geomagnetic Effects magnetometer data (54°–76° magnetic latitude) to study the solar wind control of the ionospheric equivalent current density and its time derivative ( ). We found that peaks at the premidnight and prenoon ends of the westward electrojet. The prenoon peak was most intense during fast solar wind and radial interplanetary magnetic field (IMF). The location of the peak was not affected by the IMF orientation but persisted at 8–10 magnetic local time and 70°–75° latitude, near the boundary between the westward and eastward electrojets. Sensitivity of this boundary to disturbances was suggested as a possible explanation for the persistent prenoon location of the peak. The premidnight peak was most intense during southward IMF orientation. While faster solar wind mainly resulted in more intense in the premidnight sector, stronger IMF caused the region of intense to spread to the postmidnight, dawn, and dusk sectors. A good correspondence was found between development of the nightside intensification and average substorm bulge and oval aurora as determined by Gjerloev et al. (2007). The bulge aurora covered the western end of the westward electrojet where the equivalent current also had a significant poleward component. The substorm oval aurora, on the other hand, extended eastward along the westward electrojet.

Interplanetary coronal mass ejections (ICMEs) compounded by high-speed plasma streams from coronal holes caused two intense geomagnetic storms on March 17–18, 2013 and 2015 during the current solar cycle. Ionospheric responses to the storms in the northern and southern hemispheres are compared in the context of solar wind coupling to the magnetosphere-ionosphere system. Phase scintillation is observed at high latitudes by arrays of high-rate GNSS Ionospheric Scintillation and TEC Monitors (GISTMs) and geodetic-quality GPS receivers sampling at 1 Hz. The high-rate GPS receivers are distributed in the northern and in the southern high latitudes with sparser coverage. In addition to GPS receivers, the high-latitude ionosphere dynamics is studied using arrays of ground-based instruments including HF radars, ionosondes, riometers, magnetometers, optical imagers as well as particle detectors and ultraviolet scanning imagers onboard the DMSP satellites.

The geomagnetism community is involved in the European Plate Observing System (EPOS), a European Research Infrastructure through which science communities will offer a number of services that will integrate to simplify cross-disciplinary research. The Geomagnetism community will provide data from geomagnetic observatories, from producers of geomagnetic indices and events, from geomagnetic models and from magneto-telluric observations. A number of these services (data from the INTERMAGNET network and the World Data Centre, indices and events from the International Service of Geomagnetic Indices and the access to the International Geomagnetic Reference Field and World Magnetic Model) will be integrated into EPOS systems in the first wave of services to be connected. This poster will describe the contribution from geomagnetism to EPOS. It will include a description of the data and services that the geomagnetic community will provide and also discuss how metadata will be made available from the community to the EPOS core IT systems. Finally it will describe how the provision of geomagnetic services in EPOS will be guided and governed by members of the community.

Understanding the capability of a simulation to reproduce observed features is a requirement for its use in operational space weather forecasting. We compare statistically ionospheric seasonal variations in the Grand Unified Magnetosphere-Ionosphere Coupling Simulation (GUMICS-4) global magnetohydrodynamic model with measurements. The GUMICS-4 data consist of a set of runs that was fed with real solar wind measurements and cover the period of 1 year. Ionospheric convection measurements are from the Super Dual Auroral Radar Network (SuperDARN) radars, and electric currents are derived from the magnetic field measured by the CHAMP satellite. Auroral electrojet indices are used to examine the disturbance magnetic field on ground. The signatures of electrodynamic coupling between the magnetosphere and ionosphere extend to lower latitudes in GUMICS-4 than in observations, and key features of the auroral ovals—the Region 2 field-aligned currents, electrojets, Harang discontinuity, and ring of enhanced conductivity—are not properly reproduced. The ground magnetic field is even at best about 5 times weaker than measurements, which can be a problem for forecasting geomagnetically induced currents. According to the measurements, the ionospheric electrostatic potential does not change significantly from winter to summer but field-aligned currents enhance, whereas in GUMICS-4, the electrostatic potential weakens from winter to summer but field-aligned currents do not change. This could be a consequence of the missing Region 2 currents: the Region 1 current has to close with itself across the polar cap, which makes it sensitive to solar UV conductivity. Precipitation energy and conductance peak amplitudes in GUMICS-4 agree with observations.

[1] We consider the model of the Finnish high-voltage power grid in 1978–1979, for which the accurate parameters are available for calculating geomagnetically induced currents (GIC). Moving the grid at different locations across Europe gives estimates of GIC levels in this region. For calculating the geoelectric field driving GIC, we use different layered models of the ground conductivity and 1 min geomagnetic data of the year 2003. The results show a clear concentration of large GIC in north Europe, where the peak values are about 3–5 times larger than in Central and South Europe, being up to about 200 A in this specific power grid. There are two factors contributing to this finding. First, geomagnetic variations are generally stronger in the north. Second, there are regions in the north with clearly smaller ground conductivities than typically at other areas. Both of these reasons lead to larger electric fields in the north. A very similar behavior of GIC is found in the case when a single-layered ground conductivity model is assumed everywhere. We also show that the geographic characteristics of GIC are quite insensitive to the details of the power grid model by modifying various parameters of the Finnish grid.

The geomagnetic storm of 17–18 March 2015 was caused by the impacts of a coronal mass ejection and a high-speed plasma stream from a coronal hole. The high-latitude ionosphere dynamics is studied using arrays of ground-based instruments including GPS receivers, HF radars, ionosondes, riometers, and magnetometers. The phase scintillation index is computed for signals sampled at a rate of up to 100 Hz by specialized GPS scintillation receivers supplemented by the phase scintillation proxy index obtained from geodetic-quality GPS data sampled at 1 Hz. In the context of solar wind coupling to the magnetosphere-ionosphere system, it is shown that GPS phase scintillation is primarily enhanced in the cusp, the tongue of ionization that is broken into patches drawn into the polar cap from the dayside storm-enhanced plasma density, and in the auroral oval. In this paper we examine the relation between the scintillation and auroral electrojet currents observed by arrays of ground-based magnetometers as well as energetic particle precipitation observed by the DMSP satellites. Equivalent ionospheric currents are obtained from ground magnetometer data using the spherical elementary currents systems technique that has been applied over the ground magnetometer networks in North America and North Europe. The GPS phase scintillation is mapped to the poleward side of strong westward electrojet and to the edge of the eastward electrojet region. Also, the scintillation was generally collocated with fluxes of energetic electron precipitation observed by DMSP satellites with the exception of a period of pulsating aurora when only very weak currents were observed.

We examine the role of the magnetosheath in solar wind-magnetosphere-ionosphere coupling using the Time History of Events and Macroscale Interactions during Substorms plasma and magnetic field observations in the magnetosheath together with OMNI solar wind data and auroral electrojet recordings from the International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer chain. We demonstrate that the electric field and Poynting flux reaching the magnetopause are not linear functions of the electric field and Poynting flux observed in the solar wind: the electric field and Poynting flux at the magnetopause during higher driving conditions are lower than those predicted from a linear function. We also show that the Poynting flux normal to the magnetopause is linearly correlated with the directly driven part of the auroral electrojets in the ionosphere. This indicates that the energy entering the magnetosphere in the form of the Poynting flux is directly responsible for driving the electrojets. Furthermore, we argue that the polar cap potential saturation discussed in the literature is associated with the way solar wind plasma gets processed during the bow shock crossing and motion within the magnetosheath.

We investigate the characteristics of rapid geomagnetic variations at high latitudes based on the occurrence of large time derivatives of the horizontal magnetic field (dH/dt exceeding 1 nT s−1). Analysis of IMAGE magnetometer data from North Europe in 1983–2010, covering more than two solar cycles, confirms and specifies several previous findings. We show that dH/dt activity is high around the midnight and early morning hours, and nearly vanishes at noon and early afternoon. This happens during all seasons, although the midnight maximum is nearly invisible during summer. As indicated by modelled ionospheric equivalent currents, large dH/dt values occur predominantly during westward ionospheric electrojets. Before and around midnight, dH/dt tends to be north-south oriented, whereas in the morning hours, its direction is more west-east directed. dH/dt tends to be more strictly north-south oriented during winter than other seasons. The seasonal occurrence of large dH/dt values is similar to the variation of the maximum amplitude of westward equivalent currents. The yearly fraction of east-west directed large dH/dt vectors at the Kilpisjärvi station (MLAT 65.88) varies from 31 to 47 % without any clear correlation with the general geomagnetic activity nor with the yearly averages of solar wind parameters.

The harnessing of the electric power associated with space weather storms is discussed. A simple system based on a large-scale coil located at high latitudes is devised and the power output of the system for a sample storm is calculated. It is shown that independently of the resistances of the conductors used in the coil and the number of turns in the system the maximum output of the system is restricted to the 1–10MW range, which is comparable to the maximum output of a single large wind turbine.

In this paper, we analyse in detail two famous space weather events; a railway problem on 13–14 July 1982 and a power blackout on 30 October 2003. Both occurred in Sweden during very intensive space weather storms and each of them a few years after the sunspot maximum. This paper provides a description of the conditions on the Sun and in the solar wind leading to the two GIC events on the ground. By applying modelling techniques introduced and developed in our previous paper, we also calculate the horizontal geoelectric field at the Earth's surface in southern Sweden during the two storms as well as GIC flowing in the southern Swedish 400 kV power grid during the event in October 2003. The results from the calculations agree with all measured data available. In the July-1982 storm, the geomagnetic field variation, ΔBx, reached values up to ~2500 nT/min and the geoelectric field reached values in the order of several volts per kilometer. In the October-2003 storm, the geomagnetic field fluctuations were smaller. However, GIC of some hundreds of amperes flowed in the power grid during the October-2003 event. Technological issues related to the railway signalling in July 1982 and to the power network equipment in October 2003 are also discussed.

In this work, source equivalence and computation of the re∞ected (induced) electromagnetic fleld in geophysical situations are studied. It is shown that the application of Huygens' principle allows for full generalization of Fukushima's equivalence theorem that applies only for magnetic fleld. The source equivalence is revisited for a vertical line current element, and it is shown that the equivalent charge required to replace the original source by a planar equivalent source together with the surface charge associated with the re∞ected fleld generates a purely vertical total electric fleld on the ground. Consequently, if the magnetic fleld and horizontal components of the total electric fleld on the ground are of interest, only equivalent currents need to be considered. The classical Complex Image Method (CIM) is derived from the exact image theory for planar impedance surfaces. The classical CIM is extended by considering a divergence-free source current that may have components also perpendicular to the ground plane. The extension is seen to generate a complex image charge not present in the classical CIM. Further, a generalized application of the extended CIM to geophysical situations having divergence-free volume source currents is introduced. The application involves decomposition of the source into linear current elements and rotations, translations and re∞ections of the electromagnetic fleld expressions associated with each element. The validity of the new approach is verifled for an

The Earth’s lithosphere and mantle respond to Space Weather through time-varying, depth-dependent induced magnetic and electric fields. Understanding the properties of these electromagnetic fields is a key consideration in modelling the hazard to technological systems from Space Weather. In this paper we review current understanding of these fields, in terms of regional and global-scale geology and geophysics. We highlight progress towards integrated European-scale models of geomagnetic and geoelectric fields, specifically for the purposes of modelling geomagnetically induced currents in power grids and pipelines.

In this work, methods to determine technological system parameters and the ground conductivity structure from different sets of geomagnetically induced current (GIC), magnetic field and geoelectric field observations are explored. The goal of the work is to enable optimal modeling of induced currents in any technological system experiencing GIC. As an additional product, the introduced methods can also be applied to utilize GIC observations in the imaging of the subsurface geological structures. Here a robust processing scheme and Occam’s inversion technique familiar from magnetotelluric (MT) studies are applied to the determination of the ground conductivity structure. The application of the methods to GIC data from the Finnish pipeline for a storm period of October 24-November 1, 2003 demonstrate that optimal system parameters and ground conductivity structure can be obtained using time series comprising only 8 days worth of data. Importantly, the obtained ground model is in agreement with models obtained in earlier MT studies. Furthermore, it is shown that although in an ideal case the magnetic field data used should be obtained from the immediate vicinity of the GIC observation site, some spatial separation (200–300 km) between the sites can be tolerated.

The goal of this study is to find a way to statistically estimate the Hall to Pedersen conductance ratio α from ground magnetic data. We use vector magnetic data from the CHAMP satellite to derive this relation. α is attained from magnetic satellite data using the 1-D Spherical Elementary Current Systems (SECS). The ionospheric equivalent current density can either be computed from ground or satellite magnetic data. Under the required 1-D assumption, these two approaches are shown to be equal, which leads to the advantage that the statistics are not restricted to areas covered by ground data. Unlike other methods, using magnetic satellite measurements to determine α ensures reliable data over long time sequences. The statistical study, comprising over 6000 passes between 55° and 76.5° northern geomagnetic latitude during 2001 and 2002, is carried out employing data from the CHAMP satellite. The data are binned according to activity and season. In agreement with earlier studies, values between 1 and 3 are typically found for α. Good compatibility is found, when α attained from CHAMP data is compared with EISCAT radar measurements. The results make it possible to estimate α from the east-west equivalent current density Jφ; [A/km]: α=2.07/(36.54/|Jφ|+1) for JφJφ+1) for Jφ0 (eastward). Using the same data, statistics of ionospheric and field-aligned current densities as a function of geomagnetic latitude and MLT are included. These are binned with respect to activity, season and IMF BZ and BY. For the first time, all three current density components are simultaneously studied this way on a comparable spatial scale. With increasing activity, the enhancement and the equatorward expansion of the electrojets and the R1 and R2 currents is observed, and in the nightside, possible indications of a Cowling channel appear. During southward IMF BZ, the electrojets and the R1 and R2 currents are stronger and clearer than during northward BZ. IMF BY affects the orientation of the pattern.

Motivated by the needs of the high-voltage power transmission industry, we use data from the high-latitude IMAGE magnetometer array to study characteristics of extreme geoelectric fields at regional scales. We use 10-s resolution data for years 1993–2013, and the fields are characterized using average horizontal geoelectric field amplitudes taken over station groups that span about 500-km distance. We show that geoelectric field structures associated with localized extremes at single stations can be greatly different from structures associated with regionally uniform geoelectric fields, which are well represented by spatial averages over single stations. Visual extrapolation and rigorous extreme value analysis of spatially averaged fields indicate that the expected range for 1-in-100-year extreme events are 3–8 V/km and 3.4–7.1 V/km, respectively. The Quebec reference ground model is used in the calculations.

Previous studies have demonstrated a close relationship between the time derivative of the horizontal geomagnetic field vector (dH /dt ) and geomagnetically induced currents (GIC) at a nearby location in a power grid. Similarly, a high correlation exists between GIC and the local horizontal geoelectric field (E ), typically modelled from a measured magnetic field.Considering GIC forecasting, it is not feasible to assume that detailed prediction of time series will be possible. Instead, other measures summarising the activity level over a given period are preferable. In this paper, we consider the 30-min maximum of dH /dt or E as a local activity indicator (|dH /dt |30 or |E |30 ). Concerning GIC, we use the sum of currents through the neutral leads at substations and apply its 30-min maximum as a regional activity measure (GIC30 ).We show that |dH /dt |30 at a single point yields a proxy for GIC activity in a larger region. A practical consequence is that if |dH /dt |30 can be predicted at some point then it is also possible to assess the expected GIC level in the surrounding area. As is also demonstrated, |E |30 and GIC30 depend linearly on |dH /dt |30 , so there is no saturation with increasing geomagnetic activity contrary to often used activity indices.

We present a new calculation method for solving inductive electric fields in the ionosphere. The time series of the potential part of the ionospheric electric field, together with the Hall and Pedersen conductances serves as the input to this method. The output is the time series of the induced rotational part of the ionospheric electric field. The calculation method works in the time-domain and can be used with non-uniform, time-dependent conductances. In addition, no particular symmetry requirements are imposed on the input potential electric field. The presented method makes use of special non-local vector basis functions called the Cartesian Elementary Current Systems (CECS). This vector basis offers a convenient way of representing curl-free and divergence-free parts of 2-dimensional vector fields and makes it possible to solve the induction problem using simple linear algebra. The new calculation method is validated by comparing it with previously published results for Alfvén wave reflection from a uniformly conducting ionosphere.

The method of 1D spherical elementary current systems (SECS) is a new way for determining ionospheric and field-aligned currents in spherical geometry from magnetic field measurements made by a low-orbit satellite. In contrast to earlier methods, the full ionospheric current distribution, including both divergence-free and curl-free horizontal currents, as well as field-aligned currents, can be determined. Placing infinitely many 2D SECSs of identical amplitudes at a constant latitude results in two types of 1D SECSs, which are independent of longitude, and by superposition can reproduce any ionospheric and field-aligned current system with the same property. One type of the 1D SECSs is divergence-free and toroidal with a poloidal magnetic field, and the other type is curl-free and poloidal. Associated with the divergence of the curl-free type are radial currents. The magnetic field of the combined curl-free 1D SECS and field-aligned currents is toroidal and restricted to the region above the ionosphere. Ionospheric currents are determined by placing several 1D SECSs at different latitudes and choosing their amplitudes in such a way that their combined magnetic field as closely as possible fits the one measured by the satellite. The 1D SECS method has been tested using both modeled and real data from the CHAMP satellite, and found to work excellently in 1D cases.

Rapid changes in the ionospheric current system give rise to induction currents in the conducting ground that can significantly contribute to magnetic and especially electric fields at the Earth's surface. Previous studies have concentrated on the surface fields, as they are important in, for example, interpreting magnetometer measurements or in the studies of the Earth's conductivity structure. In this paper we investigate the effects of induction fields at the ionospheric altitudes for several realistic ionospheric current models (Westward Travelling Surge, Ω-band, Giant Pulsation). Our main conclusions are: 1) The secondary electric field caused by the Earth's induction is relatively small at the ionospheric altitude, at most 0.4 mV/m or a few percent of the total electric field; 2) The primary induced field due to ionospheric self-induction is locally important, ~ a few mV/m, in some "hot spots", where the ionospheric conductivity is high and the total electric field is low. However, our approximate calculation only gives an upper estimate for the primary induced electric field; 3) The secondary magnetic field caused by the Earth's induction may significantly affect the magnetic measurements of low orbiting satellites. The secondary contribution from the Earth's currents is largest in the vertical component of the magnetic field, where it may be around 50% of the field caused by ionospheric currents. Keywords. Geomagnetism and paleomagnetism (geomagnetic induction) – Ionosphere (electric fields and currents)