Results. The PCA confirms the existence of the Gnevyshev gap (GG) for solar cycles at about 40% from the start of the cycle. The temporal evolution of sunspot area data for even cycles shows that the GG exists at least at the 95% confidence level for all sizes of sunspots. On the other hand, the GG is shorter and statistically insignificant for the odd cycles of aerial sunspot data. Furthermore, the analysis of sunspot area sizes for even and odd cycles of SC12-SC23 shows that the greatest difference is at 4.2-4.6 years, where even cycles have a far smaller total area than odd cycles. The average area of the individual sunspots of even cycles is also smaller in this interval. The statistical analysis of the temporal evolution shows that northern sunspot groups maximise earlier than southern groups for even cycles, but are concurrent for odd cycles. Furthermore, the temporal distributions of odd cycles are slightly more leptokurtic than distributions of even cycles. The skewnesses are 0.37 and 0.49 and the kurtoses 2.79 and 2.94 for even and odd cycles, respectively. The correlation coefficient between skewness and kurtosis for even cycles is 0.69, and for odd cycles, it is 0.90. Conclusions. The separate PCAs for even and odd sunspot cycles show that odd cycles are more inhomogeneous than even cycles, especially in GSN data. Even cycles, however, have two anomalous cycles: SC4 and SC6. According to the analysis of the sunspot area size data, the GG is more distinct in even than odd cycles. We also present another Waldmeier-type rule, that is, we find a correlation between skewness and kurtosis of the sunspot group cycles.

Synoptic maps of solar EUV intensities have been constructed for many decades in order to display the distribution of the different EUV emissions across the solar surface, with each map representing one Carrington rotation (i.e., one rotation of the Sun). This paper presents a new solar EUV synoptic map dataset based on full-disk images from Solar and Heliospheric Observatory/Extreme Ultraviolet Imaging Telescope (SOHO/EIT) and Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). In order to remove the significant and complicated drift of EIT and AIA EUV intensities due to sensor degradation, we construct the synoptic maps in standardized intensity scale. We describe a method of homogenizing the SOHO/EIT maps with SDO/AIA maps by transforming the EIT intensity histograms to AIA level. The new maps cover the years from 1996 to 2018 with 307 SOHO/EIT and 116 SDO/AIA synoptic maps, respectively. These maps provide a systematic and homogenous view of the entire solar surface in four EUV wavelengths, and are well suited, e.g., for studying long-term coronal hole evolution.

We study the asymmetric distribution of weak photospheric magnetic field values in the two hemispheres separately using synoptic maps from SDO/HMI, SOLIS/VSM and WSO during solar cycles 21-24. We calculate the weak-field asymmetry (shift) by fitting the distributions of weak-field values to a shifted Gaussian. Hemispheric shifts derived from the three data sets agree very well, and increase systematically when reducing the spatial resolution of the map. Shifts of the northern and southern hemisphere are typically opposite to each other. Shifts follow the evolution of the trailing flux and have a strong solar cycle variation with maxima in the early to mid-declining phase of the solar cycle. The sign of the hemispheric weak-field shift is always the same as the polarity of the polar field in the respective hemisphere and solar cycle. We also find that shifts in the south are systematically larger in absolute value than in the north.

Context. Systematic observations of magnetic field strength and polarity in sunspots began at Mount Wilson Observatory (MWO), USA in early 1917. Except for a few brief interruptions, this historical dataset continues till present. Aims. The sunspot field strength and polarity observations are critical in our project of reconstructing the solar magnetic field over the last hundred years. Here we provide a detailed description of the newly digitized dataset of drawings of sunspot magnetic field observations. Methods. The digitization of MWO drawings is based on a software package develope d by us. It includes a semi-automatic selection of solar limbs and other features of the drawing, and a manual entry of the time of observations, the measured field strength and other notes hand-written on each drawing. The data are preserved in a MySQL database. Results. We provide a brief history of the project and describe the results from digitizing this historical dataset. We also provide a summary of the final dataset, and describe its known limitations. Finally, we compare the sunspot magnetic field measurements with other instruments, and demonstrate that, if needed, the dataset could be continued using modern observations such as, for example, Vector Stokes Magnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) platform.

