Abstract: Based on the data from Chinese National Forest Inventory (CNFI) and Key Ecological Benefit Forest Monitoring plots (5075 in total) in Heilongjiang Province in 2010 and concurrent meteorological data coming from 59 meteorological stations located in Heilongjiang, Jilin and Inner Mongolia, this paper established a spatial error model (SEM) by GeoDA using carbon storage as dependent variable and several independent variables, including diameter of living trees (DBH), number of trees per hectare (TPH), elevation (Elev), slope (Slope), and product of precipitation and temperature (Rain_Temp). Global Moran’s I was computed for describing overall spatial autocorrelations of model results at different spatial scales. Local Moran’s I was calculated at the optimal bandwidth (25 km) to present spatial distribution residuals. Intrablock spatial variances were computed to explain spatial heterogeneity of residuals. Finally, a spatial distribution map of carbon storage in Heilongjiang was visualized based on predictions. The results showed that the distribution of forest carbon storage in Heilongjiang had spatial effect and was significantly influenced by stand, topographic and meteorological factors, especially average DBH. SEM could solve the spatial autocorrelation and heterogeneity well. There were significant spatial differences in distribution of forest carbon storage. The carbon storage was mainly distributed in Zhangguangcai Mountain, Xiao Xing’an Mountain and Da Xing’an Mountain where dense forests existed, rarely distributed in Songnen  Plains, while Wanda Mountain had moderatelevel carbon storage.