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.