We use the synoptic maps of the photospheric magnetic field observed at Wilcox Solar Observatory, Mount Wilson Observatory, Kitt Peak, SOHO/MDI, SOLIS/VSM, and SDO/HMI to study the distribution of weak photospheric magnetic field values in 1974-2018. We fit the histogram distribution of weak field values for each synoptic map of the six data-sets separately with a parametrized Gaussian function in order to calculate the possible shift (to be called here the weak-field asymmetry) of the maximum of the Gaussian distribution from zero. We estimate the statistical significance of the weak-field asymmetry for each rotation. We also calculate several versions of lower-resolution synoptic maps from the high-resolution maps and calculate their rotational weak-field asymmetries. We find that the weak-field asymmetries increase with decreasing map resolution. A very large fraction of weak-field asymmetries are statistically significant, with the fraction of significant weak-field asymmetries increasing with decreasing resolution. Significant weak-field asymmetries of high- and low-resolution maps mainly occur at the same times and have the same sign. Weak-field asymmetries for the different data-sets and resolutions vary quite similarly in time, and their mutual correlations are very high, especially for low-resolution maps. These results give strong evidence for weak-field asymmetries reflecting a real feature of weak field values, which is best seen in medium- and low-resolution synoptic maps and is most likely related to the supergranulation scale of the photospheric field.

The line-of-sight (LOS) component of the large-scale photospheric magnetic field has been observed since the 1950s, but the daily full-disk observations of the full vector magnetic field started only in 2010 using the SOLIS Vector Stokes Magnetograph (VSM) and the SDO helioseismic and magnetic imager (HMI). Traditionally, potential field extrapolations are based on the assumption that the magnetic field in the photosphere is approximately radial. The validity of this assumption has not been tested yet. We investigate here the structure and evolution of the three components of the solar large-scale magnetic field in 2010 - 2017, covering the ascending to mid-declining phase of solar cycle 24, using SOLIS/VSM vector synoptic maps of the photospheric magnetic field. We compare the observed VSM vector magnetic field to the potential vector field derived using the VSM LOS magnetic field observations as an input. The new vector field data allow us to derive the meridional inclination and the azimuth angle of the magnetic field and to investigate their solar cycle evolution and latitudinal profile of these quantities. SOLIS/VSM vector data show that the photospheric magnetic field is in general fairly non-radial. In the meridional plane the field is inclined toward the equator, reflecting the dipolar structure of the solar magnetic field. Rotationally averaged meridional inclination does not have significant solar cycle variation. While the vector radial component $B_r$ and the potential radial component $B_r^{PFSS}$ are fairly similar, the meridional and zonal components do not agree very well. We find that SOLIS/VSM vector observations are noisy at high latitudes and suffer from the vantage point effect more than LOS observations. This is due to different noise properties in the LOS and transverse components of the magnetic field, which needs to be addressed in future studies.

The most important parameter in the coupling between solar wind and geomagnetic activity is the $B_z$-component of the interplanetary magnetic field (IMF). However, recent studies have shown that IMF $B_y$ is an additional, independent driver of geomagnetic activity. We use here local geomagnetic indices from a large network of magnetic stations to study how IMF $B_y$ affects geomagnetic activity at different latitudes for all solar wind and, separately during coronal mass ejections (CMEs). We show that geomagnetic activity, for all solar wind, is 20% stronger for $B_y>0$ than for $B_y<0$ at subauroral latitudes of about $60^{\circ}$ corrected geomagnetic (CGM) latitude. During CMEs, the $B_y$-effect is larger, about 40%, at slightly lower latitudes of about $57^{\circ}$ (CGM) latitude. These results highlight the importance of the IMF $B_y$-component for space weather at different latitudes and must be taken into account in space weather modeling.

Over the past decades and even centuries, the astronomical community has accumulated a signif-icant heritage of recorded observations of a great many astronomical objects. Those records con-tain irreplaceable information about long-term evolutionary and non-evolutionary changes in our Universe, and their preservation and digitization is vital. Unfortunately, most of those data risk becoming degraded and thence totally lost. We hereby call upon the astronomical community and US funding agencies to recognize the gravity of the situation, and to commit to an interna-tional preservation and digitization efforts through comprehensive long-term planning supported by adequate resources, prioritizing where the expected scientific gains, vulnerability of the origi-nals and availability of relevant infrastructure so dictates. The importance and urgency of this issue has been recognized recently by General Assembly XXX of the International Astronomical Union (IAU) in its Resolution B3: "on preservation, digitization and scientific exploration of his-torical astronomical data". We outline the rationale of this promotion, provide examples of new science through successful recovery efforts, and review the potential losses to science if nothing it done.

