Considering the spatial patterns of ecosystem services plays an important role in evaluating them. Therefore, in this study, inorder to identify and determin the spatial patterns of ecosystem services, new methods of spatial statistics such as spatial correlation of Global Moran and Anselin Moran Index and GIS were used, after quantifying recreation, aesthetic value and noise reduction ecosystem services in the study area. The results showed that the Global Moran index is up to 0.82, 0.87 and 0.89 respectively, for recreation, aesthetic value and noise reduction ecosystem services. Thus, the examined ecosystem services have a positive correlation and cluster pattern. Considering ecosystem services quantifying results with Local Moran correlation analysis related to distribution of spatial correlation of ecosystem services shown regions with maximum suitability and value of ecosystem services supply have meaningful positive spatial correlation and follow high cluster pattern, while regions with no suitability of ecosystem services supply have negative spatial correlation and low cluster pattern.