Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Astron. Space Sci., 13 August 2025

Sec. Space Physics

Volume 12 - 2025 | https://doi.org/10.3389/fspas.2025.1610276

Localized geomagnetic disturbances: a statistical analysis of spatial scale

  • 1University of New Hampshire, Department of Physics and Astronomy, Durham, NH, United States
  • 2University of New Hampshire, Institute for the Study of Earth, Oceans, and Space, Durham, NH, United States
  • 3Departamento de Física, Universidad de Santiago de Chile, Santiago, Chile
  • 4Center for Interdisciplinary Research in Astrophysics and Space Sciences (CIRAS), Universidad de Santiago de Chile, Santiago, Chile
  • 5Geophysical Institute, University of Alaska Fairbanks, Fairbanks, AK, United States

Geomagnetically induced currents (GICs) pose a significant space weather hazard, driven by geomagnetic field variation due to the coupling of the solar wind to the magnetosphere-ionosphere system. Extensive research has been dedicated to understanding ground-level geomagnetic field perturbations as a GIC proxy. Still, the non-uniform aspect of geomagnetic fluctuations make it difficult to fully characterize the ground-level magnetic field across large regions of the globe. Here, we focus on localized geomagnetic disturbances (LGMDs) in the North American region and specify the degree to which these disturbances are localized. Employing the electrodynamics-informed Spherical Elementary Current Systems (SECS) method, we spatially interpolate magnetic field perturbations between ground-based magnetometer stations. In this way, we represent the ground magnetic field as a series of heatmaps at high temporal and spatial resolution. We leverage heatmaps from storm time during solar cycle 24 to automatically identify LGMDs. We build a statistical picture of the frequency with which LGMDs occur, their scale sizes, and their latitude-longitude aspect ratios. Additionally, we use an information theory approach to quantify the dependence of these three attributes on the phase of the solar cycle. We find no clear influence of the solar cycle on any of the three attributes. We offer some avenues toward explaining why LGMDs might behave broadly the same whether they arise during solar maximum or solar minimum.

1 Introduction

Detailed descriptions of geomagnetic field activity near the surface of the Earth are of acute interest to a variety of communities. Utility providers (including, but not limited to, electric power transmission companies, railway companies, and fluid pipeline operators) can wield information about ground-level geomagnetic activity as a means to mitigate disruption to their infrastructure by geomagnetically induced currents (GICs) (Molinski, 2002; Patterson et al., 2023; Boteler and Trichtchenko, 2015). It is common to derive an induced geoelectric field, and the corresponding GICs, based on the measured geomagnetic field where no direct GIC measurements are available (Kelbert, 2020). Additionally, scientists value the ground-level magnetic field as an indicator for phenomena in the atmosphere, the inner magnetosphere, and beyond (e.g., Russell et al. (2008), Madelaire et al. (2023)). The United States government maintains a National Space Weather Strategy and Action Plan to safeguard economic interests and national security against the threat of geomagnetic storms among other space weather effects (National Science and Technology Council, 2019; National Science and Technology Council, 2023).

Of particular importance to these communities is knowledge about exactly where on Earth the horizontal component of the geomagnetic field, BH, is experiencing (or will experience) a perturbation away from its baseline value. After all, disturbances to the geomagnetic field are often spatially localized, especially during periods of intense temporal variability (Dimmock et al., 2020), which is when the threat from GICs is most severe (Pirjola, 2000). Areas in which BH is markedly more perturbed from baseline than in the surroundings can be viewed as instances of this localization effect, and will hereafter be referred to as localized geomagnetic disturbances, or LGMDs, following the terminology used by Vandegriff et al. (2024). In this study, we focus on LGMDs in BN, the north-south vector component of BH. A precise characterization of an LGMD’s spatial extent could help explain, for example, why GICs might have occurred in one region’s power network but not in that of another.

Several explanations for the localization effect have been proposed. It is known that the nonuniform electric conductivity of Earth’s crust plays a role in determining the deflection of all components of the geomagnetic field, including BN, implying that localized magnetotelluric structures can help produce LGMDs (Kelbert, 2020). Substorm activity also contributes to localized magnetic and electric field enhancements by way of intensification of the westward auroral electrojet and current structures therein (Viljanen et al., 2001; Ngwira et al., 2015; 2018). Additionally, geomagnetic pulsations, such as those caused by Kelvin-Helmholtz instabilities in outer-magnetosphere plasma flow, have been shown to produce LGMDs during geomagnetic storms even in mid-latitude regions (Heyns et al., 2021). Depending on the space weather conditions and the geographic location on Earth, these phenomena and others may vary in the magnitude of their contribution to the generation of LGMDs. Though the full suite of LGMD drivers has yet to be totally and precisely characterized, the contributing factors mentioned above represent rapid progress in understanding this complicated localization effect.

Some studies have investigated LGMD spatial scale by comparing readings at multiple stations, with attention paid to the geographic distance between those stations (Boteler and van Beek, 1999; Pulkkinen et al., 2015). To give some recent examples, Dimmock et al. (2020) found that the time derivative of the horizontal component of the magnetic field could be as much as three times higher at a single magnetometer station compared to the average of all magnetometers situated across a relatively small (500 km) region of Fennoscandia. They also found that nonuniformities in the geomagnetic field at the 200 km scale can be responsible for a large percentage of the strength of GICs. Similarly, Dimitrakoudis et al. (2022) examined the BEAR Fennoscandian magnetometer network and the CARISMA North American network to determine that the correlation scale is only about 500 km between magnetometers located at similar latitudes, and this scale drops to a mere 100 km in the vicinity of the auroral oval.

However, single-point measurements from magnetometers are insufficient to fully describe the size and shape of LGMDs because the distance between magnetometer stations may be larger than the scale sizes of the LGMDs themselves. That is to say, too sparse a distribution of magnetometer stations will not allow size scales to be precisely measured. There may be a few points of measurement, but the magnetic field’s behavior between those points is unknown. To overcome this deficiency, one requires a method by which to estimate BN (or its time derivative dBN/dt) at any point of interest regardless of the presence of a nearby magnetometer station or of an LGMD; in other words, a robust interpolation method is needed. Simply using the measurement of BN from the nearest magnetometer is generally not sufficient because it could differ greatly from the actual value at the point of interest, depending on the characteristics of the surrounding geology (Kelbert, 2020) and magnetospheric conditions (Viljanen, 1998). Even using a proximity-weighted average from multiple nearby magnetometers’ readings may result in notable estimation errors in regions where magnetometers are far apart from one another (Campbell, 2017). Applying the Spherical Elementary Current Systems (SECS) interpolation method addresses this need effectively. While SECS is complicated compared to several geometric methods, its electrodynamics-based approach delivers enhanced accuracy, particularly in magnetically disturbed conditions (McLay and Beggan, 2010). Interpolation using spherical harmonic analysis – for instance, as in Upendran et al. (2022) – is a competitive alternative to SECS; however, spherical harmonic analysis is more susceptible to interpolation artifacts (Rigler et al., 2019).

The goal of this investigation is to create a dataset of LGMDs on the North American region and to characterize their sizes and shapes during intervals of enhanced activity caused by geomagnetic storms. We use SECS to generate maps of geomagnetic field perturbations with high spatial resolution from a large cache of historical magnetometer data spanning storm times from 2009 through 2019. We then analyze the LGMDs that are displayed in these maps to reveal details about the frequency with which they occur, their spatial extent, and their longitude-latitude aspect ratios.

Finally, we quantify the solar cycle dependence of our LGMD dataset based on the Kullback-Leibler divergence, an important concept in information theory. Information theoretical approaches can be powerful tools when applied to space physics and space weather. Aryan et al. (2014) used the Kullback-Leibler divergence to quantify the relative importance of various solar wind parameters in the generation of chorus waves in the magnetosphere. Wing et al. (2018) studied the memory of the solar cycle and Snelling et al. (2020) studied the causal relationship between solar flares using transfer entropy and mutual information, respectively, two quantities which are closely related to the Kullback-Leibler divergence. Readers interested in further uses of information theory in space science should refer to Balasis et al. (2023), McGranaghan (2024), and references therein.