We study the shape of sunspot cycles using the Wolf sunspot numbers and group sunspot numbers of solar cycles 1-23. We determine the most typical "model" cycles and the most asymmetric cycles, and test the validity of the two Waldmeier rules: the anti-correlation between cycle height and the length of its ascending phase (rule 1), and between cycle height and the length of the preceding cycle (rule 2). We applied the principal component analysis to sunspot cycles and studied the first two components, which describe the average cycle shape and cycle asymmetry, respectively. We also calculated their autocorrelation in order to study their recurrence properties. The best model cycles for Wolf numbers are SC12, SC14, and SC16, the successive even cycles from a long period of rather low overall solar activity. We find that the model cycles in eight different analyses using both sunspot series are almost exclusively even cycles. Correspondingly, the most asymmetric cycles are odd cycles. We find that both Waldmeier rules are valid for the whole Wolf number series of 23 cycles. Waldmeier rule 2 is also valid for group number series although its significance is weaker. Waldmeier rule 1 is not significant for the original group number series, but becomes significant for the proxy series. For separate centuries, Waldmeier rules are not always valid for Wolf numbers and very rarely for group numbers. Our results also offer a new interpretation for the Gnevyshev gap. In addition to being a local depression of solar activity, the Gnevyshev gap is a separatrix that divides cycles into two parts whose relative intensities determine the cycle asymmetry. The Gnevyshev gap is the zero value time of PC2, located approximately 33-42% into the cycle after its start.

We use recently digitized sunspot drawings from Mount Wilson Observatory to investigate the latitudinal dependence of tilt angles of active regions and its change with solar cycle. The drawings cover the period from 1917 to present and contain information about polarity and strength of magnetic field in sunspots. We identify clusters of sunspots of same polarity, and used these clusters to form ``bipole pairs''. The orientation of these bipole pairs was used to measure their tilts. We find that the latitudinal profile of tilts does not monotonically increase with latitude as most previous studies assumed, but instead, it shows a clear maximum at about 25--30 degree latitudes. Functional dependence of tilt ($\gamma$) on latitude ($\varphi$) was found to be $\gamma = (0.20\pm 0.08) \sin (2.80 \varphi) + (-0.00\pm 0.06)$. We also find that latitudinal dependence of tilts varies from one solar cycle to another, but larger tilts do not seem to result in stronger solar cycles. Finally, we find the presence of a systematic offset in tilt of active regions (non-zero tilts at the equator), with odd cycles exhibiting negative offset and even cycles showing the positive offset.

The interaction of the solar wind with the Earth's magnetic field produces geomagnetic activity, which is critically dependent on the orientation of the interplanetary magnetic field (IMF). Most solar wind coupling functions quantify this dependence on the IMF orientation with the so-called IMF clock angle in a way, which is symmetric with respect to the sign of the $B_y$ component. However, recent studies have suggested that the sign of $B_y$ is an additional, independent driver of high-latitude geomagnetic activity, leading to higher (weaker) geomagnetic activity in Northern Hemisphere (NH) winter for $B_y\! >\! 0$ ($B_y\! <\! 0$). In this paper we quantify the size of this explicit $B_y$-effect with respect to the solar wind coupling function, both for Northern and Southern high-latitude geomagnetic activity. We show that high-latitude geomagnetic activity is significantly (by about 40-50\%) suppressed for $B_y\!<\!0$ in NH winter and for $B_y\!>\!0$ in SH winter. When averaged over all months, high-latitude geomagnetic activity in NH is about 12\% weaker for $B_y\!<\!0$ than for $B_y\!>\!0$. The $B_y$-effect affects the westward electrojet strongly but hardly at all the eastward electrojet. We also show that the suppression of the westward electrojet in NH during $B_y\!<\!0$ maximizes when the Earth's dipole axis points towards the night sector, i.e., when the auroral region is maximally in darkness.

