The impact of surface effects on thermal induced spin crossover phenomenon is a subject of a broad and current interest. Using the modified Ising-like model of spin crossover solids with the ligand field as function of the molecule' positions and random component on surface by means of Metropolis Monte Carlo algorithm the thermal spin transition curves were calculated. The analysis of spin configuration during transition gives a general idea about contribution of molecules from the surface and inside the lattice into resulting magnetization of the systems. The behavior of hysteresis loop for various surface coupling and fluctuations strength has been described. (C) 2017 Elsevier B.V. All rights reserved.