A method for the stochastic simulation of (rain)drop size distributions (DSDs) in space and time using geostatistics is presented. At each pixel, the raindrop size distribution is described by a Gamma distribution with two or three stochastic parameters. The presence or absence of rainfall is modeled using an indicator field. Separable space-time variograms are used to estimate and reproduce the spatial and temporal structures of all these parameters. A simple and user-oriented procedure for the parameterization of the simulator is proposed. The only data required are DSD time series and radar rain-rate (or reflectivity) measurements. The proposed simulation method is illustrated for both frontal and convective precipitation using real data collected in the vicinity of Lausanne, Switzerland. The spatial and temporal structures of the simulated fields are evaluated and validated using DSD measurements from eight independent disdrometers.