A reconstruction of sunspot numbers for the last 1000 years was obtained using a sum of sine waves derived from spectral analysis of the time series of sunspot number Rz for the period 1700–1999. The time series was decomposed in frequency levels using the wavelet transform, and an iterative regression model (ARIST) was used to identify the amplitude and phase of the main periodicities. The 1000-year reconstructed sunspot number reproduces well the great maximums and minimums in solar activity, identified in cosmonuclides variation records, and, specifically, the epochs of the Oort, Wolf, Sporer, Maunder, and Dalton Minimums as well the Medieval and Modern Maximums. The average sunspot number activity in each anomalous period was used in linear equations to obtain estimates of the solar radio flux F10.7, solar wind velocity, and the southward component of the interplanetary magnetic field.