2 Data and data preparation

Obtaining the most complete picture of LGMD characteristics requires simultaneous data from multiple magnetometers. These data must describe long enough time scales so as to capture any variation due to the solar cycle. They must also exist at a fine enough temporal resolution to make clear how LGMD attributes evolve during periods of high geomagnetic activity. To address these requirements, we use 1-min cadence baseline-subtracted ground magnetic field data from the SuperMAG collaboration (Gjerloev, 2012). Our analysis in this work specifically examines the northward component of the perturbation vector of the interpolated magnetic field. We choose a single vector component because the northward (dBN) and eastward components do not necessarily behave the same (Viljanen and Tanskanen, 2011), and a thorough comparison between components is best reserved for future studies.

This study examines LGMD scale size over a particular region of interest (though it could, in principle, be replicated to study LGMDs across any given section of Earth’s surface). We choose the region of interest to be within the quadrangle which is aligned with the standard geographic coordinate system, and which extends from 45°N to 55°N and from 80°W to 130°W, as indicated by the orange region in Figure 1. Geographic coordinates are used here to be consistent with the use of geographic coordinates in Section 3, and offer no substantial change to the region size, shape, or data availability compared to geomagnetic coordinates. This particular region was chosen because it is of interest to contemporary GIC research: it lies at middle latitudes, is a geographically large region with good magnetometer coverage, and constitutes much of North America, which hosts a wide variety of subterranean magnetotelluric structures that influence the surface geomagnetic field. Such structures have been studied extensively, such as by Kelbert et al. (2019), and mapped as part of the National Science Foundation’s EarthScope project (Schultz, 2009). Geomagnetic storms have created GIC-related disruptions in mid-latitude North America, including the 2003 “Halloween storm” (minimum SYM-H index 432 nT) which markedly affected power systems (Kappenman, 2005; Piccinelli and Krausmann, 2014) and the March 1989 storm (minimum SYM-H 710 nT) which caused widespread blackouts in the northeastern United States and eastern Canada (Allen et al., 1989; Bolduc, 2002; Love et al., 2022). To the west of the Atlantic coast of North America, where the conductivity of the crust tends to be higher, the March 1989 storm still created hazardous ground-level geoelectric fields (Wei et al., 2013). Even storms weaker than the March 1989 storm have caused GICs in this western part of the region (e.g., Anderson et al. (1974); Cordell et al. (2024), the latter analyzing an event which occurred as recently as 2023).

Figure 1
Map of North America, indicating geographic latitude and longitude. Three concentric quadrangle regions are shown. Magnetometer data comes only from stations within the outermost blue quadrangle.. Inside it is a green quadrangle showing the geographic extent of the Spherical Elementary Current System grid. The innermost quadrangle is orange, representing the region over which analysis is performed.

Figure 1. The magnetometers whose data are used to reconstruct LGMDs within mid-latitude North America. We only use magnetometers within the area indicated by the blue quadrangle. The green quadrangle shows the geographic extent of the Spherical Elementary Current System grid used to model magnetic perturbations, and the orange quadrangle represents the region over which our analysis is performed. Lines of geographic latitude and longitude are shown in light gray.

We use only data from magnetometers that are within the larger region extending from 35°N to 70°N and from 60°W to 150°W (also in geographic coordinates), shown as the blue quadrangle in Figure 1. This allows the region of interest to be surrounded by magnetometer stations, a constraint that benefits the SECS method by reducing the magnitude of interpolation error at the edges. We further restrict the magnetometers we use based on data availability. We only use those with more than 25% of the data available each year, and less than 35% missing data averaged across all years, from 2009 through 2019. (The final selection of magnetometers is listed in Supplementary Table S1). After this step, our dataset consists of measurements from 25 magnetometers whose locations are marked in Figure 1.

We then discard timesteps wherein any magnetometer is missing a data point. Our SECS calculations described in Section 3 are only carried out when data from every magnetometer is available, which ensures consistency across timesteps with regard to SECS inversion. Therefore, we make this data availability restriction in order to preserve as many data points as possible while maximizing the number of magnetometers available for analysis.

Our choice of years across which our dataset spans – 2009 through 2019 – represents one full solar cycle (Solar Cycle 24). This allows all phases of the solar cycle to be captured, which is important because geomagnetic activity, as well as the magnetospheric processes that cause it, depend on the current solar cycle phase. We have chosen Solar Cycle 24 in particular because it is the most recent complete cycle and has the most SuperMAG stations with 1-min data available.

Finally, we examine only times during which geomagnetic storms occur. We use the same list of storms as Pinto et al. (2022), who define a storm as the period where the SYM-H index decreases below −50 nT and remains below this threshold for at least 2 h, plus 12 h of lead time before the SYM-H minimum and 24 h of recovery time afterward. Thus our final magnetometer dataset is composed of 183 storms and a total of 395,463 data points (equivalent to a total of 274.6 days). This is about 7% of the 5,784,480 data points that are in the full 2009–2019 period. We reduce our dataset in this way because we wish to study LGMDs during periods of enhanced geomagnetic activity.

3 Methods

3.1 Spherical elementary current systems

As mentioned before, extending BN measurements to regions without magnetometers can be a difficult task due to the influence of localized phenomena throughout the magnetosphere and the lithosphere. This study seeks to overcome such difficulties and precisely reconstruct geomagnetic field variations by interpolating with the Spherical Elementary Current Systems (SECS) technique (Amm and Viljanen, 1999). SECS is commonly used to estimate geoelectric fields for the purpose of calculating GICs over regional-scale infrastructure (e.g., Wik et al. (2008), Kelbert and Lucas (2020)). A summary of the technique is presented in this section; readers looking for a more in-depth explanation are encouraged to refer to Vanhamäki and Juusola (2020).

SECS interpolation operates by representing the ionospheric current distribution with a parametric model. Empirical determination of the model’s parameters at each moment in time is carried out based on magnetometer observations, which allows for the assumption that the model current distribution causes all observed ground-level geomagnetic activity. Once the model is fitted using this approach, one can apply the Biot-Savart law to the model current distribution in order to calculate the magnetic field deflection at any point within a region of interest.

The model current distribution is a superposition of curl-free and divergence-free vector functions defined on a common spherical surface, each with a different location for its pole. These basis functions are called elementary current systems. One is free to choose the placement of the poles of elementary current systems, but typically they are arranged in a latitude-longitude grid for convenience. By Helmholtz’s Theorem, any smooth, differentiable vector function can be exactly decomposed into such basis functions; using a finite number of basis functions at worst only approximates the original vector field. The divergence-free part (Jdf) and curl-free part (Jcf) of each elementary current system are defined as (Equation 1):

Jdf(ϕ,θ)=I0,df4πRIcotθ2ϕ̂,Jcf(ϕ,θ)=I0,cf4πRIcotθ2θ̂(1)

One would typically consider the coordinate system in which r is the radial distance from the center of Earth, and ϕ and θ are geographic longitude and colatitude, respectively. The “primed” spherical coordinate system (r,ϕ,θ) differs from the standard “unprimed” (r,ϕ,θ) system in that its north pole, where θ=0, is defined to be at the location of that current system. The azimuthal angle ϕ is defined accordingly. As opposed to the geographic unprimed coordinate system, we will hereafter call a primed coordinate system an egocentric coordinate system, since the pole of each current system has its own unique primed coordinate system centered on itself.

