The large variability of the drop size distribution (DSD) in space and time must be taken into account to improve remote sensing of precipitation. The ability to simulate a large number of 2D fields of DSD sharing the same statistical properties provides a very useful simulation framework that nicely complements experimental approaches based on DSD ground measurements. These simulations can be used to investigate radar beam propagation through rain and to evaluate different radar retrieval techniques. The proposed approach uses geostatistics to provide structural analysis and stochastic simulation of DSD fields. First, the DSD is assumed to follow a Gamma distribution with three parameters. Therefore, 2D fields of DSDs can be described as a three-component multivariate random function. Such 2D fields are normalized using a Gaussian anamorphosis and simulated by taking advantage of fast multivariate Gaussian simulation algorithms. Variograms and cross-variograms are used to generate fields with identical spatial structure that are consistent with the observations. To assess the proposed approach, the method is applied to data collected during intense Mediterranean rainfall. Taylor’s hypothesis is assumed to convert time series into 1D range profiles. The anisotropy of the DSD fields is derived from radar reflectivity measurements collected over the same area during the same event. A large number of DSD fields are generated and the corresponding reflectivity fields are derived. The results of the simulations are in good agreement with respect to the mean, the standard deviation and the spatial structure, demonstrating the promising potential of the proposed stochastic model of DSD fields