In the paper, the variation of the intensity of the geomagnetic field force is analysed in time and space. For the research, the data from measurements of the intensity of the geomagnetic field force at four airports (Kaunas, Klaip˙eda, Palanga andVilnius) and 6 geomagnetic field repeat stations aswell as the data from Belsk Magnetometric Observatory (Poland) were used. For the data analysis, the theory of covariance functions was applied. The estimates of the cross-covariance functions of the measured intensity of the geomagnetic field force or the estimates of auto-covariance functions of single data were calculated according to the random functions created from the force intensity measurement data arrays. The estimates of covariance functions were calculated upon varying the quantization interval on the time scale and applying the software created using Matlab package of procedures. The impact of radars of airports on the intensity of geomagnetic field variation and on changes of their covariance functions was established.

The aim of the paper is the comparison of the least squares prediction presented by Heiskanen and Moritz (1967) in the classical handbook “Physical Geodesy” with the geostatistical method of simple kriging as well as in case of Gaussian random fields their equivalence to conditional expectation. The paper contains also short notes on the extension of simple kriging to ordinary kriging by dropping the assumption of known mean value of a random field as well as some necessary information on random fields, covariance function and semivariogram function. The semivariogram is emphasized in the paper, for two reasons. Firstly, the semivariogram describes broader class of phenomena, and for the second order stationary processes it is equivalent to the covariance function. Secondly, the analysis of different kinds of phenomena in terms of covariance is more common. Thus, it is worth introducing another function describing spatial continuity and variability. For the ease of presentation all the considerations were limited to the Euclidean space (thus, for limited areas) although with some extra effort they can be extended to manifolds like sphere, ellipsoid, etc.