Here the currents are assumed to lie in an infinitely thin current shell of radius RI (which typically takes the value RE + 110 km, where RE is the radius of Earth, so as to approximately correspond with the heights of the maximum Pedersen and Hall conductivities). The factors I0,df and I0,cf are empirically determined at every moment in time in order to fit this model ionosphere to magnetometer data. The total vector field representing ionospheric currents may not actually equal the currents in the real ionosphere, especially if the associated field-aligned currents are not approximately radial. However, the elementary currents still form a basis from which the magnetic field can be later reconstructed.

Under the assumptions of this current-shell model, it can be shown that the ground magnetic effect of the curl-free current systems vanishes within the current surface (Fukushima, 1976). This allows I0,df to be rewritten more simply as I0. Thus, the horizontal component of the magnetic field due to a single current system at a point r=(RE,ϕ,θ) is solely due to Jdf and can be shown to be (Equation 2):

bH(θ)=μ0I04πREsinθRERIcosθ12REcosθRI+RERI2+cosθθ̂(2)

The constant μ0 in this equation denotes the magnetic permeability of free space. There is also a radial component to this magnetic field, but it will not be used in this analysis because we are only concerned with horizontal field components (specifically, this study eventually only examines the northward component). Then, the total horizontal vector magnetic field BH(θ) is the sum of the magnetic contributions bH(θ) from each current system.

To determine the contributions from each current system, one can define a vector Z containing the ground-truth magnetometer measurements at stations 1 through n, a vector I of scaling factors for current systems 1 through m, and a 2n×m transfer matrix T̄ (Equation 3):

Z=z1,θz1,ϕz2,θz2,ϕzn,θzn,ϕ,I=I0,1I0,2I0,m,T̄=T1,1,θT1,2,θT1,m,θT1,1,ϕT1,2,ϕT1,m,ϕT2,1,θT2,2,θT2,m,θT2,1,ϕT2,2,ϕT2,m,ϕTn,1,θTn,2,θTn,m,θTn,1,ϕTn,2,ϕTn,m,ϕ(3)

Each element Ti,j,(θ,ϕ) represents the θ or ϕ component of the magnetic effect that current system j has on magnetometer i. One must then solve the matrix equation Z=T̄I for I and plug the resulting scaling factors into Equation 2 to obtain interpolated magnetic field values. This is typically done using singular value decomposition and setting small singular values to 0 to properly condition the matrix. One usually considers “small singular values” to mean those that are less than the maximum singular value multiplied by some cutoff coefficient ϵ. The best value of ϵ depends on the specifics of the SECS model in question; refer to Section 3.2 for information on how this value was determined for the present case.

When using Equation 2 to find the contribution of each current system to the magnetic field, it is important to notice again the difference between the geographic and the egocentric coordinate systems. These contributions are most easily added together when recognizing that a contribution vector bego=(bθ,bϕ) at geographic colatitude and longitude (θi,ϕi) due to a current system at geographic colatitude and longitude (θj,ϕj) may be transformed from the egocentric to the geographic coordinate system with Equations 4, 5:

θego=arccoscosθicosθj+sinθisinθjcos(ϕiϕj),(4)
bH=bθbϕ=bNbE=bθcoscosθjcosθicosθegosinθisinθegobθsinsinθjsin(ϕjϕi)sinθego(5)

The total instantaneous magnetic field in the geographic coordinate system is then the sum of each contribution vector bH, which we calculate component-wise to find the northward component, BN, and the eastward component, BE. Since BN is this study’s variable of interest, we focus solely on BN without using BE.

3.2 Model details

We define our model current distribution as the superposition of a set of current systems whose poles are placed on a quadrangular latitude-longitude grid extending from 40.5°N to 60.5°N and from 70.5°W to 140.5°W (the green quadrangle in Figure 1) in geographic coordinates. The current system array is offset half a degree from whole-number placement to avoid colocating the current systems with points at which the magnetic field is evaluated, as this would cause a singularity to arise in Equation 2. There is also the danger of a singularity if a current system is colocated with a magnetometer; however, this becomes obvious with a cursory inspection of the interpolation outputs. If one places the current systems such that this is not a problem, one can rest assured that no singularities will arise later because the locations of current systems and magnetometers do not change over time.

Poles are placed every 2°, making this a 10-by-35 array of current systems. The performance of SECS interpolation is not especially sensitive to small changes to the number and spacing of the current systems, so the 2-degree separation was chosen based on its use in McLay and Beggan (2010). We further validated this choice by performing joint Bayesian optimization of the number of rows and columns of current systems (while keeping the boundaries of the array fixed) using the Optuna package for Python (Akiba et al., 2019). The objective of this optimization was to minimize the absolute error between the model outputs at the location of each magnetometer and the measurements of the magnetometers themselves, averaged across magnetometers in each timestep, and then averaged across all timesteps in the dataset. The errors at each magnetometer were calculated using a leave-one-out scheme in which the model was fit using data from all magnetometers except the one in question, meaning a new model was fit for each magnetometer, all identical except for which magnetometer was omitted. The optimization found that the error was minimized for an 8-by-39 array of current systems. This agrees with our similar choice of array dimensions, whose error is less than 1% higher than the optimum. As such, we kept the 10-by-35 array so as to maintain a rounded 2-degree separation, as well as for possible future comparison to studies such as McLay and Beggan (2010).

Together with the number of rows and columns of current systems, the above procedure also jointly optimized a third quantity: the cutoff coefficient ϵ. We use the best value determined through this optimization, ϵ=0.0932. This value is comparable in magnitude to what has been used in previous related studies (McLay and Beggan, 2010; Weygand et al., 2011).

When creating a SECS-based model, it is sometimes helpful to introduce an additional current sheet located at a negative altitude (that is, internal to the planet). Doing so may allow one to delineate which part of geomagnetic perturbations arise from ionosphere-associated phenomena and which perturbations arise from currents induced in the subsurface (Pulkkinen et al., 2003). However, the goal of the SECS model as it is used here is only to reproduce geomagnetic perturbations, not to explain them in terms of their contributing equivalent currents. As a result, the additional current sheet has not been employed in this study.

3.3 Identification and characterization of LGMDs

Measurements of dBN are interpolated to give estimates of this field component at a set of points within the region of interest. The points could have any spatial density since they are sampled from the output of the SECS method, which is a continuous function. We limit ourselves to making estimations on a regularly spaced 14-by-32 grid of points to keep computation time under control. Each of these “heatmaps” (e.g., Figure 2) is associated with a set of simultaneous magnetometer measurements, so exactly one heatmap is created per timestep in our dataset.

Figure 2
Map displaying magnetic field interpolations over North America, from an example heatmap. The color ranges from blue, representing southward perturbations, to red, representing northward perturbations. The scale on the colorbar to the right shows values from -40 to 100 nanoteslas. Three black contours are shown on the map to denote the perimeters where LGMDs have been found.

Figure 2. An example heatmap showing the result of interpolating dBN measurements across the region of interest. Automatically drawn contours that surround areas of intense disturbance are shown heavy black lines. Contours that are closed are considered to demarcate LGMDs. The region on the left is not considered an LGMD because it cannot be enclosed by a contour, whereas the middle and right regions are considered LGMDs.

Once the heatmaps have been created, one can visually inspect them to view the estimated extent of LGMDs. However, we need to define a method to locate LGMD in a given heatmap. We do so here by granting a quantitative aspect to the definition of LGMD given in Section 1; that is, by defining a threshold value dBthreshold. Any point in a heatmap whose absolute value exceeds that of the threshold (i.e., |dBN|>dBthreshold), is considered part of an LGMD. Points identified this way that are adjacent, including diagonally, are considered to be part of the same LGMD. We then draw a contour line around each LGMD in a heatmap using the marching squares contour-finding algorithm implemented in the measure.find_contours function of the scikit-image Python package (Van Der Walt et al., 2014).

