Estimation of sill matrices in the linear model of coregionalization

When using least squares to fit the linear model of coregionalization to multivariate geostatistical data, the sill matrices for the different regions must be estimated, subject to the constraint that they be non-negative definite. In 1992, Goulard and Voltz proposed and empirically examined an iterative algorithm for doing this. Although no proof was given for its convergence or for the uniqueness of the solution to the problem, the algorithm has subsequently been extensively and successfully used. In this paper, we prove that the minimization problem, in fact, has a unique solution and that the algorithm is guaranteed to converge to it from any starting point. We also discuss the effect of the starting point on the speed of convergence.