Defining scale thresholds for geomagnetic storms through statistics

Abstract. Geomagnetic storms, as part of the Sun-Earth relations, are continuously monitored with different indices and scales. These indices usually have some associated scale thresholds to quantify the severity or risk of geomagnetic disturbances. However, the most usual scale thresholds are arbitrarily chosen. In this work we aim to quantify the range of the thresholds through a new method. These new thresholds are based on statistical distribution fitting.

We used different geomagnetic indices, as Dst, SYM-H, and K_p, since they are relevant for space weather purposes. The first two indices have been discriminated between their negative values and the whole dataset. We considered two periods: a short-term one, comprising data from 1997 to 2012; and long-term ones, which are from 1957–2012 for Dst and 1932–2012 for K_p.

We look for the best fit for different statistical continuous distributions applied to these indices. The best fits and the data distribution functions yield to intersects that can be used to define thresholds. The best fit distribution functions are more coincidental between them when considering determined similar datasets, as non-central f-distribution for negative values, meaningful for geomagnetic disturbances; or non-central Student's-t, when the whole dataset is taken. The method yields different values for thresholds depending on the index. Thresholds for geomagnetic storms can be chosen by common values of SYM-H and Dst, as −75 nT for moderate storms; −150 nT for intense storms, and −330 nT for extreme storms. For the case of K_p, the value equal to 5 may mark the departure from quiet time to stormy time.

The obtained values are close to those usually considered as thresholds for, typically, Dst and K_p; therefore the thresholds defined here may provide criteria to assess the vulnerability to geomagnetic activity on design or mitigation purposes.