To identify LGMDs in this study, we define our perturbation threshold as dBthreshold = 25.95 nT, which corresponds to the absolute value of the median of the distribution of the minimum values from each heatmap. It is not a particularly strict threshold, but that is because we are interested in all LGMDs, not only anomalously intense ones. LGMDs that extend past the edge of the region of interest, and thus are not associated with a closed contour, cannot be characterized with confidence, and are thus completely omitted from our analysis. Examples of both closed and open contours drawn with our method can be seen in Figure 2.

All our results in Section 4 are based on the heatmaps created using the procedure outlined in this section. Some are, more specifically, based on the perimeters of LGMDs identified in the heatmaps, each of which is calculated by adding up the distances between successive vertices of the contour that surrounds it. One must respect the spheroidal shape of Earth when computing these distances. Given two contour vertices with (latitude, longitude) coordinates (λ1,ϕ1) and (λ2,ϕ2) in radians and the radius of Earth RE = 6,378 km, we find the haversine distance d between vertices (Equation 6) rather than the Euclidean distance as

d=2REarcsinsin2λ2λ12+cosλ1cosλ2sin2ϕ2ϕ12.(6)

In this work, we focus on perimeter as a measure of LGMD size. While there is no fundamental problem with investigating the behavior of LGMDs with different measures – for example, diameter – concentrating on the perimeter of LGMDs removes a degree of arbitrariness from the study. In the example of studying LGMD diameter, one would be forced to confront the question of along which direction to calculate that diameter, especially considering that LGMDs can vary widely in shape (an effect which is displayed in Figure 4). This is not an insurmountable issue; we hope that future work will extend the simpler baseline focus of this study to include other size metrics.

We also wish to study the shape of LGMDs, and pursue analysis of their aspect ratios as a way of doing so. We define the aspect ratio A for a given LGMD as (Equation 7):

A=(ϕEϕW)cosλE+λW2λNλS,(7)

with λ denoting the latitude and ϕ denoting the longitude of the easternmost, westernmost, northernmost, and southernmost points (subscripts E, W, N, and S, respectively) on the LGMD’s perimeter. This quantity describes the ratio of the LGMD’s longitudinal extent to its latitudinal extent, with the longitudinal extent being measured along the line of latitude halfway between the easternmost and westernmost points. Section 4 contains the results of the application of this definition to our set of computed contours.

4 Results and discussion

4.1 Full-solar-cycle distributions

After applying Spherical Elementary Current Systems and our contour-finding algorithm, we can study the behavior of LGMDs by examining relevant statistics of the resultant LGMD perimeter dataset. For that, we first consider the distribution of the number of LGMDs identified in each of our 395,463 heatmaps shown in Figure 3a. The distribution has a significant skew to the right, and no heatmap was ever observed to contain more than six LGMDs. A large proportion – about 58% – of heatmaps exhibit no LGMDs, though the exact amount is somewhat sensitive to the dBN threshold used to discriminate LGMDs from background geomagnetic activity.

Figure 3
Two histogram charts show data on identified Localized Geomagnetic Disturbances (LGMDs). Chart a displays the frequency of occurrence by number of LGMDs, decreasing from over 100,000 for zero LGMDs to just under 1,000 for six LGMDs. Chart b shows frequency by LGMD perimeter in kilometers, with higher frequencies at lower perimeters, decreasing steadily and with minor fluctuations as perimeter increases to 12,000 kilometers.

Figure 3. (a) The distribution of how frequently LGMDs are present in dBN heatmaps as a function of how many separate LGMDs are indeed present. (b) The distribution of the perimeters of all LGMDs over the entire solar cycle.

We can also consider the perimeters of the LGMDs, as shown in Figure 3b. This figure describes the overall distribution of perimeter values. It is agnostic toward simultaneity and simply includes all LGMDs; that is, any two LGMDs contribute their perimeters toward this distribution even if they reside in the same heatmap. This distribution has a significant positive skew and shows that the mean, median, and mode LGMD perimeters are 1,620 km, 1,372 km, and 550 km respectively. Also of note is that the largest LGMD was observed to be 12,024 km in perimeter. As a point of comparison, the region of interest measures 9,344 km in perimeter, but it is possible for an LGMD perimeter to exceed this value if its associated contour forms a concave shape. However, most LGMDs are smaller than 9,000 km in perimeter, with only a handful of cases in the entire dataset close to the size of the region of interest.

Importantly, the typical perimeter of LGMDs as per these results is of the order that one would expect from estimates of the few-hundred-kilometer spatial scale derived using other methods, such as those mentioned in Section 1 (e.g., Dimmock et al. (2020); Dimitrakoudis et al. (2022)). This approximate correspondence between perimeter and diameter holds unless the LGMD in question is extremely oblong, which is quite rare according to the results contained in Figure 4.

Figure 4
Histogram of LGMD aspect ratios. The x-axis represents the logarithm base ten of aspect ratios ranging from approximately -0.5 to 1.25. The y-axis shows the frequency of occurrence on a logarithmic scale from 10^0 to 10^4. The distribution peaks around 0.25 and decreases log-symmetrically towards both ends.

Figure 4. The distribution of the aspect ratios of all LGMDs over the entire solar cycle. Note that the distribution is centered on a positive number, meaning that LGMDs tend to be wider longitudinally than they are latitudinally.

Next, we examine the distribution of the aspect ratios of observed LGMDs. The distribution shown in Figure 4 is roughly symmetric about its central value at log(A) = 0.29 (A = 1.95). This value is markedly positive, meaning that we observe LGMDs with a longitudinal extent greater than their latitudinal extent more frequently than those with the opposite orientation. In fact, only 5.7% of the mass of the log aspect ratio distribution lies below zero. This means that over 94% of observed LGMDs are at least slightly oblong along the east-west direction. The cause of this asymmetry around log(A) = 0 is not entirely clear, but it may be related to the fact that current systems which induce perturbations to BN, such as the auroral electrojets, are also extended longitudinally.

One must note that our choice of the region of interest could be biasing this distribution in the positive direction by “cutting off” some data points on the negative end of the distribution. This is because the region of interest is wider in the longitudinal direction than it is in the latitudinal direction. To illustrate this, an LGMD which has been found to have a longitudinal extent of 1,400 km might be counted among those in Figure 4, but one with a latitudinal extent of 1,400 km would surpass the boundaries of the region of interest; its perimeter contours would not be closed, and thus it would not be represented in Figure 4. If this bias does exist, it could be corrected by extending the region of interest to have the same latitudinal extent as its longitudinal extent. However, for our purposes, this is not practical, because there are not sufficient magnetometers with available data to support such an expansion without reducing interpolation quality. Instead, we examined the distribution of aspect ratios of only LGMDs less than 1,000 km in perimeter. The LGMDs in this subset are too small to be cut off by the boundaries of the region of interest, even if they are oblong in the north-south direction. When normalized, the distribution of this subset is almost identical to the full distribution. We take this to be sufficient evidence that the dimensions of the chosen region of interest do not significantly impact the distribution of LGMD aspect ratios, at least for statistical purposes.

We have inspected the distributions of three attributes: the number of LGMDs per heatmap, the perimeter of LGMDs, and the aspect ratio of LGMDs. For readers interested in additional details concerning these LGMD attributes, plots of the distributions of perimeter and aspect ratio, separated by the number of simultaneous LGMDs, are available in the Supplementary Material. Additional statistics, whose significance is not analyzed in this work, are also available in the Supplementary Figure S3S8.

4.2 Impact of solar cycle phase

We are also interested in whether the phase of the solar cycle has any relationship to the statistics explored in the previous sections. Therefore, we now proceed with an analysis of two different subsets of our dataset: the set of timesteps surrounding solar maximum (14 October 2012 through 31 January 2015) and those surrounding solar minimum (1 January 2009 through 31 October 2009 and 4 March 2018 through 31 December 2019) for Solar Cycle 24. The dates defining these intervals are set following Reyes et al. (2021), who precisely defined the phase boundaries for that solar cycle.