Coronal holes (CH) are regions of open magnetic field lines in the solar corona and the source of fast solar wind. Understanding the evolution of coronal holes is critical for solar magnetism as well as for accurate space weather forecasts. We study here the extreme ultraviolet (EUV) synoptic maps at three wavelengths (195A/193A, 171A and 304A) measured by Solar and Heliospheric Observatory/Extreme Ultraviolet Imaging Telescope (SOHO/EIT) and Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) instruments. The two datasets are first homogenized by scaling the SDO/AIA data to the SOHO/EIT level by means of histogram equalization. We then develop a novel automated method to identify CHs from these homogenized maps by determining the intensity threshold of CH regions separately for each synoptic map. This is done by identifying the best location and size of an image segment, which optimally contains portions of coronal holes and the surrounding quiet Sun allowing us to detect the momentary intensity threshold. Our method is thus able to adjust itself to the changing scale size of coronal holes and to temporally varying intensities. To make full use of the information in the three wavelengths we construct, a composite CH distribution, which is more robust than distributions based on one wavelength. Using the composite CH dataset we discuss the temporal evolution of CHs during the solar cycles 23 and 24.

We use the principal component analysis (PCA) to study geomagnetic activity at annual resolution using a network of 26 magnetic stations in 1966-2015, and an extended network of 40 stations in 1980-2015. The first principal component (PC1) describes the long-term evolution of global geomagnetic activity, and has an excellent correlation with indices like the Kp/Ap index. The two networks give identical results for PC1. The second principal component (PC2) is highly correlated with the annual percentage of high-speed streams (HSS). The extended network has a slightly higher sensitivity to HSSs. We verify the non-trivial latitudinal distribution of the second empirical orthogonal function (EOF2). We find that the amplitude of the 22-year variation of geomagnetic activity has a closely similar latitudinal distribution as EOF2. This verifies that the 22-year variation of geomagnetic activity is related to HSSs. The most likely cause is the Russell-McPherron mechanism.

High speed solar wind streams (HSSs) are very efficient drivers of geomagnetic activity at high latitudes. In this paper we use a recently developed $\Delta{H}$ parameter of geomagnetic activity, calculated from the night-side hourly magnetic field measurements of the Sodankyl\"a observatory, as a proxy for solar wind (SW) speed at monthly time resolution in 1914-2014 (solar cycles 15-24). The seasonal variation in the relation between monthly $\Delta{H}$ and solar wind speed is taken into account by calculating separate regressions between $\Delta{H}$ and SW speed for each month. Thereby, we obtain a homogeneous series of proxy values for monthly solar wind speed for the last 100 years. We find that the strongest HSS-active months of each solar cycle occur in the declining phase, in years 1919, 1930, 1941, 1952, 1959, 1973, 1982, 1994 and 2003. Practically all these years are the same or adjacent to the years of annual maximum solar wind speeds. This implies that the most persistent coronal holes, lasting for several solar rotations and leading to the highest annual SW speeds, are also the sources of the highest monthly SW speeds. Accordingly, during the last 100 years, there were no coronal holes of short duration (of about one solar rotation) that would produce faster monthly (or solar rotation) averaged solar wind than the most long-living coronal holes in each solar cycle produce.

Solar coronal holes are sources of high-speed solar wind streams, which cause persistent geomagnetic activity especially at high latitudes. Here we estimate seasonal solar wind speeds at 1 AU for the last 100 years using high-latitude geomagnetic measurements and show that they give information on the long-term evolution of important structures of the solar large-scale magnetic field, such as persistent coronal holes. We find that the centennial evolution of solar wind speed at 1 AU is different for equinoxes and solstices, reflecting differences in the evolution of polar coronal hole extensions and isolated low-latitude coronal holes. Equinoctial solar wind speeds had their centennial maximum in 1952, during the declining phase of solar cycle 18, verifying that polar coronal holes had exceptionally persistent extensions just before the peak of the Grand Modern Maximum of solar activity. On the other hand, solstice speeds had their centennial maximum during the declining phase of solar cycle 23 due to large low-latitude coronal holes. A similar configuration of seasonal speeds as in cycle 23 was not found earlier, not even during the less active cycles of early 20th century. Therefore the exceptional occurrence of persistent, isolated low-latitude coronal holes in cycle 23 is not related to the absolute level of sunspot activity but, most likely, to the demise of the Grand Modern Maximum.

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.

Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.

