Spatial interpolation using nonlinear mathematical programming models for estimation of missing precipitation records

New mathematical programming models are proposed, developed and evaluated in this study for estimating missing precipitation data. These models use nonlinear and mixed integer nonlinear mathematical programming (MINLP) formulations with binary variables. They overcome the limitations associated with spatial interpolation methods relevant to the arbitrary selection of weighting parameters, the number of control points within a neighbourhood, and the size of the neighbourhood itself. The formulations are solved using genetic algorithms. Daily precipitation data obtained from 15 rain gauging stations in a temperate climatic region are used to test and derive conclusions about the efficacy of these methods. The developed methods are compared with some nave approaches, multiple linear regression, nonlinear least-square optimization, kriging, and global and local trend surface and thin-plate spline models. The results suggest that the proposed new mathematical programming formulations are superior to those obtained from all the other spatial interpolation methods tested in this study.