We investigate the size distribution of the maximum areas and the instantaneous distribution of areas of sunspot groups using the Greenwich sunspot group record spanning the interval 1874-1976. Both distributions are found to be well described by log-normal functions. Using a simple model we can transform the maximum area distribution into the instantaneous area distribution if the sunspot area decay rates are also distributed log-normally. For single-valued decay rates the resulting snapshot distribution is incompatible with the observations. The current analysis therefore supports the results of Howard (1992) and MartinezPillet (1993). It is not possible to distinguish between a linear and a quadratic decay law, however, with the employed data set.

Customary two-dimensional flux transport models for the evolution of the magnetic field at the solar surface do not account for the radial structure and the volume diffusion of the magnetic field. When considering the long-term evolution of magnetic flux, this omission can lead to an unrealistic long-term memory of the system and to the suppression of polar field reversals. In order to avoid such effects, we propose an extension of the flux transport model by a linear decay term derived consistently on the basis of the eigenmodes of the diffusion operator in a spherical shell. A decay rate for each eigenmode of the system is determined and applied to the corresponding surface part of the mode evolved in the flux transport model. The value of the volume diffusivity associated with this decay term can be estimated to be in the range 50--100 km^2/s by considering the reversals of the polar fields in comparison of flux transport simulations with observations. We show that the decay term prohibits a secular drift of the polar field in the case of cycles of varying strength, like those exhibited by the historical sunspot record.