Estimation of sill matrices in the linear model of coregionalization

Samuel D. Oman, Bella Vakulenko-Lagun

Research output: Contribution to journalArticlepeer-review


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.

Original languageEnglish
Pages (from-to)15-27
Number of pages13
JournalMathematical Geosciences
Issue number1
StatePublished - 2009
Externally publishedYes


  • Algorithm
  • Convergence
  • Direct and cross semivariograms
  • Linear model of coregionalization
  • Quadratic functions

ASJC Scopus subject areas

  • Mathematics (miscellaneous)
  • General Earth and Planetary Sciences


Dive into the research topics of 'Estimation of sill matrices in the linear model of coregionalization'. Together they form a unique fingerprint.

Cite this