A latent Gaussian Markov random-field model for spatiotemporal rainfall disaggregation

Rainfall data are often collected at coarser spatial scales than required for input into hydrology and agricultural models. We therefore describe a spatiotemporal model which allows multiple imputation of rainfall at fine spatial resolutions, with a realistic dependence structure in both space and time and with the total rainfall at the coarse scale consistent with that observed. The method involves the transformation of the fine scale rainfall to a thresholded Gaussian process which we model as a Gaussian Markov random field. Gibbs sampling is then used to generate realizations of rainfall efficiently at the fine scale. Results compare favourably with previous, less elegant methods.