Mössbauer spectroscopy with ⁵⁷Fe (¹¹⁹Sn) probe layers is a useful method to study the local magnetic structures at buried interfaces. However interface alloying, which always exists in the real samples, have to be taken into account for accurate interpretation of experimental data. We developed an algorithm, which describes the interface intermixing in the multilayers. Substituting deposited atoms by atoms of substrate and floating of deposited atoms in the upper layers during epitaxial growth leads to the formation of asymmetric chemical and magnetic interfaces. This asymmetry in the M₁/M₂ superlattices can explain the difference between magnetic responses from M₁ on M₂ and M₂ on M₁ interfaces which were observed in experiments. Applying this intermixing model to the systems with probe layers located at different distances from the interfaces gives the natural explanation of hyperfine fields distributions on probe atoms and helps us clarify some discrepancies reported in the literature.