dependent of ray sampling, but was not less than 2 km. Processing consisted of several iterations; each iteration contained a source location in the 3D velocity model, matrix calculation and inversion. The inversion was performed simultaneously for the P and S velocities, source corrections (four parameters for each source) and station corrections.Free parameters for the inversion