Abstract
A new method of spherical harmonic analysis which utilizes both scalar and vector observations is presented and applied to two distinct samples of data obtained between 1900 and 1965. The oblateness of the earth is taken into account and time derivatives up to 2nd order are determined along with the main field. One of the resulting models has been employed in the officially adopted International Geomagnetic Reference Field.