In the declining phase of the solar cycle, when the new-polarity fields of the solar poles are strengthened by the transport of same-signed magnetic flux from lower latitudes, the polar coronal holes expand and form non-axisymmetric extensions toward the solar equator. These extensions enhance the occurrence of high-speed solar wind streams (HSS) and related co-rotating interaction regions in the low-latitude heliosphere, and cause moderate, recurrent geomagnetic activity in the near-Earth space. Here, using a novel definition of geomagnetic activity at high (polar cap) latitudes and the longest record of magnetic observations at a polar cap station, we calculate the annually averaged solar wind speeds as proxies for the effective annual occurrence of HSS over the whole Grand Modern Maximum (GMM) from 1920s onwards. We find that a period of high annual speeds (frequent occurrence of HSS) occurs in the declining phase of each solar cycle 16-23. For most cycles the HSS activity clearly maximizes during one year, suggesting that typically only one strong activation leading to a coronal hole extension is responsible for the HSS maximum. We find that the most persistent HSS activity occurred in the declining phase of solar cycle 18. This suggests that cycle 19, which marks the sunspot maximum period of the GMM, was preceded by exceptionally strong polar fields during the previous sunspot minimum. This gives interesting support for the validity of solar dynamo theory during this dramatic period of solar magnetism.

In this paper, we study two sets of local geomagnetic indices from 26 stations using the principal component (PC) and the independent component (IC) analysis methods. We demonstrate that the annually averaged indices can be accurately represented as linear combinations of two first components with weights systematically depending on latitude. We show that the annual contributions of coronal mass ejections (CMEs) and high speed streams (HSSs) to geomagnetic activity are highly correlated with the first and second IC. The first and second ICs are also found to be very highly correlated with the strength of the interplanetary magnetic field (IMF) and the solar wind speed, respectively, because solar wind speed is the most important parameter driving geomagnetic activity during HSSs while IMF strength dominates during CMEs. These results help in better understanding the long-term driving of geomagnetic activity and in gaining information about the long-term evolution of solar wind parameters and the different solar wind structures.

Context. The relation between solar surface rotation and sunspot activity still remains open. Sunspot activity has dramatically reduced in solar cycle 24 and several solar activity indices and flux measurements experienced unprecedentedly low levels during the last solar minimum. Aims. We aim to reveal the momentary variation of solar surface rotation, especially during the recent years of reducing solar activity. Methods. We used a dynamic, differentially rotating reference system to determine the best-fit annual values of the differential rotation parameters of active longitudes of solar X-ray flares and sunspots in 1977-2012. Results. The evolution of rotation of solar active longitudes obtained with X-ray flares and with sunspots is very similar. Both hemispheres speed up since the late 1990s, with the southern hemisphere rotating slightly faster than the north. Earlier, in 1980s, rotation in the northern hemisphere was considerably faster, but experienced a major decrease in the early 1990s. On the other hand, little change was found in the southern rotation during these decades. This led to a positive asymmetry in north-south rotation rate in the early part of the time interval studied. Conclusions. The rotation of both hemispheres has been speeding up at roughly the same rate since late 1990s, with the southern hemisphere rotating slightly faster than the north. This period coincides with the start of dramatic weakening of solar activity, as observed in sunspots and several other solar, interplanetary and geomagnetic parameters.

Aims. Sunspot number is a benchmark series in many studies, but may still contain inhomogeneities and inconsistencies. In particular, an essential discrepancy exists between the two main sunspot number series, Wolf (WSN) and group (GSN) sunspot numbers, before 1848. The source of this discrepancy has so far remained unresolved. However, the recently digitized series of solar observations in 1825-1867 by Samuel Heinrich Schwabe, who was the primary observer of the WSN before 1848, makes such an assessment possible. Methods. We construct sunspot series, similar to WSN and GSN, but using only Schwabe's data. These series, called WSN-S and GSN-S, respectively, were compared with the original WSN and GSN series for the period 1835-1867 to look for possible inhomogeneities. Results. We show that: (1) The GSN series is homogeneous and consistent with the Schwabe data throughout the entire studied period; (2) The WSN series decreases by roughly ~20% around 1848 caused by the change of the primary observer from Schwabe to Wolf and an inappropriate individual correction factor used for Schwabe in the WSN; (3) This implies a major inhomogeneity in the WSN, which needs to be corrected by reducing its values by 20% before 1848; (4) The corrected WSN series is in good agreement with the GSN series. This study supports the earlier conclusions that the GSN series is more consistent and homogeneous in the earlier part than the WSN series.

Samuel Heinrich Schwabe made 8486 drawings of the solar disk with sunspots in the period from November 5, 1825 to December 29, 1867. We have measured sunspot sizes and heliographic positions on digitized images of these drawings. A total of about 135,000 measurements of individual sunspots are available in a data base. Positions are accurate to about 5% of the solar radius or to about three degrees in heliographic coordinates in the solar disk center. Sizes were given in 12 classes as estimated visually with circular cursor shapes on the screen. Most of the drawings show a coordinate grid aligned with the celestial coordinate system. A subset of 1168 drawings have no indication of their orientation. We have used a Bayesian estimator to infer the orientations of the drawings as well as the average heliographic spot positions from a chain of drawings of several days, using the rotation profile of the present Sun. The data base also includes all information available from Schwabe on spotless days.

