Spatial auto-regressive (SAR) models are widely used in geosciences for spatial analyses; their main feature is the presence of weight (W) matrices, which define the neighboring relationships between the spatial units. The statistical properties of parameter and forecast estimates strongly depend on the structure of such matrices. The least squares (LS) method is the most flexible and can estimate systems of large dimensions; however, it is not consistent in the presence of multilateral (sparse) matrices. Instead, the unilateral specification of SAR models provides triangular weight matrices which allow good statistical properties to LS and enable the implementation of sequential prediction functions. In this paper we show the better performance in out-of-sample forecasting of unilateral SAR and LS with respect to multilateral and maximum likelihood (ML) methods. This conclusion is supported by extensive numerical simulations and applications to real geological data, both on regular lattice and irregular polygonal form.