The invention relates to a method of determining a susceptibility gradient map, the method comprising acquiring (300) a magnetic resonance k- space dataset, said dataset comprising magnetic resonance echo data, wherein a recursive iteration is used to determine the susceptibility gradient map. Further, the invention relates to a method of determining a susceptibility gradient map, the method comprising acquiring a magnetic resonance k-space dataset, said dataset comprising magnetic resonance echo data, wherein the acquired dataset comprises magnetic resonance echo data of multiple echoes with different echo times, wherein a susceptibility gradient map is determined separately for each echo time resulting in a temporary susceptibility gradient map, wherein the method further comprises calculating a total susceptibility gradient map by combining all determined temporary susceptibility gradient maps.