Because of the lack of reliable sunspot observation, the quality of sunspot number series is poor in the late 18th century, leading to the abnormally long solar cycle (1784--1799) before the Dalton minimum. Using the newly recovered solar drawings by the 18--19th century observers Staudacher and Hamilton, we construct the solar butterfly diagram, i.e. the latitudinal distribution of sunspots in the 1790's. The sudden, systematic occurrence of sunspots at high solar latitudes in 1793--1796 unambiguously shows that a new cycle started in 1793, which was lost in traditional Wolf's sunspot series. This finally confirms the existence of the lost cycle that has been proposed earlier, thus resolving an old mystery. This letter brings the attention of the scientific community to the need of revising the sunspot series in the 18th century. The presence of a new short, asymmetric cycle implies changes and constraints to sunspot cycle statistics, solar activity predictions, solar dynamo theories as well as for solar-terrestrial relations.

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 100 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 100 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 1900s to 1960, a dramatic dropout in 1960s and a (mostly weaker) increase thereafter. Moreover, at all stations, the activity at the end of the 20th 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 100 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.

Recent studies of the heliospheric magnetic field (HMF) have detected interesting, systematic hemispherical and longitudinal asymmetries which have a profound significance for the understanding of solar magnetic fields. The in situ HMF measurements since 1960s show that the heliospheric current sheet (HCS) is systematically shifted (coned) southward during solar minimum times, leading to the concept of a bashful ballerina. While temporary shifts can be considerably larger, the average HCS shift (coning) angle is a few degrees, less than the $7.2^{\circ}$ tilt of the solar rotation axis. Recent solar observations during the last two solar cycles verify these results and show that the magnetic areas in the northern solar hemisphere are larger and their intensity weaker than in the south during long intervals in the late declining to minimum phase. The multipole expansion reveals a strong quadrupole term which is oppositely directed to the dipole term. These results imply that the Sun has a symmetric quadrupole S0 dynamo mode that oscillates in phase with the dominant dipole A0 mode. Moreover, the heliospheric magnetic field has a strong tendency to produce solar tilts that are roughly opposite in longitudinal phase. This implies is a systematic longitudinal asymmetry and leads to a "flip-flop" type behaviour in the dominant HMF sector whose period is about 3.2 years. This agrees very well with the similar flip-flop period found recently in sunspots, as well as with the observed ratio of three between the activity cycle period and the flip-flop period of sun-like stars. Accordingly, these results require that the solar dynamo includes three modes, A0, S0 and a non-axisymmetric mode. Obviously, these results have a great impact on solar modelling.

The extension of the sunspot number series backward in time is of considerable interest for dynamo theory, solar, stellar, and climate research. We have used records of the Be-10 concentration in polar ice to reconstruct the average sunspot activity level for the period between the year 850 to the pr esent. Our method uses physical models for processes connecting the Be-10 concentration with the sunspot number. The reconstruction shows reliably that the period of high solar activity during the last 60 years is unique throughout the past 1150 years. This nearly triples the time interval for which such a statement could be made previously

In a recent paper (Usoskin et al., 2002a), we have reconstructed the concentration of the cosmogenic $^{10}$Be isotope in ice cores from the measured sunspot numbers by using physical models for $^{10}$Be production in the Earth's atmosphere, cosmic ray transport in the heliosphere, and evolution of the Sun's open magnetic flux. Here we take the opposite route: starting from the $^{10}$Be concentration measured in ice cores from Antarctica and Greenland, we invert the models in order to reconstruct the 11-year averaged sunspot numbers since 850 AD. The inversion method is validated by comparing the reconstructed sunspot numbers with the directly observed sunspot record since 1610. The reconstructed sunspot record exhibits a prominent period of about 600 years, in agreement with earlier observations based on cosmogenic isotopes. Also, there is evidence for the century scale Gleissberg cycle and a number of shorter quasi-periodicities whose periods seem to fluctuate in the millennium time scale. This invalidates the earlier extrapolation of multi-harmonic representation of sunspot activity over extended time intervals.