In this section, we reproduce each histogram from the previous section to include the solar maximum and solar minimum data subsets alongside the full dataset. These subsets each contain different numbers of data points because each represents a different range of instances in time. Therefore, the distributions shown have been normalized such that the area under each is 1. This is in contrast to the distributions in Section 4.1; those essentially represent the relevant attributes’ probability distributions, whereas the plots in this section show probability density functions (PDFs) because they have been normalized in this way. We show PDFs so that the data subsets may be directly compared. The distribution of the number of LGMDs per heatmap, though it represents a discrete variable as opposed to a continuous variable, can be represented by this normalization as a PDF because the “bin” width is 1.

The effects of separating solar minimum and solar maximum from the overall dataset are shown on the number of LGMDs (Figure 5), their perimeter (Figure 6), and their longitudinal-latitudinal aspect ratios (Figure 7). For each of these three statistics, the trends of the solar maximum and solar minimum subsets follow each other very closely. This motivates the question: are the differences that do exist between the subsets simply due to aleatoric variability (that is, do they only exist by chance?), or are they caused by some higher-order physical phenomena? To answer this question, we can compute a statistical divergence between the distributions corresponding to the solar maximum and solar minimum subsets. If this divergence is sufficiently large, that will be evidence against the hypothesis that the distributions differ only because of aleatoric effects.

Figure 5
Bar chart titled

Figure 5. The blue bars represent the discrete probability density function (PDF) of the number of LGMDs per 1-min timestep over the entirety of Solar Cycle 24 and are equivalent to the histogram shown in Figure 3a. The orange and green bars show the same statistic for the solar minimum and solar maximum data subsets, respectively.

Figure 6
Histogram depicting LGMD Sizes showing probability density versus LGMD perimeter in kilometers. Three distributions are shown: full solar cycle (blue), solar minimum (orange), and solar maximum (green). The probability density decreases with increasing perimeter across all three distributions.

Figure 6. The blue trace represents the probability density function (PDF) of LGMD perimeter over the entirety of Solar Cycle 24 and is equivalent to the histogram shown in Figure 3b. The orange and green traces show the same statistic for the solar maximum and solar minimum data subsets, respectively.

Figure 7
Histogram showing LGMD aspect ratios with probability density on the y-axis and log10 of aspect ratio on the x-axis. Three datasets are compared: Full Solar Cycle (blue), Solar Minimum (orange), and Solar Maximum (green). Each dataset shows distinct distributions of aspect ratios.

Figure 7. The blue trace represents the probability density function (PDF) of LGMD aspect ratio over the entirety of Solar Cycle 24 and is equivalent to the histogram shown in Figure 4. The orange and green traces show the same statistic for the solar minimum and solar maximum data subsets, respectively.

We first quantify the distance between a PDF associated with the solar maximum dataset, X, and that associated with the solar minimum dataset, N, using JKL, the symmetrized form of the Kullback-Leibler divergence DKL (Kullback and Leibler, 1951) (Equation 8):

JKL(X,N)=DKL(X,N)+DKL(N,X)=iwXilnXiNi+iwNilnNiXi(8)

Here, the sums are over the “bins” of the discretized PDFs. The width of each bin is w, which is written without subscript because, for this study, each bin within the same PDF has the same width. We use the symmetrized form of the divergence to avoid confusion that may arise due to DKL(X,N) generally not being equal to DKL(N,X).

From an information theory perspective, JKL measures how much information would be absent if one PDF were to be used to approximate the other. As the distance between X and N increases, JKL(X,N) will also increase. However, there is some ambiguity as to the threshold above which JKL can be said to indicate that the differences between X and N are due to more than just aleatoric variability. Some clarification appears as one takes JKL to various limits. First, the symmetrized divergence between a distribution and itself is zero. Second, while there is no upper bound on JKL even for PDFs, one can create a “benchmark” PDF U known a priori to be sufficiently distant from X. If JKL(X,N) is much closer to JKL(X,U) than it is to 0, that may point to a phenomenological difference between X and N.

For this comparison, we simply define U as the PDF describing a normalized uniform distribution, which is fundamentally different from any of the PDFs shown in Figures 57. The JKL(X,N) values for each of the statistics presented in this section are shown in Table 1. Also shown for comparison are the symmetrized divergences between X and U and between N and U for LGMD number, perimeter, and the base-10 logarithm of aspect ratio, log10(A). Because any conclusion drawn about log10(A) has a linear equivalent that can be drawn about A, we will refer to the “base-10 log aspect ratio” just as the “aspect ratio” for the sake of simplicity.

Table 1
www.frontiersin.org

Table 1. Values for the symmetrized Kullback-Leibler divergence between various probability density functions (PDFs) formed from data from solar maximum (X) and solar minimum (N). Divergences from U, the PDF of a uniform distribution, are shown for comparison.

It is clear from the reported values that the PDFs of each statistic shows much less divergence between solar maximum and solar minimum than they do between a uniform PDF and either solar maximum and solar minimum. The values of JKL(X,N) are considerably larger for perimeter and aspect ratio than they are for LGMD number, but this likely due to the larger number of bins in those two quantities’ PDFs accentuating the divergences between PDFs; JKL(X,N) for perimeter and aspect ratio is in each case more than an order of magnitude smaller than either JKL(X,U) or JKL(N,U). The fact that JKL(X,N) is so small in comparison does not suggest than these three LGMD attributes exhibit any dependence on the phase of the solar cycle.

The information theory perspective makes the JKL criterion interpretable. However, there are some disadvantages to this approach. The effect of the logarithms in the definition of JKL (Equation 8) emphasizes the relative importance of the tails of the PDF to the symmetrized divergence. This could raise the concern that outliers in each of the attributes’ distributions may be influencing the conclusions drawn from Table 1. Additionally, bins in which one of the PDFs is equal to zero must be ignored to keep JKL well defined. This means that some small number of bins in the perimeter and aspect ratio PDFs are not taken into account by this criterion. These drawbacks do not invalidate the use of JKL in this context, but they do suggest another criterion may be needed to verify the preliminary conclusions about the effect of the solar cycle. We have used the total variation distance DTV (Equation 9) to supplement the use of JKL:

DTV(X,N)=12iwXiNi(9)

The total variation is computed outside of “logarithm space”, which avoids both drawbacks mentioned regarding JKL. The DTV values for each of the statistics presented in this section are shown in Table 2. The divergence between X and N remains the smallest divergence by over one order of magnitude for the number of LGMDs. For their perimeter and their aspect ratio, it is smaller by about a factor of eight. Using the total variation distance as a criterion reinforces the results obtained with the initial Kullback-Leibler analysis, which indicate insufficient evidence to suggest a phenomenological difference between solar maximum and solar minimum LGMD aspect ratios.

Table 2
www.frontiersin.org

Table 2. Values for the total variation distance between various probability distributions formed from data from solar maximum (X) and solar minimum (N). Distances from U, the PDF of a uniform distribution, are shown for comparison.

One might expect the solar cycle to have an effect on the PDFs of these three LGMD attributes. After all, the solar cycle has an effect on the general geoeffectiveness of the solar wind–strong geomagnetic storms, storms due to coronal mass ejections, and strong ground-level geomagnetic activity are all more common during solar maximum, for instance (Webb, 1991; Richardson and Cane, 2012; Borovsky and Denton, 2006). However, as we move away from global size scales toward the smaller scales of LGMDs, we can not reject the hypothesis that the solar maximum and solar minimum distributions differ only because of aleatoric effects.

