Geostatistics and remote sensing techniques are frequently used in analyzing the spatial variability of terrestrial ecosystems. We analyzed spatial variation of moss families by geostatistics and Landsat imagery in a typical semi-arid landscape in North Central Anatolia, Turkey. We sampled 49 sites, chosen based on elevation, slope steepness, and slope aspect. Moss families were determined in a 10-m(2) representative area at each sampling site. The samples were transported to a laboratory and identified for moss families. In total, 10 families were found. Semivariogram analysis was performed to analyze the spatial structure of these families. The semivariogram analysis showed that the moss families were spatially dependent within 117m in the study area. Thirteen thematic classes were categorized by Landsat Thematic Mapper (TM) image in the study area. The classification resulted in an overall kappa statistic of 0.8535, producer accuracy of 74.29, and user accuracy of 86.67. The family with the lowest classification accuracy was Orthotrichaceae (kappa of 0.6379, producer accuracy of 64.52, and user accuracy of 66.67). The moss families and the other classes were identified with a 0.78 kappa statistic value and an 80.74 % accuracy level by using the Landsat TM. The classification showed that Brachytheciaceae, Pottiaceae, Bryaceae, and Grimmiaceae were themost abundantmoss families in this semi-arid environment.