Our work proposes a new combined method for interpreting data on planetary magnetic fields, which includes the construction of regional modified S-approximations considering ellipticity and local approximations of a nonstationary field, which is the solution of some parabolic type of equation with constant coefficients. The main provisions of the method of linear integral representations are developed as applied to time-dependent differential operators, which is of fundamental importance in solving many inverse problems of mathematical geophysics.