This does not necessarily mean that the solar cycle has no influence on the properties of LGMDs, only that any influence that does exist could not be detected with the methods used here. This is perhaps due to the fairly strict choice of the uniform distribution’s PDF as a benchmark PDF. One distribution being more similar to a uniform distribution than to the other solar cycle phase’s distribution – the criterion necessary to reject the baseline hypothesis – is a high bar to meet. There still may be some much smaller solar cycle effect that could require a more sensitive test to extract. We have used the uniform distribution’s PDF as a benchmark even still because other tests would include further amounts of arbitrariness. For instance, one could repeat this analysis on LGMD aspect ratio, replacing U with the PDF of a log-normal distribution, but it would be difficult to justify a choice of mean and variance for this new benchmark. Because of this, first investigating the influence that other phenomena have on LGMD attributes would likely be more fruitful than identifying a more sensitive test.

Though this numerical analysis suggests that the distribution of aspect ratios does not strongly depend on the phase of the solar cycle, it does not reveal the physical processes that do determine the shape of the PDFs shown. This motivates a cause-and-effect analysis of LGMD aspect ratio, which likely deserves a study all to itself. Such a study could explore phenomena inside the magnetosphere, where the effect of the solar wind is modulated by a multitude of current systems, and it could connect these phenomena to the LGMD attributes lain out here. For instance, it could investigate the shape of localized structures embedded within the auroral electrojet and how they correspond to LGMDs. Similarly, heeding the results of Vandegriff et al. (2024) wherein points on Earth’s surface were magnetically mapped to spatially dispersed locations in the magnetotail, it may be of interest to examine the influence of large-scale features in the magnetosphere on LGMD aspect ratio.

Better establishing the empirical relationship between the solar wind, magnetospheric and ionospheric processes, and the characteristics of LGMDs will provide knowledge about the Sun-to-Earth phenomenological chain, and help to describe the link between LGMDs’ attributes and their drivers.

5 Conclusion

Understanding the complete spatial picture of ground-level geomagnetic disturbances is key when it comes to protecting GIC-susceptible infrastructure and developing our knowledge of magnetospheric physics. We use the Spherical Elementary Current Systems method to interpolate readings from 25 magnetometer stations in a way consistent with physics that drives activity in the magnetosphere. We densely sample the interpolations to create heatmaps showing the behavior of dBN within a wide region of mid-latitude North America.

From hundreds of thousands of these heatmaps, representing nearly all storm time in Solar Cycle 24, we analyze three statistics: the number of LGMDs per heatmap, the perimeters of LGMDs, and the aspect ratios of LGMDs. We develop frequency distributions of these statistics, finding that the typical scale size of LGMDs reconstructed via our methods agrees with scale sizes presented in related studies.

Finally, we subsample these frequency distributions to obtain ones that just represent data collected near solar maximum and ones that represent the same for solar minimum. For each statistic, the distributions from different phases of the solar cycle follow much the same trends; however, there are small differences between LGMDs that occur near solar maximum and ones that occur near solar minimum. These differences are determined to be insignificant by calculating the respective distributions’ probability density functions (PDFs), then comparing their Kullback-Leibler divergences and their total variation distances. In other words, the phase of the solar cycle does not strongly direct the attributes of LGMDs during storm time.

This study opens the path toward continued investigation of LGMDs and their drivers. Some possible areas of future study include investigating any correlation between LGMD attributes and the intensity of the storm during which those LGMDs occur. On the other hand, since we found no influence from the solar cycle, perhaps magnetospheric and ionospheric drivers are considerably more important than solar wind effects in determining the characteristics of LGMDs, and perhaps substorm activity in particular is more relevant than a focus on geomagnetic storms. This would agree with the findings of Ngwira et al. (2018), Ngwira et al. (2025) (among others), which emphasize the behavior of ionospheric currents and the mesoscale magnetospheric structures to which those currents systems are magnetically connected. It may, therefore, be constructive to correlate LGMD attributes with measurements of the ionosphere and aurorae above them, or, taking inspiration from Vandegriff et al. (2024), to correlate LGMD attributes with structures in the magnetotail in physics-based simulations.

Though we have investigated only the northward perturbation component dBN in the present work, it should be noted that our methods could be applied to the eastward (dBE) or horizontal (dBH) components of perturbation. We acknowledge, however, that these quantities do not necessarily behave the same, and the dominant east-west flow of the auroral electrojets may have implications for small current structures driving dBN LGMDs which might not affect current structures that drive dBE LGMDs. Comparing the components’ behaviors as they pertain to this work would make a good topic for future study.

Similarly, our methods have applicability to the time derivative of magnetic field perturbations dB/dt, though in this study we have investigated only the deflection from the magnetic field’s baseline value (dBN, in this case). For comparison to the statistics on dBN LGMDs presented in Section 4.2, we show corresponding distributions for dBN/dt LGMDs in the Supplementary Material. There, LGMD identification is based on intense rate of change instead of on intense perturbation magnitudes. A more thorough analysis of dBN/dt LGMDs has not been performed because of the noise which is introduced by taking a discrete derivative over dBN to obtain dBN/dt.

There is no perfect interpolation method since all such methods are mere models of the physics governing the geomagnetic field. However, one cannot ignore the utility of interpolation in enabling advanced analysis techniques that could not be applied to single-point measurements. When its limitations are understood and accounted for, spatial interpolation has the power to facilitate geospace dynamics research and to enable the development of more effective and actionable predictive models.

Data availability statement

The magnetometer data which was used for this work is from the SuperMAG collaboration, and can be found at https://supermag.jhuapl.edu/. Software created for this work can be found in Mukundan (2025).

Author contributions

RM: Formal Analysis, Writing – review and editing, Methodology, Writing – original draft, Data curation, Investigation, Conceptualization. AK: Supervision, Investigation, Writing – review and editing, Funding acquisition. JM: Investigation, Writing – review and editing, Formal Analysis. VP: Investigation, Formal Analysis, Writing – review and editing. MC: Formal Analysis, Investigation, Writing – review and editing. DH: Writing – review and editing, Funding acquisition.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by NSF EPSCoR Award OIA-1920965 and NSF Award 2331527. VP thanks ANID grant SA772100112 and DICYT Regular grant 042431PA.

Acknowledgments

We thank all members of the MAGICIAN team at UNH and UAF who participated in the discussions leading to this article. We also thank the SuperMAG collaboration for providing the data used in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2025.1610276/full#supplementary-material

References

Akiba, T., Sano, S., Yanase, T., Ohta, T., and Koyama, M. (2019). “Optuna: a next-generation hyperparameter optimization framework,”. New York, NY, USA: Association for Computing Machinery, KDD ’, 2623–2631. doi:10.1145/3292500.3330701

CrossRef Full Text | Google Scholar

Allen, J., Sauer, H., Frank, L., and Reiff, P. (1989). Effects of the march 1989 solar activity. Eos, Trans. Am. Geophys. Union 70, 1479–1488. doi:10.1029/89EO00409

CrossRef Full Text | Google Scholar

Amm, O., and Viljanen, A. (1999). Ionospheric disturbance magnetic field continuation from the ground to the ionosphere using spherical elementary current systems. Earth, Planets Space 51, 431–440. doi:10.1186/BF03352247

CrossRef Full Text | Google Scholar

Anderson, C. W., Lanzerotti, L. J., and MacLennan, C. G. (1974). Outage of the l4 system and the geomagnetic disturbances of 4 August 1972. Bell Syst. Tech. J. 53, 1817–1837. doi:10.1002/j.1538-7305.1974.tb02817.x

CrossRef Full Text | Google Scholar

Aryan, H., Yearby, K., Balikhin, M., Agapitov, O., Krasnoselskikh, V., and Boynton, R. (2014). Statistical study of chorus wave distributions in the inner magnetosphere using Ae and solar wind parameters. J. Geophys. Res. Space Phys. 119, 6131–6144. doi:10.1002/2014JA019939

CrossRef Full Text | Google Scholar

