Abstract

Gumbel’s first distribution is applied to smoothed monthly mean sunspot numbers for solar cycles 10 to 24. According to that, the next minimum for solar cycle 24-25 transition would be the deepest solar minimum of the last 150 years. This study provides an additional insight about changes in the Sun which are currently happening.

1. Introduction

The sunspot cycle is characterized by its variability as much as its regularity. Sunspots are regions of strong photospheric magnetic fields that have the polarities north and south. The spots change with an approximate 11-year period and their magnetic cycle is twice as long. Physical processes contribute to produce the Sun’s 11-year activity cycle as the occurrence of coronal holes. Many models of the solar dynamo exist but none is complete. Part of the solution related to flows in the solar convection zone is studied by helioseismology and large-scale numerical models. By using dynamo theory the cycle's activity is predicted with one parameter of physical analysis [1].

The Sun shows a roughly 11-year cycle of activity, from stormy to quiet cycles and back again (as it is shown in Figure 1).

The solar wind velocity and the strength and direction of the interplanetary magnetic field affect the coupling of energy from the solar wind into the magnetosphere. Some part of this energy is transitorily stored in the magnetosphere and then it is finally deposited in the upper atmosphere in auroral processes and current-driven heating. Properties of the solar wind and interplanetary magnetic field are correlated with the geomagnetic activity (see [24]).

The last solar minimum, currently ending solar cycle 23, is singular. Prior five solar minima occurred in 2008, 1996, 1986, 1976, and 1966. Current cycle had its minimum in December 2008.

The space missions are optimized to study solar maximum conditions and not solar minimum states. The proceedings of the twenty-third series of international workshops related to the SOHO mission show a discussion about if the Sun is at the beginning of a new grand minimum of activity, similar to the historical Dalton and Maunder minima.

The Dalton minimum (1795–1830), the Maunder minimum (1635–1705), and the Spörer minimum (1450–1550) happened when the Sun was quite different from now (e.g., [5]).

Deep solar minima appear randomly and the Sun gives back fast. Last solar minimum was one of the deep solar minima that appear every century [6].

Possibly, the solar minimum is related to the basic understanding of the solar dynamo and its effects on the photosphere and the corona.

According to a long-term analysis [7] it was estimated that the galactic cosmic ray aspect during the recent solar minimum was the largest of previous minima and the increase was more pronounced for locations outside the magnetosphere.

The expanded last solar minimum has produced even more novel conditions in the solar wind and in the Earth’s magnetosphere when compared in opposition to the records obtained since the beginning of the space age. The observed solar wind value of 370 km s−1 is depressed in average value. The predicted values such as the 300 km s−1 in 1912-1913 are not good enough to be compared with the cycle 23-24 values. The median speed of the solar wind is 450 km s−1 and consists of fast streams which rise from 350 km s−1 to 600 km s−1, and when these fast streams vanish, the radiation belts measured by SAMPEX rapidly decline in flux [6]. The last three solar cycles (1976–2008) had enough solar wind measurements to have last solar minimum with a near term perspective; these solar cycles had greater solar activity than normal. This period has been called a grand solar maximum (e.g., [8]).

The magnetic field production was weak in solar cycle 23, and along solar cycles 21 and 22, strong polar fields were produced by this flux transport. The solar magnetic dynamo can change on decadal scales and should reveal how the observed solar field is generated. This behavior was not usual over the last 30–40 years because we did not have uniform records before 1970. The Sun is still behaving strangely in comparison with the last few solar minima [6].

By using the dynamo theory Schatten et al. [1] predicted the cycle’s activity using one parameter of physical analysis.

By considering polar field measurements available for four solar cycles, Svalgaard et al. [9] predicted that the approaching solar cycle 24 will be potentially the smallest cycle in the last 100 years.

Duhau [10] presented evidence that solar activity is in a declining time from the year 1993.

Hathaway [11] studied international sunspot counts and found that the average period of a solar cycle is 131 months with a standard deviation of 14 months. Decaying solar cycle 23 has so far lasted 142 months and is not abnormally depressed. The available 13-smoothed monthly sunspot (MSSM) was 5.70. This is larger than 12 of the last 23 solar minimum values.

To explain changes in space weather, it is necessary to study data and their connections between the effects and the involved physics.

Before last solar minimum, the sunspots and records of solar activity manifested that these years were times of high solar wind; then the Sun will return to less active conditions and will do for a relatively long period of time, according to the fact that it is necessary to monitor carefully the polar regions of the Sun. However, there is much to be learned about the Sun, even at quiet times. There are magnetic regions appearing; there are coronal holes, fast and slow streams, and time variations in the solar wind. There are reasons to think that the low field intensities happened during last solar minimum presage of a low-activity cycle.

To analyze SMMS, Gumbel’s asymptotic distributions [12] are considered. They were applied by Siscoe [13, 14] and Silbergleit et al. [15, 16] to study sunspots in each solar cycle and the largest geomagnetic storms. Analytic expressions are obtained for extreme value populations from which the statistical parameters—mode, median, mean, and standard deviation—are derived.

