A method for computing amplitude independent gradient for seismic velocity inversion in a frequency domain includes receiving seismic data associated with a region. The region comprises one or more earth subsurface layers represented by a plurality of points, and each point is associated with a seismic velocity. Seismic velocities at the plurality of points are determined by iteratively updating the seismic velocities based on a plurality of gradient values, where each gradient value corresponds to a point and is determined by evaluating a gradient of an objective function at a location of the point. A seismic image of the one or more earth subsurface layers is displayed based on the determined seismic velocities.