Balasis, G., Balikhin, M. A., Chapman, S. C., Consolini, G., Daglis, I. A., Donner, R. V., et al. (2023). Complex systems methods characterizing nonlinear processes in the near-earth electromagnetic environment: recent advances and open challenges. Space Sci. Rev. 219, 38. doi:10.1007/s11214-023-00979-7

CrossRef Full Text | Google Scholar

Bolduc, L. (2002). GIC observations and studies in the hydro-québec power system. J. Atmos. Solar-Terrestrial Phys. 64, 1793–1802. doi:10.1016/S1364-6826(02)00128-1

CrossRef Full Text | Google Scholar

Borovsky, J. E., and Denton, M. H. (2006). Differences between CME-Driven storms and CIR-Driven storms. J. Geophys. Res. Space Phys. 111. doi:10.1029/2005JA011447

CrossRef Full Text | Google Scholar

Boteler, D. H., and Trichtchenko, L. (2015). “Telluric influence on pipelines,” in Oil and gas pipelines (John Wiley and Sons, Ltd), 275–288. doi:10.1002/9781119019213.ch21

CrossRef Full Text | Google Scholar

Boteler, D. H., and van Beek, G. J. (1999)). August 4, 1972 revisited: a new look at the geomagnetic disturbance that caused the L4 cable system outage. Geophys. Res. Lett. 26, 577–580. doi:10.1029/1999GL900035

CrossRef Full Text | Google Scholar

Campbell, K. M. (2017). A thesis submitted to the graduate faculty of the North Dakota state university of agriculture and applied science. Master’s thesis. Fargo, ND: North Dakota State University.

Google Scholar

Cordell, D., Mann, I. R., Parry, H., Unsworth, M. J., Cui, R., Clark, C., et al. (2024). Modeling geomagnetically induced currents in the Alberta power network: Comparison and validation using hall probe measurements during a magnetic storm. Space weather. 22. doi:10.1029/2023SW003813

CrossRef Full Text | Google Scholar

Dimitrakoudis, S., Milling, D. K., Kale, A., and Mann, I. R. (2022). Sensitivity of ground magnetometer array elements for GIC applications I: resolving spatial scales with the BEAR and CARISMA arrays. Space weather. 20, e2021SW002919. doi:10.1029/2021SW002919

CrossRef Full Text | Google Scholar

Dimmock, A. P., Rosenqvist, L., Welling, D. T., Viljanen, A., Honkonen, I., Boynton, R. J., et al. (2020). On the regional variability of dB/dt and its significance to GIC. Space weather. 18, e2020SW002497. doi:10.1029/2020SW002497

CrossRef Full Text | Google Scholar

Fukushima, N. (1976). Generalized theorem for no ground magnetic effect of vertical currents connected with pedersen currents in the uniform-conductivity ionosphere. Rep. Ionos. Space. Res. Jpn. 30, 35–40.

Google Scholar

Gjerloev, J. W. (2012). The SuperMAG data processing technique. J. Geophys. Res. Space Phys. 117. doi:10.1029/2012JA017683

CrossRef Full Text | Google Scholar

Heyns, M. J., Lotz, S. I., and Gaunt, C. T. (2021). Geomagnetic pulsations driving geomagnetically induced currents. Space weather. 19, e2020SW002557. doi:10.1029/2020SW002557

CrossRef Full Text | Google Scholar

Kappenman, J. G. (2005). An overview of the impulsive geomagnetic field disturbances and power grid impacts associated with the violent sun-earth connection events of 29–31 October 2003 and a comparative evaluation with other contemporary storms. Space weather. 3. doi:10.1029/2004SW000128

CrossRef Full Text | Google Scholar

Kelbert, A. (2020). The role of global/regional Earth conductivity models in natural geomagnetic hazard mitigation. Surv. Geophys. 41, 115–166. doi:10.1007/s10712-019-09579-z

CrossRef Full Text | Google Scholar

Kelbert, A., Bedrosian, P. A., and Murphy, B. S. (2019). “The first 3D conductivity model of the contiguous United States,” in Geomagnetically induced currents from the sun to the power grid (American Geophysical Union AGU), 127–151. doi:10.1002/9781119434412.ch8

CrossRef Full Text | Google Scholar

Kelbert, A., and Lucas, G. M. (2020). Modified GIC estimation using 3-D Earth conductivity. Space weather. 18, e2020SW002467. doi:10.1029/2020SW002467

CrossRef Full Text | Google Scholar

Kullback, S., and Leibler, R. A. (1951). On information and sufficiency. Ann. Math. Statistics 22, 79–86. doi:10.1214/aoms/1177729694

CrossRef Full Text | Google Scholar

Love, J. J., Lucas, G. M., Rigler, E. J., Murphy, B. S., Kelbert, A., and Bedrosian, P. A. (2022). Mapping a magnetic superstorm: march 1989 geoelectric hazards and impacts on United States power systems. Space weather. 20, e2021SW003030. doi:10.1029/2021SW003030

CrossRef Full Text | Google Scholar

Madelaire, M., Laundal, K., Gjerloev, J., Hatch, S., Reistad, J., Vanhamäki, H., et al. (2023). Spatial resolution in inverse problems: the EZIE satellite mission. J. Geophys. Res. Space Phys. 128, e2023JA031394. doi:10.1029/2023JA031394

CrossRef Full Text | Google Scholar

McGranaghan, R. M. (2024). Complexity heliophysics: a lived and living history of systems and complexity science in heliophysics. Space Sci. Rev. 220, 52. doi:10.1007/s11214-024-01081-2

CrossRef Full Text | Google Scholar

McLay, S. A., and Beggan, C. D. (2010). Interpolation of externally-caused magnetic fields over large sparse arrays using spherical elementary current systems. Ann. Geophys. 28, 1795–1805. doi:10.5194/angeo-28-1795-2010

CrossRef Full Text | Google Scholar

Molinski, T. S. (2002). Why utilities respect geomagnetically induced currents. J. Atmos. Solar-Terrestrial Phys. 64, 1765–1778. doi:10.1016/S1364-6826(02)00126-8

CrossRef Full Text | Google Scholar

Mukundan, R. (2025). SEC-sizes: code for characterization of LGMD spatial scale. Release Publ. Front. Astronomy Space Sci. (1.2.1). doi:10.5281/zenodo.13227147

CrossRef Full Text | Google Scholar

National Science and Technology Council (2019). National space weather strategy and action plan

Google Scholar

National Science and Technology Council (2023). Implementation plan of the national space weather strategy and action plan

Google Scholar

Ngwira, C. M., Nishimura, Y., Weygand, J. M., Engebretson, M. J., Pulkkinnen, A., and Schuck, P. W. (2025). Observations of localized horizontal geomagnetic field variations associated with a magnetospheric fast flow burst during a magnetotail reconnection event detected by the THEMIS spacecraft. J. Geophys. Res. Space Phys. 130, e2024JA032651. doi:10.1029/2024JA032651

CrossRef Full Text | Google Scholar

Ngwira, C. M., Pulkkinen, A. A., Bernabeu, E., Eichner, J., Viljanen, A., and Crowley, G. (2015). Characteristics of extreme geoelectric fields and their possible causes: localized peak enhancements. Geophys. Res. Lett. 42, 6916–6921. doi:10.1002/2015GL065061

CrossRef Full Text | Google Scholar

Ngwira, C. M., Sibeck, D., Silveira, M. V. D., Georgiou, M., Weygand, J. M., Nishimura, Y., et al. (2018). A study of intense local dB/dt variations during two geomagnetic storms. Space weather. 16, 676–693. doi:10.1029/2018SW001911

CrossRef Full Text | Google Scholar

Patterson, C. J., Wild, J. A., and Boteler, D. H. (2023). Modeling the impact of geomagnetically induced currents on electrified railway signaling systems in the United Kingdom. Space weather. 21. doi:10.1029/2022SW003385

