Issue 
J. Space Weather Space Clim.
Volume 6, 2016



Article Number  A23  
Number of page(s)  8  
DOI  https://doi.org/10.1051/swsc/2016018  
Published online  13 May 2016 
Research Article
Sfe: waiting for the big one
Observatori de l’Ebre, (OE) CSIC – Universitat Ramon Llull, Roquetes, Spain
^{*} Corresponding author: jjcurto@obsebre.es
Received:
26
October
2015
Accepted:
19
April
2016
Accurate measurements of the radiation delivered during the two largest solar flares ever observed are unavailable. In the case of the Carrington event (1858) the Xray and UV radiation was not recorded, while in the case of the big flare which happened after the storm of 29–31 October 2003 we will call from now on as Halloween event (2003) the radiation saturated the Xray radiometer. Despite many studies, a consensus regarding the real values of these events at the moment of maximum radiation has never been reached.
In this paper, we used an alternative approach to try and determine these values. We estimated the values from the perturbations they produced in the Earth’s magnetism – these are known as Solar Flare Effects (Sfe). Firstly, we established an empirical relationship between the variation in the radiation (cause) and its effect on the magnetism (consequence). Then, using the inverse function, we estimated the energy flux of both events. We found that both flares can actually be classified as being larger than X45.
Finally, we also calculated the return period for a Carringtonlike flare. Assuming that this event had an intensity of about X45 – according to our calculations – we estimated the return period to be 90 ± 60 years.
Key words: Sfe / Solar flares / Carrington event / Return period / Statistics
© J.J. Curto et al., Published by EDP Sciences 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Introduction
Solar flare effects (Sfe) are rapid magnetic variations which are related to the enhancement of the amount of radiation produced during Solar flare events (Curto et al. 1994a). Xray and EUV emissions are the main electromagnetic radiation which cause variations on the electronic density in the ionospheric layers. From the F to the D regions, there are electron density enhancements during solar flares and on Earth the magnetic signature of a flare is visible in the illuminated hemisphere.
Interest in the occurrence and frequency of solar flares has increased in the field of Space Weather because of the perturbations they produce on these variables – the electron density in the ionosphere or the earth’s magnetic field. Both of these are used either actively or passively by key technological systems such as the GPS positioning/guidance system, HF communications, satellite communications, etc. (Lanzerotti 1979, 1983). The impact of severe Space Weather events on domestic and international networks can lead to huge economic costs (Cannon 2013; Schulte in den Bäumen et al. 2014). In this paper, we will concentrate on the radiative aspect of the flares producing Sfe and not on the ejected particles producing geomagnetic storms.
Although Carrington and Halloween events were at the top of the lists of outstanding solar flare effects identified in the scientific literature (Cliver & Svalgaard 2004), there is no consensus regarding the real value of the two largest flares ever observed at their maximum points. The Carrington event (1858) was not registered, and the Halloween event (2003) saturated the Xray detectors (Thomson et al. 2005). In fact, over the last few days of October 2003 a sequence of events happened, but the biggest one as regards radiation was that of 4 November, and, in this paper, it will be known as the Halloween event. This lack of consensus has led to different speculations regarding the true size of these events and we have decided to carry out new estimations from a different viewpoint.
The first observed solar flare, on 1 September 1859 (Carrington 1860; Hodgson 1860), is actually associated with the largest space weather event ever recorded (Cliver 2006). As reported by Cliver & Svalgaard (2004), the Carrington event, as it is commonly called, is at the top of sizeordered lists of solar flare effects (Sfe), or magnetic crochet amplitudes as these observations used to be called.
Increasing interest in extreme Space Weather events for both practical and theoretical reasons has led to the reexamination of various aspects of the Carrington event, especially the flare size (Kiplinger & Garcia 2004; Tranquille et al. 2009; Clarke et al. 2010; Cliver & Dietrich 2013). As most of the mentioned studies agree, the Carrington flare may have been a large one, about the size of a X45–X50 solar flare.
The second flare event we study, that of 4 November 2003, is generally considered to be the most intense SXR event recorded during the satellite period, and occurred during the Halloween event sequence of flares (Gopalswamy et al. 2005). The GOES 18 Å emission in this event saturated at an SXR classification of X17.4. NOAA extrapolated the curve of the register and obtained a maximum of X28 (NOAA 2003). However, according to some authors it may have been a Carringtonlike flare. Thus, Kiplinger & Garcia (2004) reconstructed the light curve to estimate a peak classification of X30.6. Thomson et al. (2004, 2005) and Brodrick et al. (2005) estimated an SXR classification of X45 ± 5 and X34–X48, respectively, but years later Tranquille et al. (2009) deduced a flare classification of X24.8 ± 12.6.
In this study we will assess the magnitude of these two flares by taking into account the magnetic signatures of the Sfe they created.
First, we shall look for an empirical relationship that can express the Sfe amplitude as a function of the flare energy. Then, using the inverse function, we can estimate the energy flux of the flares from the magnetic measurements of the Sfe, when the flare energy could not be measured.
Finally, we will calculate the return period of a Carringtonlike flare assuming that this event had an intensity of about X45 according to our calculations.
1. Data set
First, we identified a series of relevant events with simultaneously available radiation data from satellites and magnetic variations from magnetic observatories on Earth. Original data of Xray flares from GOES satellites belonging to the National Oceanic and Atmospheric Administration (NOAA) have been used (ftp://ftp.ngdc.noaa.gov/STP/spaceweather/solardata/solarfeatures/solarflares/xrays/goes/). The database contains information of flares from 1 September 1975 to the present. In our case, the last date used is 21 May 2015.
As regards UV data, they belong to the CELIAS/SEM experiment on the Solar Heliospheric Observatory (SOHO) spacecraft and they were downloaded from the Space Sciences Center in the University of Southern California which hosts the data (http://www.usc.edu/dept/space_science/semdatafolder/long/15_sec_avg/).
Finally, magnetic data were downloaded from the INTERMAGNET database (http://www.intermagnet.org/datadonnee/downloadeng.php) except for the Carrington case where we digitized the magnetogram of Kew Observatory (UK) for that day using our measuring technique (Curto et al. 1996).
For our data, we chose a sample rate of one minute.
2. Sizing the biggest flares: what
During a solar flare a huge amount of energy is liberated. When this reaches the higher levels of the earth’s atmosphere it produces a large variation in the electron density for a short period of time. These sudden changes in electrical conductivity produce an increase in currents which, in turn, induce magnetic variations on Earth. Their sudden appearance makes them easily identifiable.
2.1. Use of a mathematical model
Our initial objective is to find a simple empirical relationship which relates the Xray and UVray energy flux from a solar flare, with the value of the perturbation produced on the earth’s magnetic field.
There are several issues involved which mean that determining the basic magnitudes of the problem is not easy. One of these is that, on a morphological level, the magnetic variation shape depends heavily on the spectral balance of the group of radiations making up the flare (Curto et al. 2016). The main flux (X and UVray) vary over time and are different from one event to another, thus meaning it is difficult to make a prediction when we do not have direct measurements.
Another problem is that, regarding both radiation but especially magnetic variation, the signal of the event we aim to analyse (during the flare time) represents an increase in the value of a factor which is itself already variable in time and which depends on various, not always quantifiable, factors such as the neutral winds or the electric conductivity in the dynamo region which modulates electric currents producing the magnetic signature of the diurnal variation, Sq (Curto et al. 1994b). It is difficult to separate the variation of the event itself from that which would be seen anyway without this perturbation. Hence, it seems easier to work with the derivatives from the functions as the rate of diurnal variation Sq (T = 24 h) is much lower than that of the Sfe (T approx. minutes) – in this way, we do not need to search for the base levels, a laborious task which involves a certain amount of subjectivity (Curto et al. 1994a). The relations ruling the physics of these phenomena are well known and it is possible to model the Standard behaviour of the Sfe (Curto et al. 1994b). However, when we are dealing with a specific event of which we have few data of the variations in the ionosphere, these models do not apply and we have to resort to a simpler approach. Hence, in this work, for the sake of simplicity and taking into account that the contribution of the Xray and the UVray to ionisation is accumulative and does not interact because they act at different levels of the ionosphere, we have supposed that the ordinary differential equation (ODE) governing the process will be of this type:(1)
Where B(t) is the earth’s magnetic field and X(t − a) and U(t − b) are the functions of the energy flux of Xrays and UVrays, respectively. They are all dependent on time, t. The three (dimensional) parameters α, β and σ will be determined statistically for each flare using multiple linear fits. For the sake of simplicity, we shall only choose magnetic observations above the same meridian as the vortex of the system of currents, guided by the position of the subsolar point. In this way, the magnetic contribution of the Sfe will be mainly made up of the North component of the magnetic field, X. Moreover, we will choose observatories at a medium distance (2000–3000 Km) from the focus so that their magnetic signatures could be expected to be of the same order of magnitude and only dependent on the intensity of the flare. The term σ (the intercept) is a constant to which the variation in the magnetic field caused by other phenomena is attributed, and which is considered to be constant during the flare. A priori, we will choose events which have happened within intervals of time calm enough for any rapid variations with magnetospheric origin (which could distort the calculations) to be considered to be negligible. The parameters a and b are the quantities of time which explain the inertia of the ionosphere to respond to external excitations via ionization radiation which, in our case, synchronize cause and effect.
Next, we apply our model to several Sfe events with different intensities to obtain a representation covering the main range of variation and, thus, to be able to study the behaviour of the parameters in Eq. (1). For this, we chose events where we had all the data (UVray, Xray and magnetism) available and where, specifically, the earth’s magnetic field was measured at a point on Earth. The events were selected in order to have samples in each range of energies producing Sfe. So few samples were needed. After this we calculated the numerical derivative of X(t + a), U(t + a) and B(t), and applied a multiple lineal regression with the earth’s magnetic field as the dependent variable. In this way, we calculated all the parameters α (nT m^{2}/W), β (nT m^{2}/W), σ (nT/min.), a (min.) and b (min.) (Table 1).
Parameters calculated for seven different solar flares. The last column indicates the acronym of the magnetic observatory used in each event.
Inertial parameters have 012 minute values. As in Le et al. (2011), the α and β parameters vary with the intensity of the flare, decreasing as the flare has more power. Regression coefficients, R ^{2}, are quite high which confirms that our mathematical model is sufficiently accurate to reproduce the event. However, according to our experience gained in revising many cases, the fitting works well when the event is pure in the sense that only an excess of radiation is present at that moment. However, if the event was preceded by others, then effects of particle precipitation appear and interfere with the Sfe signal at the time of the flare creating magnetic variations that do not follow our simple relationship.
2.2. Test of validity of the model
As a test of the procedure we have modelled different events where we know all the data. The comparison between real values and modelled ones for two of these events is shown in Figure 1. The dark blue points correspond to the real magnetic field as measured by the observatory, while the light blue points represent the magnetic field calculated from the differential Eq. (1) using the values of Xray and UVray flux measured by satellite. The red points represent the Xray flux contribution to the magnetic field, and the purple ones the UVrays contribution. The calculated values reproduce the measurements very well and we can conclude that the method is valid for making relevant estimations.
Fig. 1. Comparison between real/experimental magnetic field and modelled/theoretical magnetic field calculated using Eq. (1) for two solar flares: (a) 8 November 2013 (X1.1; Observatory: IRT) and (b) 4 June 2007 (M9.0; Observatory: KAK). 
As we can see, the contribution of the Xray component is greater than that of the UVrays in the first case (Fig. 1a), especially in the first part of the perturbation. In the second case (8/11/2013) (Fig. 1b) it is the UVray component which dominates.
2.3. Inverse modelling
As we have stated before, the objective of this work is to be able to estimate the flux of Xrays starting out from the value of magnetic field variations and UVray radiation. To do that, we have to use the same Eq. (1) but, this time, the dependent and independent variables interchange their role:(2)
Using this tool we can carry out the modelling of the case we are especially interested in: the Halloween event. This event has been one of the largest events ever recorded according to most lists of events classified by intensity, such as that of Cliver (2006). However, as previously mentioned, the power of the radiation of this flare has not been accurately measured yet. During the Halloween event, the GOES satellite sensor saturated and it was not possible to measure the maximum value of the flare (NOAA estimated it as X28 in a later reanalysis).
For the Carrington event, the direct use of Eqs. (1) and (2) is not possible because it occurred long before the use of satellites and it was, obviously, impossible to take radiation measurements. However, the earth’s magnetic field was measured during this event and so the sizes of the Sfe causing the event are known. This is what we will use to estimate the intensity of the flares which caused these perturbations. To solve this particular case, we will use UV and Xray relationships found by Le et al. (2011) who, after a statistical analysis of a large set of flare samples, established linear regressions of UV and Xrays (with A and B parameters) which correlate quite well when the flare distance to the central solar meridian (central meridian distance, CMD) is taken into consideration.
Finally, we can find the following relationship linking the Xray and Sfe variations:(3) (4)
According to Le et al. (2011), parameters A and B depend on the CMD and the radiation intensity and, thus, they must be calculated taking into account the position the flare had on the Sun.
3. Results concerning the power of the Halloween and Carrington flares
Although both events have many characteristics in common, their different data availability forced us to use different strategies to obtain their Xray flux. To show the way we proceeded more clearly, a diagram comparing both procedures follows:
Halloween event:
(a1) With data during the time of nosaturation (X_{Halloween} + U_{Halloween} + B_{Halloween}) + equ. 1 > α, β and σ.
(b1) Then, with α, β and σ from (a1) + (U_{Halloween} + B_{Halloween}) + equ. 2 > X_{Halloween} for the whole period.
Carrington event:
(a2) X_{Carrington1} = X_{Halloween} + equ. 5 > α1, β1 and σ1.
(b2) With α1, β1 and σ1 from (a2) + equ. 4 > X_{Carrington2}.
(c2) X_{Carrington2} from (b2) + equ. 5 > α2, β2 and σ2.
(b2) With α2, β2 and σ2 from (c2) + equ. 4 > X_{Carrington3}.
3.1. The Halloween event
At this point it is worth remembering that for this event, the Xray flux sensor produced data during the event except for some minutes of the maxima (between 19:43 and 19:56). The calculation of the Halloween event will be done in two stages: (a) first, with the available information (B, X and U) during the time of nosaturation we will calculate the α, β and σ parameters with the direct relationship 1; (b) later, using the geomagnetic field variation and UV flux data information of the whole event, the parameters calculated in the first stage and the inverse relationship 2 we will find the evolution of the radiation. To calculate the Halloween flare we used magnetic data from the Canadian observatory, Victoria (Fig. 2), which was at the north of the vortex of the current system at that moment.
Fig. 2. Sfe corresponding to 4 November 2003 as recorded at Victoria magnetic observatory in Canada. This event reached its maximum at 19:48 and is popularly known as the Halloween Sfe. 
If we do a multiple linear fit with the data we have (throughout the event except for 14 points), we obtain the results expressed in Table 2.
Parameters for the Halloween event.
With the parameters calculated, we can now model the whole event (Fig. 3).
Fig. 3. Complete temporal evolution of Xray flux for the Halloween event on 4 November 2003 according to our mathematical method. 
The maximum point was at 19:46 and reached an intensity of X46.3 ± 1.8.
On the other hand, taking into account results from Tables 1 and 3, a relationship between power and the values of the respective parameters can be seen (Fig. 4).
Fig. 4. Values of the parameters α (a) and β (b) as a function of the power of the flare. The red line is a fit following a power law. 
Values of the different constants of the statistical fit of parameters α and β.
Both parameters fit well with an exponential function:(5)the values of the constants for each parameter are given in Table 3. Y is the value of the parameter and P the intensity of the flare adjusting the constants Y_{0}, A and t.
Once again we have a high correlation according to the coefficient values of regression R^{2}.
3.2. The Carrington event
To model the Carrington event, we only have the geomagnetic data available, so we cannot calculate the parameters directly with Eq. (1). Hence, we have to use a slightly different strategy.
The data of the magnetic variation (Fig. 5) were obtained by digitalizing the records of the magnetic observatory at Kew (UK) for that day using our measuring technique (Curto et al. 1996).
Fig. 5. Sfe corresponding to 1 September 1858 as recorded at the magnetic observatory in Kew (UK). This event had a maximum peak at 11:17 and is popularly known as the Carrington Sfe. 
The Carrington event is very similar to the Halloween one as regards the magnetic variation amplitude, thus we can suspect that the delivered energy in both cases was similar, so we will use the same value as an initial trial. In fact, the estimation of the Carrington event power was done in an iterative process. First we applied to Eq. (4) the value of the parameters derived from Eq. (5), taking as the initial exploratory value of the Carrington power a value equal to that of the Halloween event. Thus we obtained a second estimated value of the Carrington power. With this second value we initiated another iteration, giving a third value. As this third result coincided with the second one, we stopped the iteration and took this to be the definitive one (Fig. 6).
Fig. 6. Xray flux of Carrington event on 1 September 1858 modelled using our mathematical method, with a maximum of X45.7 ± 2.2 at 11:15. 
The maximum happened at 11:15, reaching an intensity of X45.7 ± 2.2.
Following our mathematical model, we obtain the result that the Halloween event was a X46.3 ± 1.8 flare class and the Carrington event was a X45.7 ± 2.2 flare class. To sum up, at the maximum point, according to our calculations both the Carrington and the Halloween events were higher than X45, doubling the intensities of the next most powerful flares below them on the list of greatest flares ever recorded, as well as being way higher than any of the other recorded flares. This confirms the exceptionality of these events.
4. Return period: when
As humankind are users of a multitude of evermorecomplex technological systems, which are also increasing in vulnerability, we are interested in knowing when we can expect another event so large that it could cause significant perturbations in telecommunication networks, electricity supplies or damage to satellites which offer essential services for a society as technologydependent as ours. We are particularly interested in knowing when a Carringtontype major event will occur again.
Calculation
To calculate the return period of a Carringtonlike flare we used NOAA flare lists which cover the whole period: 1 October 1975 to 21 May 2015, separating the flares into bins of 5 flare classification units starting with C1, i.e. C1–C5, C5.1–C9.9, M1–M5, M5–M9.9, X1–X5, X5–X9.9,… (Fig. 7a). We used the original data except for the case of the Halloween event for which we used our calculations instead of the NOAA data. In general, the number of flares follows a power law with respect to their energy. Like a fractal structure, this happens in each subband we take.
Fig. 7. (a) Number of flares for each bin. Each bin contains 5 standard units. (b) Number of Xclass flares for each bin. Each bin contains 5 standard units. 
Due to a large “weight” of the C and M flares compared to the X flares (Fig. 7a), any attempt to adjust a function will underestimate the Xray population, resulting in a nonrealistic return period calculation for Xrays. Thus, we only consider Xtype flares (Fig. 7b) and, therefore, the functions we use to adjust data are more representative for this range. This allowed us to find reliable results that match with our expectations.
To calculate the return periods for the Xtype flares, we used the Generalized Extreme Value distribution (GEV) (also known as the FisherTippett distribution) because of its versatility. GEV is perfect for our aims because it is very appropriate to model the smallest or largest value among a large set of observations with values randomly distributed, as in our case. GEV combines three simpler distributions – the Gumbel (Gumbel 1954), Fréchet (Fréchet 1927) and Weibull (Weibull 1951) functions – into a single form allowing a continuous range of possible shapes that includes all three of the simpler distributions. Thus, the generalized extreme value distribution allows us to “let the data decide” which distribution is appropriate. Although we tried other distributions, only the GEV distribution adjusted the flare distribution properly.
4.1. GEV fit and results
The adjusted GEV distribution (Ibáñez Rosales 2011) can be written as:(6)where.
And k is the total number of flares considered. In our case, the number of Xray (class X) flares reported in the NOAA database is k = 491. We must multiply the function by this factor because the original GEV function (without the 491 factor) is normalized to 1.
Using this method to fit the GEV distribution to the data, the coefficient values are: σ = 0.198 ± 0.017, μ = 0.799 ± 0.008 and ξ = 0.66 ± 0.10. Figure 8a shows the adjusted curve. The zoom in the Carringtonlike energy range (Fig. 8b) reveals how specific the occurrence of such an event is.
Fig. 8. (a) Number of observed Xray flares with GEV regression. (b) Same representation as in Figure 8a, but with a zoom in the energy range of the Carringtonlike flare. 
The coefficient of determination is R ^{2} = 0.99987. Hence, we can see we had a very precise adjustment.
With the parameters obtained, the return value of a flare with an intensity compressed between a 5unit bin can be calculated as:(7)
Where τ is the period of time since the beginning of the data until the end and x belongs to the bin studied. In our case τ = 40 years. For a flare between X45 and X50, x = 10. This expression gives us a return period of T = 90 ± 60 years.
5. Conclusions
In this work we have calculated the Xray intensity for the Carrington and Halloween events. We then estimated the return period of a Carringtonlike event.
The fact that there has always been controversy regarding the maximum values reached by the two largest events ever measured encouraged us to carry out a new estimation based on corresponding magnetic data on the Earth (Sfe). We first looked for a relationship between the intensity of the radiation incident at earth level (cause) and the amplitude of the magnetic variation produced (consequence). Using this relationship we used inverse modelling to estimate the radiation values at the maximum moment. These were similar for both cases and were above X45. Despite this, it has to be remembered that there was a significant difference between these two flares: Carrington reported seeing, with the naked eye, two intensely white and bright patches of light on the image of the Sun’s disk projected onto a plate of glass. The Carrington flare emitted such a large quantity of radiation in the visible range that it became the first flare ever reported and the first to be linked to a terrestrial response (Cliver 2006).
Afterwards we calculated the Sfe occurrence statistics for the period where we have data from satellites (1973–2015) and we looked for a GEVtype function which, when adjusted to the distribution of the flare happenings in the X band, allowed us to calculate the return periods for extreme events such as the Carrington event.
As mentioned previously, we obtained a return period T_{GEV} = 90 ± 60 years for a Carrington event to happen and a Carringtonlike flare (X45) was seen in 2003, 144 years after the original Carrington flare (1859).
This study allowed us to estimate the return period of a Carringtonlike flare as being 90 years. Having said that, though, we have to bear in mind that our data set is very reduced because radiation data were only available for the last 40 years and the fact that, before this period, only one event, that of the Carrington flare in 1859, had been observed. To obtain more reliable results, a study with more years of data should be carried out. The ideal case would be with hundreds of years of data, – in which case the sample would be more representative (e.g. more great flares (>X25) would have been registered). However, for now, with this limited amount of data, the return period of T = 90 ± 60 years is the best estimate of the real return period that we were able to obtain.
Recent papers looking at several other parameters, including observed CME speeds and the strength of the equatorial current system in the Earth’s magnetosphere, concluded that the risk of an extreme space weather episode could be as high as 12% in a decade (Riley 2012). From this result, the UK Royal Academy of Engineering issued a report on extreme space weather in 2013. They deduced a return period or recurrence interval of 79 years for a superstorm. The usual 95% confidence interval implies we might only have to wait two years for a superstorm, or we might have to wait 300 years. This is a consequence of the nature of randomness. On the other hand, Thomson et al. (2011) quantified the probability of having extreme behaviour in geomagnetic activity as ranging from 100 to 200 years. Thus, different approaches converge in the conclusion that extreme events are very rare and we do not expect one to happen in the near future.
We can conclude that, according to our estimate, humanity will probably have to wait another century to see a big flare located in the Sun’s visible hemisphere with the capacity to produce impacts on Earth via a new Carringtonlike Sfe.
Acknowledgments
This research has been partially supported by Spanish government project CTM201452182C31P of MINECO.
The authors also wish to thank the NOAA for the satellite system which collects valuable information regarding solar radiation and for publishing these data. We acknowledge the CELIAS/ SEM experiment on the Solar Heliospheric Observatory (SOHO) spacecraft (SOHO is a joint European Space Agency, United States National Aeronautics and Space Administration mission). We also thank INTERMAGNET and all the collaborating observatories which provide highquality magnetic data. Finally, we thank the International Service of Rapid Magnetic Variations and their collaborating observatories for providing lists of Sfe events, essential for a study such as ours.
The editor thanks two anonymous referees for their assistance in evaluating this paper.
References
 Brodrick, D., S. Tingay, and M. Wieringa. Xray magnitude of the 4 November 2003 solar flare inferred from the ionospheric attenuation of the galactic radio background. J. Geophys. Res., 110, A09S36, 2005, DOI: 10.1029/2004JA010960. [Google Scholar]
 Cannon, P.S. Extreme Space Weather – a report published by the UK Royal Academy of Engineering. Space Weather, 11, 138–139, 2013, DOI: 10.1002/swe.20032. [CrossRef] [Google Scholar]
 Carrington, R.C. Description of a singular appearance seen in the Sun on September 1, 1859. Mon. Not. Roy. Astron. Soc., 20, 13–14, 1860. [Google Scholar]
 Clarke, E., C. Rodger, M. Clilverd, T. Humphries, O. Baillie, and A. Thompson. An estimation of the Carrington flare magnitude from solar flare effects (sfe) in the geomagnetic records 2010, Royal Astron. Soc. National Astron. Meeting, University of Glasgow, UK, 2010. [Google Scholar]
 Cliver, E.W. The 1859 space weather event: then and now. Adv. Space Res., 38, 119–129, 2006. [CrossRef] [Google Scholar]
 Cliver, E.W., and W.F. Dietrich. The 1859 space weather event revisited: limits of extreme activity. J. Space Weather Space Clim., 3, A31, 2013, DOI: 10.1051/swsc/2013053. [CrossRef] [EDP Sciences] [Google Scholar]
 Cliver, E.W., and L. Svalgaard. The 1859 solarterrestrial disturbance and the current limits of extreme space weather activity. Sol. Phys., 224, 407–422, 2004. [Google Scholar]
 Curto, J.J., C. AmoryMazaudier, J.O. Cardús, J.M. Torta, and M. Menvielle. Solar flare effects at Ebre: regular and reversed solar flare effects, statistical analysis (1953 to 1985), a global case study and a model of elliptical ionospheric currents. J. Geophys. Res., 99 (A3), 3945–3954, 1994a. [Google Scholar]
 Curto, J.J., C. AmoryMazaudier, J.O. Cardús, J.M. Torta, and M. Menvielle. Solar flare effects at Ebre: unidimensional physical, integrated model. J. Geophys. Res., 99 (A12), 23289–23296, 1994b. [Google Scholar]
 Curto, J.J., E. Sanclement, and J.M. Torta. Automatic measurement of magnetic records on photographic paper. Comput. Geosci., 22 (4), 359–368, 1996. [CrossRef] [Google Scholar]
 Curto, J.J., L.F. Alberca, and J. Castell. Dynamic aspects of the Solar flare effects and their impact in the detection procedures. J. Ind. Geophys. Union, 2, 99–104, 2016. [Google Scholar]
 Fréchet, M. Sur la loi de probabilité de l’écartmaximum. Ann. Soc. Polon. Math., 6, 93, 1927. [Google Scholar]
 Gopalswamy, N., L. Barbieri, E.W. Cliver, G. Lu, S.P. Plunkett, and R.M. Skoug. Introduction to violent SunEarth connection events of October–November 2003. J. Gepohys. Res., 110, A09S15, 2005, DOI: 10.1029/2004JA010958. [Google Scholar]
 Gumbel, E.J. Statistical theory of extreme values and some practical applications, Applied Mathematics Series 33, 1st ed., U.S. Department of Commerce, National Bureau of Standards, 1954. [Google Scholar]
 Hodgson, R. On a curious appearance seen in the Sun. Mon. Not. Roy. Astron. Soc., 20, 15, 1860. [Google Scholar]
 Ibáñez Rosales, A. Análisis estadístico de valores extremos y aplicaciones, Report of Universidad de Granada, 2011. [Google Scholar]
 Kiplinger, A.L., and H.A. Garcia. Soft Xray parameters of the great flares of active region 486. Bull. Am. Astron. Soc., 36, 739, 2004. [Google Scholar]
 Lanzerotti, L.J. Geomagnetic influences on manmade systems. J. Atmos. Terr. Phys., 41, 787–796, 1979. [Google Scholar]
 Lanzerotti, L.J. Geomagnetic induction effects in groundbased systems. Space Sci. Rev., 34, 347–356, 1983. [Google Scholar]
 Le, H., L. Liu, H. He, and W. Wan. Statistical analysis of solar EUV and Xray flux enhancements induced by solar flares and its implication to upper atmosphere. J. Geophys. Res., 116 (A11), A11301, 2011, DOI: 10.1029/2011JA016704. [Google Scholar]
 NOAA. Space environment center preliminary report and forecast of solar geophysical data, Rep. 1471, Silver Spring, Maryland, 2003. [Google Scholar]
 Riley, P. On the probability of occurrence of extreme space weather events. Space Weather, 10, S02012, 2012, DOI: 10.1029/2011SW000734. [NASA ADS] [CrossRef] [Google Scholar]
 Schulte in den Bäumen, H., D. Moran, M. Lenzen, I. Cairns, and A. Steenge. How severe space weather can disrupt global supply chains. Nat. Hazards Earth Syst. Sci., 14, 2749–2759, 2014. [Google Scholar]
 Thomson, A.W.P., E.B. Dawson, and S.J. Reay. Quantifying extreme behavior in geomagnetic activity. Space Weather, 9, S10001, 2011, DOI: 10.1029/2011SW000696. [CrossRef] [Google Scholar]
 Thomson, N.R., C.J. Rodger, and R.L. Dowden. Ionosphere gives size of greatest solar flare. Geophys. Res. Lett., 31, L06803, 2004, DOI: 10.1029/2003GL019345. [Google Scholar]
 Thomson, N.R., C.J. Rodger, and M.A. Clilverd. Large solar flares and their ionospheric Dregion enhancements. J. Geophys. Res., 110, A06306, 2005, DOI: 10.1029/2005JA011008. [Google Scholar]
 Tranquille, C., K. Hurley, and H.S. Hudson. The Ulysses catalog of solar hard Xray flares. Sol. Phys., 258, 141–166, 2009. [CrossRef] [Google Scholar]
 Weibull, W. A statistical distribution function of wide applicability. J. Appl. Mech.Trans. ASME, 18 (3), 293–297, 1951. [Google Scholar]
Cite this article as: Curto JJ, Castell J & Moral FD. Sfe: waiting for the big one. J. Space Weather Space Clim., 6, A23, 2016, DOI: 10.1051/swsc/2016018.
All Tables
Parameters calculated for seven different solar flares. The last column indicates the acronym of the magnetic observatory used in each event.
All Figures
Fig. 1. Comparison between real/experimental magnetic field and modelled/theoretical magnetic field calculated using Eq. (1) for two solar flares: (a) 8 November 2013 (X1.1; Observatory: IRT) and (b) 4 June 2007 (M9.0; Observatory: KAK). 

In the text 
Fig. 2. Sfe corresponding to 4 November 2003 as recorded at Victoria magnetic observatory in Canada. This event reached its maximum at 19:48 and is popularly known as the Halloween Sfe. 

In the text 
Fig. 3. Complete temporal evolution of Xray flux for the Halloween event on 4 November 2003 according to our mathematical method. 

In the text 
Fig. 4. Values of the parameters α (a) and β (b) as a function of the power of the flare. The red line is a fit following a power law. 

In the text 
Fig. 5. Sfe corresponding to 1 September 1858 as recorded at the magnetic observatory in Kew (UK). This event had a maximum peak at 11:17 and is popularly known as the Carrington Sfe. 

In the text 
Fig. 6. Xray flux of Carrington event on 1 September 1858 modelled using our mathematical method, with a maximum of X45.7 ± 2.2 at 11:15. 

In the text 
Fig. 7. (a) Number of flares for each bin. Each bin contains 5 standard units. (b) Number of Xclass flares for each bin. Each bin contains 5 standard units. 

In the text 
Fig. 8. (a) Number of observed Xray flares with GEV regression. (b) Same representation as in Figure 8a, but with a zoom in the energy range of the Carringtonlike flare. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.