Ecological scholars have shown spatial scale dependence with biodiversity. However, because there is no single scale appropriate for measuring variation in resource distributions, the relationship between spatial distributions of limiting resources and patterns of plant diversity has eluded plant ecologists. Thus, a multiscale method is imperative for understanding spatial distribution patterns. The objectives of this study were: (1) to quantify species diversity with selected species indices at different spatial scales, (2) to analyze the spatial distribution patterns of species diversity for the same spatial scales using multifractal theory (a multiscale, spatial structure analytical tool) and (3) to determine relationships between them. A 100m×100m(1hm2.) plot, located in TianMu Mountain National Natural Preservation Areas, in the western of Zhejiang Province, was investigated by precision measuring instruments. For spatial pattern analysis of species based on multiscale idea, the plot was divided into 10 subsets of 10 different spatial scales, namely 10m×10m, 20m×20m……100m×100m. In each scale, species diversity indices, including the Shannon-Wiener index (H), the Margalef index (K), and the Evenness index (E), were first calculated, and then multifractal parameters, such as the singularity index (Lipshitz-Hlder exponents) α and its fractal dimensions f(α) (multifractal spectrums), were computed. All visual algorithm programs were designed with the MatLab language (version 6.5 and above). Then, variations of species diversity and spatial distribution patterns for different scales were analyzed with linear and nonlinear regression and their relationships developed. Results of nonlinear regression analysis with power functions fit to curves of increasing spatial scale versus diversity indices showed strong relationships with H increasing (R2 = 0.71) but E (R2 = 0.89) and K (R2 = 0.88) decreasing indicating that species diversity existed with scale dependence. From the multifractal method, results indicated that species spatial distribution had multifractal features. For instance, through spatial scale changes the multifractal spectrum of species showed a large difference with f(α) being the subset of species with singularity α. So, an infinite hierarchy of dimensions using multifractal formality can characterize the internal spatial pattern. Also, when the spatial scale increased, αmin decreased while αmax increased meaning that large-scale spatial patterns of species were aggregated compared with small-scale patterns. In addition, the f(α)-α spectrum range (SR) and its symmetrical (Dist) in different scales showed that larger spatial scales had a wider SR and less symmetry meaning the spatial distribution of species for larger scales was nonuniform. Moreover, for the diversity indices H, E, and K, strong linear relationships with αmin (R2 = 0.67-0.94) and clear nonlinear associations with power functions for SR (R2 = 0.58-0.89) and Dist (R2 = 0.62-088) were found. Since interactions of species diversity and spatial characteristics are very complex, the above relationships need further validation along with precise explanations of any correlations among their ecological processes. Overall, the importance of spatial scale in fractal theory, patterns, processes, and scales can be linked to fractals in ecology with some innovative achievements occurring when traditional research methods are combined with fractal theory.