Therefore, it is important to review the sunspot series in the last centuries as a manner to evaluate changes in the Sun. Figure 1 shows the sunspot numbers related to the last fifteen solar cycles. To analyze the statistics of solar minima, the period from 1848 to 2013 is considered due to the goodness of the data (see [17]), because the sunspots records are considered unsatisfactory during 1700–1748, questionable during 1749–1817, and certain since 1848.

2. Data and Prediction Technique

The data (smoothed monthly mean sunspots, SMMS) were downloaded from ftp://ftp.ngdc.noaa.gov/, that is the web side of the National Geophysical Data Center.

Extreme value statistics are useful to study the tails of series (the lowest or highest values) for data normally distributed. While it is not often important to know the distribution of variables which describe a population, by studying their tails (that appear to be normally distributed) useful information about the distribution itself can be yielded (on the basis of extreme statistics [12]).

For a given minima observation, the probability that this value could be less than is , where Gumbel’s first asymptotic distribution is described by . The probability that this value could be equal to or greater than s is (see [12]).

The theory of extremes gives the mathematical expression of : whereandare constants obtained by a linear square fit.

As the probability function for the maximum amplitude for each solar cycle is unknown, the observed values of are deduced according to [12]. For observed extreme values (where is the total of data), the relationship between and s is obtained.

Extreme magnitudes are calculated by considering the solar cycle data between the years 1848 and 2013 and by using the minimum sunspot numbers per cycle (s) as they are shown in Table 1.

In order to use all fifteen observations, a procedure useful for small samples was introduced by [18]. With the limited objective, a plotting rule [19] for the th of ordered data from a double exponential distribution is obtained as where is the number of observations, is the order of the observation from the minima, and is the estimate of the cumulative frequency of the th term [19]. The related values are defined by The and values are tabulated in columns 3 and 4 of Table 1.

Figure 2 shows the straight line which represents the best fitting adjustment. The slope of the regression line is equal to in units of SMMS−1.

The return periods and are calculated by using the expressions: Equations (4) and (5) define and . These values are plotted in Figure 3. represents the expected times required to find one solar cycle minimum amplitude with the extreme equal to or exceeding abscissa ( branch); shows the expected periods to have one cycle with extreme less than abscissa ( branch). These functions also show the increase in the expected range with a growing number of cycles. and mean that, for a given number of solar cycles (e.g., ), it is expected that values will be confined to a bounded interval. The bounds are estimated by (4) and (5), being one case higher (or lower) than the bound obtained through (4) (or (5)). For both bounds coincide at the median value (mv) which is the midpoint of the distribution (as it is shown in Figure 3).

According to this theory, for sixteen cycles fourteen of them are expected to be in the range , one case lower than the lowest bound and the other one equal to or greater than the highest bound. In reference to Table 1, fourteen solar minima remain within the interval 12.3–0.6; one of them is equal to the highest bound (i.e., 12.3); then, the unknown next minimum should have a value lower than the lowest bound (i.e., 0.6).

As regards [13] the arithmetic mean (am), which is the average of all results, is related to the standard deviation Sd defined by [13] as Sd = 1.2825 a.

The measure of the relative dispersion Rd is obtained to divide Sd by the mode. The statistical characteristic parameters of extreme values obtained are the mode (), the median (mv = 5.2), and the mean (am = 5.6), all in units of SMMS with Sd = 4.0 SMMS and Rd = 1.0.

3. Results, Discussion, and Conclusion

Figure 2 shows the best fitting adjustment obtained considering Gumbel’s first asymptotic distribution. The values of the parameters in (1) areandSMMS.

The statistical characteristic parameters of extreme values obtained are mo = 3.9 sunspots, mv = 5.2, and am = 5.6 SMMS, with Sd = 4.0 SMMS and Rd = 1.0.

Figure 3 shows the return periods versus SMMS. The upper and lower bounds are shown for 16 solar cycles. According to that, fourteen minima are in the interval defined by 0.6 and 12.3 SMMS. One data (in Table 1 data number 15 for solar cycle 22) is on the upper bound; then it is expected that the next sunspot minimum would have a lower value than the lowest bound.

According to our results, the solar cycle 25 minimum could probably be the smallest of the last fifteen cycles. This result gives evidence of sunspots decrease and it agrees with prior articles written by Kane [17], Svalgaard et al. [9], and Duhau [10] who predicted that solar activity is in a declining time and the current solar cycle will be the smallest cycle in the last 100 years. Sunspot numbers and every other indicator of solar activity continued depressed, as they did along 2008 and through mid-2009. According to the present study, the solar cycle 25 minimum would be almost a 60% smaller than the deeper cycle detected before (see Table 1). Thus, we conclude that we should be near the deepest minimum (as measured by smoothed monthly mean sunspots), if the statistics were the same as those for the last fifteen cycles.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Thanks are due to the unknown reviewers for their valuable suggestions and comments. This research was partially supported by PIP N° 11420090100258 (CONICET) and PIP N° 20020090100131 (UBACYT), Argentina.