Abstract: Micro-reclamation technology had become a new research direction of land reclamation in coal mining areas recently. The man-made ecosystem had been established after micro-reclamation for several years,the spatial heterogeneity of soil nutrients under shrubbery in the micro-reclamation demonstration base had been paid close attention to the further research. The microbial reclamation area in the coal mining subsidence area of Daliuta,Shendong,was selected as the research area. The total nitrogen,organic matter,and available phosphorus in shrub forests at microbial reclamation area in the coal mining subsidence area were studied using classical statistics,geostatistics,and coordinated kriging interpolation methods to obtain the spatial heterogeneity and its influencing factors. The results showed that the average values of soil total nitrogen,organic matter,and available phosphorus content in the study area were 0. 50 g / kg, 11. 93 g / kg and 2. 24 mg / kg respectively,all of which exhibited a moderate-intensity variation and spatial autocorrelation. The soil mycorrhizal hyphal density had been chosen as the key factors of the arbuscular mycorrhizal fungi characteristics. The results showed that the hyphal density has very significant positive correlations with soil total nitrogen,organic matter and available phosphorus. The best fitting models for the cooperative variation functions of soil total nitrogen,organic matter,and available phosphorus with hyphal density as an auxiliary variable are gaussian model,spherical model and gaussian model,and their coefficients of determination are 0. 709,0. 757 and 0. 467,which are higher than those of single-variable soil total nitrogen,organic matter,and available phosphorus. The soil total nitrogen,organic matter and available phosphorus all showed strong spatial autocorrelation with hyphal density as auxiliary variable,which can significantly improve the prediction accuracy of soil total nitrogen,organic matter and available phosphorus content. The spatial distribution of soil total nitrogen was block-shaped,and had a strong dependence on space. Soil available phosphorus was irregularly patchy and highly fragmented. Cooperative Kriging interpolation with mycelium density as an auxiliary variable can significantly improve the prediction accuracy of soil total nitrogen organic matter and available phosphorus. There are different spatial heterogeneity and spatial correlations of soil total nitrogen and organic matter in shrub forests at the microbial reclamation zone of coal mining subsidence area.