Covariance determination as the heart of Least Squares Collocation gravity field modeling is based on fitting an analytical covariance to the empirical covariance, which is stemmed from gravimetric data. The main objective of this study is to process different local covariance strategies over four regions with different topography and spatial data distribution in Iran. For this purpose, Least Squares Collocation based on Remove-Compute-Restore technique is implemented. In the Remove step, gravity reduction in regions with a denser distribution and a rougher topography is more effective. In the Compute step, the assessment of the Collocation estimates on the gravity anomaly control points illustrates that data density is more relevant than topography roughness to have a good covariance determination. Moreover, among the different attempts of localizing the covariance estimation, a recursive approach correcting the covariance parameters based on the agreement between Least Squares Collocation estimates and control points shows better performance. Furthermore, we could see that covariance localization in a region with sparse or bad distributed observations is a challenging task and may not necessarily improve the Collocation gravity modeling. Indeed, the geometrical fitness of the empirical and analytical covariances-which is usually a qualitative test to verify the precision of the covariance determination-is not always an adequate criterion.