Abstract: The distribution patterns of the Pinus tabulaeformis population in the National Nature Reserve of Lingkong Mountain in Shanxi Province were studied using geostatistics. According to the data of 4 hm² sample plots, the semivariances of the P. tabulaeformis diameter grades were calculated, semivariogram models were simulated, and Kriging interpolation was performed to draw Kriging maps showing different diameter grades. Results showed that the semivariograms of diameter gradesⅠ, Ⅲ and Ⅳ well fit the exponential model, while diameter grade Ⅱ fit the spherical model, suggesting aggregated spatial distribution patterns in all four diameter grades. In addition, all diameter grades exhibited medium spatial correlation with the range of spatial autocorrelation of 23.4 m, 15.2 m, 11.1 m and 24.9 m, respectively, and the fractal dimensions were ranked as: diameter grade Ⅱ (1.999) > diameter grade Ⅲ (1.995) > diameter grade Ⅳ (1.973) > diameter grade Ⅰ (1.969). The semivariances of height were basically the same as the semivariances of diameter at breast height (DBH), which further verified the accuracy of our analysis. The Kriging map drawn by the Surfer software demonstrated the patch aggregation of the P. tabulaeformis population. The semivariance analysis and Kriging interpolation methods filled the loopholes of traditional population spatial distribution analyses, and were advantageous in revealing the changes in spatial pattern and aggregate intensity for more precise studies of plant patterns in the landscape.