CrossRef Full Text | Google Scholar

Piccinelli, R., and Krausmann, E. (2014). Space weather and power grids – a vulnerability assessment. Joint Research Centre: Institute for Environment and Sustainability. doi:10.2788/20848

CrossRef Full Text | Google Scholar

Pinto, V. A., Keesee, A. M., Coughlan, M., Mukundan, R., Johnson, J. W., Ngwira, C. M., et al. (2022). Revisiting the ground magnetic field perturbations challenge: a machine learning perspective. Front. Astronomy Space Sci. 9. doi:10.3389/fspas.2022.869740

CrossRef Full Text | Google Scholar

Pirjola, R. (2000). Geomagnetically induced currents during magnetic storms. IEEE Trans. Plasma Sci. 28, 1867–1873. doi:10.1109/27.902215

CrossRef Full Text | Google Scholar

Pulkkinen, A., Amm, O., and Viljanen, A.BEAR working group (2003). Separation of the geomagnetic variation field on the ground into external and internal parts using the spherical elementary current system method. Earth, Planets Space 55, 117–129. doi:10.1186/BF03351739

CrossRef Full Text | Google Scholar

Pulkkinen, A., Bernabeu, E., Eichner, J., Viljanen, A., and Ngwira, C. (2015). Regional-scale high-latitude extreme geoelectric fields pertaining to geomagnetically induced currents. Earth, Planets Space 67, 93. doi:10.1186/s40623-015-0255-6

CrossRef Full Text | Google Scholar

Reyes, P. I., Pinto, V. A., and Moya, P. S. (2021). Geomagnetic storm occurrence and their relation with solar cycle phases. Space weather. 19, e2021SW002766. doi:10.1029/2021SW002766

CrossRef Full Text | Google Scholar

Richardson, I. G., and Cane, H. V. (2012). Solar wind drivers of geomagnetic storms during more than four solar cycles. J. Space Weather Space Clim. 2, A01. doi:10.1051/swsc/2012001

CrossRef Full Text | Google Scholar

Rigler, E. J., Fiori, R. A. D., Pulkkinen, A. A., Wiltberger, M., and Balch, C. (2019). “Interpolating geomagnetic observations,” in Geomagnetically induced currents from the sun to the power grid (American Geophysical Union AGU), 15–41. doi:10.1002/9781119434412.ch2

CrossRef Full Text | Google Scholar

Russell, C. T., Chi, P. J., Dearborn, D. J., Ge, Y. S., Kuo-Tiong, B., Means, J. D., et al. (2008). THEMIS ground-based magnetometers. Space Sci. Rev. 141, 389–412. doi:10.1007/s11214-008-9337-0

CrossRef Full Text | Google Scholar

Schultz, A. (2009). EMScope: a Continental scale magnetotelluric observatory and data discovery resource. Data Sci. J. 8, IGY6–IGY20. doi:10.2481/dsj.SS_IGY-009

CrossRef Full Text | Google Scholar

Snelling, J. M., Johnson, J. R., Willard, J., Nurhan, Y., Homan, J., and Wing, S. (2020). Information theoretical approach to understanding flare waiting times. Astrophysical J. 899, 148. doi:10.3847/1538-4357/aba7b9

CrossRef Full Text | Google Scholar

Upendran, V., Tigas, P., Ferdousi, B., Bloch, T., Cheung, M. C. M., Ganju, S., et al. (2022). Global geomagnetic perturbation forecasting using deep learning. Space weather. 20, e2022SW003045. doi:10.1029/2022SW003045

CrossRef Full Text | Google Scholar

Vandegriff, E. M., Welling, D. T., Mukhopadhyay, A., Dimmock, A. P., Morley, S. K., and Lopez, R. E. (2024). Exploring localized geomagnetic disturbances in global MHD: physics and numerics. Space weather. 22, e2023SW003799. doi:10.1029/2023SW003799

CrossRef Full Text | Google Scholar

Van Der Walt, S., Schönberger, J. L., Nunez-Iglesias, J., Boulogne, F., Warner, J. D., Yager, N., et al. (2014). scikit-image: image processing in python. PeerJ 2, e453. doi:10.7717/peerj.453

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanhamäki, H., and Juusola, L. (2020). “Introduction to spherical elementary current systems,” in Ionospheric multi-spacecraft analysis tools: approaches for deriving ionospheric parameters. Editors M. W. Dunlop, and H. Lühr (Cham: Springer International Publishing), 5–33. doi:10.1007/978-3-030-26732-2_2

CrossRef Full Text | Google Scholar

Viljanen, A. (1998). Relation of geomagnetically induced currents and local geomagnetic variations. IEEE Trans. Power Deliv. 13, 1285–1290. doi:10.1109/61.714497

CrossRef Full Text | Google Scholar

Viljanen, A., Nevanlinna, H., Pajunpää, K., and Pulkkinen, A. (2001). Time derivative of the horizontal geomagnetic field as an activity indicator. Ann. Geophys. 19, 1107–1118. doi:10.5194/angeo-19-1107-2001

CrossRef Full Text | Google Scholar

Viljanen, A., and Tanskanen, E. (2011). Climatology of rapid geomagnetic variations at high latitudes over two solar cycles. Ann. Geophys. 29, 1783–1792. doi:10.5194/angeo-29-1783-2011

CrossRef Full Text | Google Scholar

Webb, D. F. (1991). The solar cycle variation of the rates of CMEs and related activity. Adv. Space Res. 11, 37–40. doi:10.1016/0273-1177(91)90086-Y

CrossRef Full Text | Google Scholar

Wei, L. H., Homeier, N., and Gannon, J. L. (2013). Surface electric fields for North America during historical geomagnetic storms. Space weather. 11, 451–462. doi:10.1002/swe.20073

CrossRef Full Text | Google Scholar

Weygand, J. M., Amm, O., Viljanen, A., Angelopoulos, V., Murr, D., Engebretson, M. J., et al. (2011). Application and validation of the spherical elementary currents systems technique for deriving ionospheric equivalent currents with the north American and Greenland ground magnetometer arrays. J. Geophys. Res. Space Phys. 116. doi:10.1029/2010JA016177

CrossRef Full Text | Google Scholar

Wik, M., Viljanen, A., Pirjola, R., Pulkkinen, A., Wintoft, P., and Lundstedt, H. (2008). Calculation of geomagnetically induced currents in the 400 kV power grid in southern Sweden. Space weather. 6. doi:10.1029/2007SW000343

CrossRef Full Text | Google Scholar

Wing, S., Johnson, J. R., and Vourlidas, A. (2018). Information theoretic approach to discovering causalities in the solar cycle. Astrophysical J. 854, 85. doi:10.3847/1538-4357/aaa8e7

CrossRef Full Text | Google Scholar

Keywords: space weather, geomagnetically induced currents, localized geomagnetic disturbance, ground based magnetometer, spherical elementary current systems

Citation: Mukundan R, Keesee AM, Marchezi JP, Pinto VA, Coughlan M and Hampton D (2025) Localized geomagnetic disturbances: a statistical analysis of spatial scale. Front. Astron. Space Sci. 12:1610276. doi: 10.3389/fspas.2025.1610276

Received: 11 April 2025; Accepted: 21 July 2025;
Published: 13 August 2025.

Edited by:

Georgios Balasis, National Observatory of Athens, Greece

Reviewed by:

Stavros Dimitrakoudis, National and Kapodistrian University of Athens, Greece
Adamantia Zoe Boutsi, National Observatory of Athens, Greece

Copyright © 2025 Mukundan, Keesee, Marchezi, Pinto, Coughlan and Hampton. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Raman Mukundan, cmFtYW4ubXVrdW5kYW5AdW5oLmVkdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.