Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-active (storm-like) days and (substorm-like) hours in each year to accuracies of 50-60 pc. The results are the first physics-based quantification of the space weather conditions in both the Dalton and Maunder minima. We predict terrestrial disturbance levels in future repeats of these minima, allowing for the weakening of Earth's dipole moment.

Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data recorded by near-Earth interplanetary spacecraft and by solar magnetograms over the past 53 years. In this paper, we use these reconstructions to quantify power input into the magnetosphere over the past 400 years. For each year, both the annual mean power input is computed and its distribution in daily means. This is possible because the distribution of daily values divided by the annual mean is shown to maintain the same lognormal form with a constant variance. This study is another important step towards the development of a physics-based, long-term climatology of space weather conditions.