Did Open Solar Magnetic Field Increase during the Last 1 Years: A Reanalysis of Geomagnetic Activity arxiv:astro-ph/1167v1 16 Nov K. Mursula Department of Physical Sciences, University of Oulu, Finland; email: Kalevi.Mursula@oulu.fi D. Martini Department of Physical Sciences, University of Oulu, Finland, and GGKI, SOPRON, Hungary A. Karinen Department of Physical Sciences, University of Oulu, Finland Abstract. Long-term geomagnetic activity presented by the aa index has been used to show that the heliospheric magnetic field has more than doubled during the last 1 years. However, serious concern has been raised on the long-term consistency of the aa index and on the centennial rise of the solar magnetic field. Here we reanalyze geomagnetic activity during the last 1 years by calculating the recently suggested IHV (Inter-Hour Variability) index as a measure of local geomagnetic activity for seven stations. We find that local geomagnetic activity at all stations follows the same qualitative long-term pattern: an increase from early 19s to 196, a dramatic dropout in 196s and a (mostly weaker) increase thereafter. Moreover, at all stations, the activity at the end of the th century has a higher average level than at the beginning of the century. This agrees with the result based on the aa index that global geomagnetic activity, and thereby, the open solar magnetic field has indeed increased during the last 1 years. However, quantitatively, the estimated centennial increase varies greatly from one station to another. We find that the relative increase is higher at the high-latitude stations and lower at the low and mid-latitude stations. These differences may indicate that the fraction of solar wind disturbances leading to only moderate geomagnetic activity has increased during the studied time interval. We also show that the IHV index needs to be corrected for the long-term change of the daily curve, and calculate the corrected IHV values. Most dramatically, we find the centennial increase in global geomagnetic activity was considerably smaller, only about one half of that depicted by the aa index. Keywords: geomagnetic activity, solar open field, long-term change 1. Introduction One of the most interesting and important questions in solar-terrestrial physics is whether the magnetic activity of the Sun has indeed greatly increased during the last 1 years. A significant increase in solar magnetic activity is indicated, e.g., by the well known fact that the average amplitude of sunspot cycles during the latter part of the th century is higher than in the beginning. The increasing sunspot activity leads, c 18 Kluwer Academic Publishers. Printed in the Netherlands. IHV_paper_3.9..tex; /11/18; 8:56; p.1
K. Mursula, D. Martini, and A. Karinen Table I. Information on stations used. Magnetic coordinates are calculated using the IGRF model. Station IAGA Geographic Geomagnetic Midnight Data Code Lat Long Lat Long UT hour start Sodankylä SOD 67.7 6.6 63.96 1.5 191 Sitka SIT 57.5.67 6.33 79.79 9 19 Eskdalemuir ESK 55.3 356.8 57.86 83.85 1911 Niemegk NGK 5.7 1.68 51.89 97.69 3 191 Cheltenham CLH 38.73 83.16 9.1 353.71 5 191 Fredericksburg FRD 38. 8.63 8.59 353.11 5 1956 Tucson TUC 3.5 9.17.6 315.63 7 199 Honolulu HON 1.31 1.91 1.57 69.37 1 19 according to a simple model presented by Solanki et al. (), to a corresponding long-term increase in the total solar magnetic field, as well as in the open solar magnetic field, i.e., in the heliospheric magnetic field (also called the interplanetary magnetic field). Heliospheric magnetic field is known to be the main modulator of cosmic rays and a significant cause of geomagnetic activity. Lockwood et al. (1999) used the long-term geomagnetic aa index to show that the heliospheric magnetic field is now more than twice stronger than 1 years ago. Recently, using long-term measurements of cosmogenic isotopes and a sophisticated chain of physical models, it was shown (Usoskin et al., 3) that the solar cycle averaged sunspot activity since 19s is higher than at any other time during the last 1 years. Despite the versatility and uniformity of these important results, serious concern has recently been raised on the centennial rise of solar magnetic activity. E.g., the long-term consistency of the geomagnetic aa index on which perhaps the most reliable evidence on centennial rise is based, has been seriously questioned (Clilverd et al., ). Because of these problems with the aa index, Svalgaard et al. (3) introduced the IHV index as a more straightforward and homogeneous measure of local long-term geomagnetic activity. Using the data from the Cheltenham/Fredricksburg station pair, they found no evidence for an increase in the corresponding IHV index during the last 1 years. In this paper we calculate the IHV index for several geomagnetic stationsinordertoobtainamoreglobal andreliableviewonthelong-term development of global geomagnetic activity and in order to conclude whether or not there is reliable evidence for a centennial increase of IHV_paper_3.9..tex; /11/18; 8:56; p.
Geomagnetic activity during the last 1 years 3 Daily variation [nt] 1 5 5 CLH/FRD 1 193 15 19 1986 5 1 15 5 a 1 1 CLH/FRD 1917 198 199 3 5 1 15 5 b Daily variation [nt] SOD 193 19 1986 6 5 1 15 5 hours (UT) c 1 5 5 SOD 1917 198 199 1 5 1 15 5 hours (UT) Figure 1. Daily curves for CLH/FRD and SOD stations for some sunspot minimum (a and c) and maximum (b and d) years. Time is given in UT hours. d geomagnetic activity and, thereby, of the open solar magnetic field. Moreover, we study in detail the daily variation of the magnetic field, its long-term change and its effect on the IHV index.. Data and daily variation We use here data from six stations and one station pair which have a long record of magnetic observations from early 19s. The codes, coordinates, local midnight UT hours, and start years of these stations are depicted in Table I. (The stations are listed in the order of decreasing latitude). CLH and FRD form a station pair and is included here in order to allow a comparison with the results by Svalgaard et al. (). Figure 1 depicts the yearly averaged daily variation curves of the H-component for CLH/FRD and SOD for some sunspot minimum and maximum years.(the daily average was subtracted in each case in order to remove the effect of the secular variation of the magnetic field). One can see that while the daily curves at the two stations are very different (even if shifted so as to coincide in LT time), they remain fairly similar inform over thenearly 1 year timeinterval. We remindthat the daily curve is, especially at low and mid-latitudes like CLH/FRD, mostly due to the dayside ionospheric S q (quiet day) current. IHV_paper_3.9..tex; /11/18; 8:56; p.3
K. Mursula, D. Martini, and A. Karinen 15 SOD a SIT b 1 3 5 5 1 6 5 CLH/FRD c 15 HON d 3 1 1 5 Figure. Yearly averages of the daily range at stations (left axis in nt units, thick lines) together with the aa index (panels a and b; right axis in nt units, thin lines) and sunspot numbers (panels c and d; right axis, thin lines). Figure 1 shows that although the form of the daily curve remains roughly the same each year, its range (difference between daily maximum and minimum) varies greatly, especially between the minimum and maximum sunspot years. Figure depicts the yearly averaged range of the daily curve for a number of stations over the whole time interval. The variation of the daily range with solar cycle is apparent. Moreover, one can see in Figure that the similarity between the daily range and sunspot number curves (and the correlation between the corresponding values) is excellent for the low and mid-latitude stations (Figure c and d). This is in accordance with the fact that the S q variation is mainly due to solar UV radiation, and thereby, is closely correlated with sunspot numbers. However, at high latitudes (Figure a and b), where current systems other than the S q current start dominating the daily curve, the daily range curves correlate better with geomagnetic activity than sunspot numbers. IHV_paper_3.9..tex; /11/18; 8:56; p.
Geomagnetic activity during the last 1 years 5 6 SOD a 3 SIT k 1 15 ESK c 1 5 aa e 8 TUC g 6 1 8 NGK d 6 8 CLH/FRD f 6 6 HON h Figure 3. Yearly averages of the uncorrected IHV index (in nt) for the 7 stations together with the aa index (panel e). 3. Original IHV index Figure 1 shows that the CLH/FRD daily curve has a smooth section (almost a plateau) in the local pre-midnight hours from to 6 UT prior to an increase at about 1 UT and a large depletion at local noon. Because of this smoothness (and because this time coincides with the geomagnetically active pre-midnight local sector) Svalgaard et al. () used the CLH/FRD data in their analysis and defined the IHV index as an average of the six absolute differences of the successive hourly values of the H component between and 6 UT (19-1 LT). We have calculated the IHV index (from the H component) for all the eight stations of Table I to form seven long-term series of local geomagnetic activity. (CLH and FRD IHV values were calculated separately and joined together in the same way as in Svalgaard et al., ). Figure 3 depicts the seven series of yearly averaged IHV values and the aaindexforcomparison. Theabsolutevaluesof theihv indicesvary, as expected, greatly with the magnetic latitude of the station so that the values at the highest SOD station are roughly an order of magnitude IHV_paper_3.9..tex; /11/18; 8:56; p.5
6 K. Mursula, D. Martini, and A. Karinen Table II. Values of the IHV, IHV-q and IHV-cor indices in the beginning (191-19) and at the end (1979-) of the last century. Station IHV IHV Rel. IHV-q IHV-q Rel. IHV-cor IHV-cor Rel. start end change start end change start end change SOD 8.65. % 1.6 18.81 9% 16.1 1. 3% SIT 8.59 13.89 6%.1 5.37 153% 6.8 8.5 31% ESK 3.96 6.77 71% 1.51 1. -18%.5 5.5 15% NGK 5.11 6.19 1%.69.8 19%. 5.37 1% CLH/FRD.7.73 6%.6.7 %.1.6 6% TUC.63 5. 17% 1.8 1.35 5% 3.55.6 1% HON 3.37.15 3%.87 1.19 37%.5.96 18% larger than at the lowest HON station. Despite this difference, all IHV curves depict a clear solar cycle variation which more closely follows the cyclicity of the aa index than of sunspots. In fact, the IHV index seems to present a fairly good proxy of local geomagnetic activity, as suggested by Svalgaard et al. (). Moreover, all the seven IHV series, as well as the aa index, depict the same qualitative long-term pattern during the last 1 years: There is an increase of activity from the beginning of the th century until 196, then a dramatic dropout in early 196s, and (in most stations) weaker increasing activity therafter. We have underlined this pattern in Figure 3 for each station by including the best fitting line for the period until 196 and another line for 1963-. Because of the (often overall) maximum in 196 there is no uniform increase in geomagnetic activity during the last 1 years and a two-line fit presents this steplike behaviour better than a one-line fit over the full time interval. However, the centennial increase is evident and in all the seven series the IHV index attains a higher average level at the end of the century than in the beginning. We have quantified this increase by calculating the average values of the IHV index during the last (1979-) and first (191-19) years of the century. (Because of different start years, the stations cover a slightly different fraction of the first years). We have depicted these average levels as well as the implied percentual change (increase) of local geomagnetic activity in Table II. This qualitatively uniform development in all the seven series proves that geomagnetic activity has indeed increased during the last 1 years, as has earlier been claimed on the basis of the aa index. Note, however, that the relative increase (column of Table II) depicted by the seven series varies greatly, being higher at the highlatitude stations (up to 7% at ESK) and lower at the low to midlatitude stations. While this difference can not yet be fully understood, it may indicate that the fraction of those disturbances in the solar IHV_paper_3.9..tex; /11/18; 8:56; p.6
Geomagnetic activity during the last 1 years 7 wind that cause only moderate geomagnetic activity (like substorms) observed mainly at high latitudes has increased during this time interval. Also, the CLH/FRD IHV series is exceptional in depicting by far the smallest increase of about 6% only, less than half of the second smallest increase at TUC and a tenth of the increase at ESK and SIT. This exceptionally weak increase of local geomagnetic activity at CLH/FRD made Svalgaard et al. () to erroneously doubt the increase of global geomagnetic activity during the last 1 years. Although we find that there is no doubt that global geomagnetic activity has increased during the last 1 years, the exact amount of this increase is not completely unambiguous. The corresponding relative increase in the aa index is 65%, i.e., close to the percentual increase at SIT and ESK. However, a considerably lower increase is found at the mid-latitude stations (NGK, CLH/FRD; about 5 ) that are closest in magnetic latitude to the aa index stations. We can now make a conservative estimate on the centennial increase in global geomagnetic activity based on the results of Table II. Selecting one high-latitude (SIT), one mid-latitude (NGK) and one low-latitude (HON) station from among those stations whose observations start either in 191 or 19, we obtain an average centennial increase of 35%. (We neglect here CLH/FRD since, as mentioned by Svalgaard et al.,, the values before 1915 may be overestimated). This is also very close to the overall averaged increase of all stations of about 3%. Accordingly, we find that the centennial increase is considerably lower, roughly by a factor of two, than depicted by the aa index.. Corrected IHV index As discussed above (see Figures 1 and ), the range of daily curve varies strongly with solar activity (and geomagnetic activity at high latitudes). Since the sunspot cycles depict a roughly similar long-term evolution (increase of cycle amplitudes until SC 19 in 1957; dropout and increase) as geomagnetic activity, it is clear that the daily range is also affected by this long-term change, as can also be seen in Figure. Because of its very definition,theihv index dependson therangeof the daily curve during the 7 pre-midnight hours (19-1 LT). Therefore, the long-term change of the daily range is expected to contribute to the long-term change of the IHV index. Moreover, since geomagnetic activity is defined as a deviation from the quiet time daily curve, the long-term variation in the daily range has to be removed from the IHV index. We have done this as follows. We first calculated the yearly averaged daily curves for each station in order to obtain a proxy for the IHV_paper_3.9..tex; /11/18; 8:56; p.7
8 K. Mursula, D. Martini, and A. Karinen 6 5 SOD a 5 SIT b 3 15 1 1 5 7 6 CLH/FRD c 6 5 HON d 5 3 3 1 Figure. Yearly averages of the corrected (thick lines) and uncorrected (thin lines) IHV index for stations. Best fitting lines to IHV-cor are also included. quiet time daily variation in each year. Then we calculated the yearly quiet-time IHV-q value from these smooth yearly curves in the same way as the IHV index is calculated for each day from the raw data. The averaged values of the IHV-q index at the start (191-19) and end (1979-) of the 1-year time interval (and the relative change) are included in Table II. Note that the smallest IHV-q values are found for CLH/FRD, in agreement with the flatness of the daily curve during the IHV hours (see Figure 1). Considerably larger IHV-q values are found in many other stations, in particular in SOD where the daily curve is very steep in the local pre-midnight hours (see Figure 1). Note that in most stations the relative long-term change in IHV-q is positive and of the same order of magnitude as the relative change in IHV, except at ESK where it is negative. We have formed the corrected IHV-cor index by subtracting the yearly IHV-q values from the original daily IHV index. Figure depicts the yearly averaged IHV-cor indices for a few stations together with the original IHV index. We have also depicted the best fitting lines to the IHV-cor series in Figure. Note that the correction seems to be IHV_paper_3.9..tex; /11/18; 8:56; p.8
Geomagnetic activity during the last 1 years 9 Table III. The slope (multiplied by 1 ) of the best fitting line for normalized IHV (original and corrected) values and aa index. Mean and standard deviation are obtained from IHV values only. Data SOD SIT ESK NGK aa CLH/FRD TUC HON Mean St.dev. IHV/aa 55 7 61 11 6 3 16 IHV-cor 35 3 77 1 19 3 1 more important at high and low latitudes while at mid-latitudes the correction is smaller. The averaged values of the IHV-cor index at the start and end of the 1-year time interval (and the relative change) are included in Table II. A faster long-term increase in IHV-q than in IHV leads to a smaller increase in the corrected IHV-cor index. This is the situation in SOD, TUC, and HON where the centennial increase in IHV-cor is some 15-% smaller than in IHV. Even a larger reduction is found for SIT where the increase in the corrected IHV-cor is nearly one half of the original IHV value. In the other (interestingly, all midlatitude) stations the reverse is true and there is a larger centennial increase in IHV-cor than in the original IHV. The extreme is, because of the exceptional change of the daily curve, the ESK station where the IHVcor depicts some 6% larger increase than IHV. Thus, the net result of the daily curve correction is that the average centennial increase in the IHV-cor remains nearly the same, at about 35%, as in the original IHV index. However, if the centennial increase is calculated from the same three longest stations as above (SIT, NGK, HON), a somewhat smaller value of about % is obtained for the IHV-cor. Moreover, in order to further study the centennial increase in the various stations we have normalized both the original and the corrected IHV series at each station by their respective averages and depicted the slopes of the best fitting lines in Table III. This Table verifies the above results of the higher centennial increase at high latitudes and the different effect of the daily curve correction at high and low latitudes vs. mid-latitudes. Moreover, Table III also gives the average slope of the centennial increase (and its standard deviation) at the seven stations. Note that this average slope is roughly the same, about.3, bothfor the original and corrected IHV index. This average slope corresponds to the average centennial increase of about 38%, in a very good agreement with the average increase estimated above using the -year time intervals, but remains clearly below the corresponding slope (about.61) for the aa index. These results further underline our main result that the centennial increase is indeed real but considerably smaller than (roughly one half of) that depicted by the aa index. IHV_paper_3.9..tex; /11/18; 8:56; p.9
1 K. Mursula, D. Martini, and A. Karinen 5. Conclusions We have reanalyzed geomagnetic activity during the last 1 years by calculating the recently suggested IHV index as a measure of local geomagnetic activity for 7 stations (or station pairs). We find that local geomagnetic activity in all stations follows the same qualitative longterm pattern as depicted by the geomagnetic aa index: an increase from early 19s to 196, a dramatic dropoutin 196s and a (mostly weaker) increase thereafter. In all stations the average local geomagnetic activity during the last two decennia of the th century is stronger than during the first two decennia. This verifies the result based on the aa index that global geomagnetic activity, and thereby, the open solar magnetic field has indeed increased during the last 1 years. However, we find that the amount of centennial increase varies greatly from one station to another. In particular, a very small increase is found in the CLH/FRD series, which was erroneously interpreted in evidence of no centennial increase. We find that the relative increase is higher at the high-latitude stations and lower at the low to midlatitude stations. These differences may indicate that the fraction of solar wind disturbances leading to only moderate geomagnetic activity has increased during the studied time interval. We have also shown that the IHV index needs to be corrected for the long-term change of the daily curve, and we have calculated the corrected IHV values. In four out of seven stations this correction reduces the estimated centennial increase of local geomagnetic activity. Summarizing, our results give strong evidence that the centennial increase indeed took place during the last century. However, the increase of about 65% depicted in the aa index is roughly twice larger than our more global estimate of about -35% based on seven series of long-term observations. References Clilverd, M. A., T. D. G. Clark, E. Clarke, H. Rishbeth, and T. Ulich:, J. Geophys. Res., 17, SSH -1, CiteID 11, DOI 1.19/1JA51. Lockwood, M., Stamper, R., and Wild, M. N.: 1999, Nature, 399, 37. Solanki, S. K., Schussler, M., and Fligge:, Nature, 8, 5. Svalgaard, L., E. W. Cliver, and P. Le Sager:, Adv. Space Res., 3, 3639. Usoskin, I. G., S. Solanki, M. Schussler, K. Mursula, and K. Alanko: 3, Phys. Rev. Lett., 91, 1111. IHV_paper_3.9..tex; /11/18; 8:56; p.1