Ethene metathesis proceeding on Mo-methylidene centers on silica is investigated with density functional theory, applying the cluster approach. Three different locations of the active sites are considered, in which the Mo center replaces a pair of geminal silanols, two silanols from adjacent geminal pairs and two single silanols, respectively. It is shown that metathesis activity of the Mo-methylidene sites strongly depends on their location on silica. Different reactivity of the centers toward alkene is explained by differences in their geometrical and electronic structure parameters. The calculated C-H stretching vibrations of the proposed Mo-methylidene, Mo-ethylidene, and molybdacyclobutane surface complexes are well consistent with the reported IR spectra for the corresponding species generated on real molybdena-silica catalysts. On the basis of the obtained results it is proposed that among the studied cases, the Mo centers replacing two silanols from adjacent geminal pairs of silica surface are the most adequate models of the real active sites.