Recently Pop ({\it Solar Phys.} {\bf 276}, 351, 2012) identified a Laplace (or double exponential) distribution in the number of days with a given absolute value in the change over a day, in sunspot number, for days on which the sunspot number does change. We show this phenomenological rule has a physical origin attributable to sunspot formation, evolution, and decay, rather than being due to the changes in sunspot number caused by groups rotating onto and off the visible disc. We also demonstrate a simple method to simulate daily sunspot numbers over a solar cycle using the \cite{2012SoPh..276..351P} result, together with a model for the cycle variation in the mean sunspot number. The procedure is applied to three recent solar cycles. We check that the simulated sunspot numbers reproduce the observed distribution of daily changes over those cycles.

A Bayesian method for forecasting solar cycles is presented. The approach combines a Fokker--Planck description of short--timescale (daily) fluctuations in sunspot number (\citeauthor{NobleEtAl2011}, 2011, \apj{} \textbf{732}, 5) with information from other sources, such as precursor and/or dynamo models. The forecasting is illustrated in application to two historical cycles (cycles 19 and 20), and then to the current solar cycle (cycle 24). The new method allows the prediction of quantiles, i.e. the probability that the sunspot number falls outside large or small bounds at a given future time. It also permits Monte Carlo simulations to identify the expected size and timing of the peak daily sunspot number, as well as the smoothed sunspot number for a cycle. These simulations show how the large variance in daily sunspot number determines the actual reliability of any forecast of the smoothed maximum of a cycle. For cycle 24 we forecast a maximum daily sunspot number of $166\pm 24$, to occur in March 2013, and a maximum value of the smoothed sunspot number of $66\pm5$, indicating a very small solar cycle.

Sunspot numbers exhibit large short-timescale (daily-monthly) variation in addition to longer timescale variation due to solar cycles. A formal statistical framework is presented for estimating and forecasting randomness in sunspot numbers on top of deterministic (including chaotic) models for solar cycles. The Fokker-Planck approach formulated assumes a specified long-term or secular variation in sunspot number over an underlying solar cycle via a driver function. The model then describes the observed randomness in sunspot number on top of this driver function. We consider a simple harmonic choice for the driver function, but the approach is general and can easily be extended to include other drivers which account for underlying physical processes and/or empirical features of the sunspot numbers. The framework is consistent during both solar maximum and minimum, and requires no parameter restrictions to ensure non-negative sunspot numbers. Model parameters are estimated using statistically optimal techniques. The model agrees both qualitatively and quantitatively with monthly sunspot data even with the simplistic representation of the periodic solar cycle. This framework should be particularly useful for solar cycle forecasters and is complementary to existing modeling techniques. An analytic approximation for the Fokker-Planck equation is presented, which is analogous to the Euler approximation, which which allows for efficient maximum likelihood estimation of large data sets and/or when using difficult to evaluate driver functions.