The theory of surface magnetism of iron borate, FeBO3, has been extended by taking into consideration a crystal field contribution to the surface magnetic anisotropy energy. For this purpose, a model of distortion of the six-fold oxygen environment of iron ions in the near-surface layer of iron borate has been put forward. The spin Hamiltonian parameters for isolated Fe3+ ions in the distorted environment of the near-surface layer have been calculated using the Newman's superposition model. The crystal field contribution to the surface magnetic anisotropy energy has been calculated in the framework of the perturbation theory. The model developed allows concluding that the distortions of the iron environment produce a significant crystal field contribution to the surface magnetic anisotropy constant. The results of experimental studies of the surface magnetic anisotropy in iron borate can be described assuming the existence of relative contractions in the near-surface layer of the order of 1 %.