Next Article in Journal
Use of Diabetes Treatment Satisfaction Questionnaire in Diabetes Care: Importance of Patient-Reported Outcomes
Next Article in Special Issue
A New TS Algorithm for Solving Low-Carbon Logistics Vehicle Routing Problem with Split Deliveries by Backpack—From a Green Operation Perspective
Previous Article in Journal
Patient-Related Factors Associated with Adherence to Recommendations Made by a Fracture Liaison Service: A Mixed-Method Prospective Study
Previous Article in Special Issue
Will Happiness Improve the Psychological Integration of Migrant Workers?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 15
Issue 5
10.3390/ijerph15050946
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Pattern and Population Structure of Artemisia ordosica Shrub in a Desert Grassland under Enclosure, Northwest China

by
Jiankang Liu
1 and
Kebin Zhang
1,2,*
1
School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China
2
Key Laboratory of State Forestry Administration on Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(5), 946; https://doi.org/10.3390/ijerph15050946
Submission received: 4 April 2018 / Revised: 2 May 2018 / Accepted: 2 May 2018 / Published: 9 May 2018
(This article belongs to the Special Issue Green Environment, Green Operations and Sustainability)
Browse Figures
Versions Notes

Abstract

:
Enclosure is an effective practice for restoring and rehabilitating the degraded grassland ecosystem caused by overgrazing. Shrub species, which are dominant in most desert grasslands in arid and semiarid regions, have some beneficial ecological functions for grassland restoration. However, how the population structure and spatial pattern of the Artemisia ordosica shrub changes in a grassland ecosystem under enclosed practice is not well understood. This study, conducted in the Mu Us desert in northwest China, was designed to measure the A. ordosica population according to the chronosequence of enclosure (enclosure periods ranged from 5 years, 10 years, 15 years, and 25 years), contrasting this with an adjacent continuously grazed grassland. The results showed that the enclosed grasslands had a higher number of individuals of different age classes (seedling, adult, aging, and dead group) and greater population coverage, but shrubs had significant lower (p < 0.05) crown diameter and height in comparison with those in continuously grazed grassland. Further, enclosed grasslands had a significantly higher (p < 0.05) Shannon-Wiener index (H) and Evenness index (E), but a significantly lower (p < 0.05) Richness index (R) than continuously grazed grassland. The crown of A. ordosica showed a significant linear positive correlation with height in all plots across succession, indicating that it was feasible to analyze the age structure by crown. The crown-class distribution structure of the A. ordosica population approximated a Gaussian distribution model in all survey plots. Within the population, seedling and adult groups exhibited aggregated spatial distribution at small scales, while aging and dead A. ordosica groups showed random distribution at almost all scales in different plots. The seedling A. ordosica group showed a positive correlation with adults at small scales in all plots except in 10 years of enclosure. However, it showed independent correlation with aging and dead groups at almost all scales. In long-term enclosed plots, the mortality rate of the A. ordosica population increased, therefore assistance management practices, such as fertilization, mowing, interval grazing, and seasonal grazing, must be employed to maintain population stability after long-term enclosure. This study can improve understanding and clarify the effects of enclosures in the desert grasslands of northwest China.

1. Introduction

Grassland, which occupies approximately 37% of the earth’s terrestrial surface [1], is one of the world’s largest vegetation types. It is an important and essential natural ecosystem, which plays a key role in biodiversity conservation, socioeconomic development, greenhouse gas balance, food supply, and regional ecosystem stability [1,2,3]. However, due to climatic variation (global warming, droughts) [4,5,6] and human activities (plowing, deforestation, overcutting, and overgrazing) [7,8,9], grassland ecosystem degradation (characterized by the reduction of vegetation coverage, productivity, biodiversity and ecosystem services, and soil nutrient degradation) has become a major environmental issue all over the world in the recent decades [10,11,12,13,14]. Previous studies indicated that 49.3% of the world’s grasslands has been degraded to different degrees of severity (cf. 31.8% in China) [4,7]. Mu Us desert, as an important ecological barrier in northwest China, plays an important role in the national protection of wind erosion and sand storms [15]. However, due to increasing anthropogenic activities in this region, the grassland ecosystem has seriously deteriorated and become a crucial problem to be solved in recent decades [16,17].
Although a variety of practices, such as reseeding, plowing, fertilizing, establishing artificial grasslands, controlling grazing intensity, and rational management of the grassland, have been implemented to restore degraded grassland in Mu Us desert, enclosure has been widely considered as the most useful measure for low cost, simple implementation and quick effects [18,19,20]. Enclosure can increase plant biomass (above- and below-ground biomass) [21,22], improve plant diversity and soil physicochemical properties [19,23,24], promote ecosystem resilience [25] through livestock exclusion, and improve the livelihood of environmental conditions [26,27].
Vegetation restoration of degraded ecosystems is the joint action progress of plant–plant interactions [28,29], plant–soil feedback responses [30,31], and the facilitation of or interference between natural communities and environment [32,33]. The population structure and spatial pattern are the specific performances of comprehensive interactions among intra- and inter-specific interactions, species characteristics, interference, and environmental factors over the long term [34,35,36,37]. The population structure not only reflects the configuration and population dynamics of individuals, but also provides an important scientific basis for predicting the succession and development trend of populations [37,38]. The spatial pattern and correlation of species which represent the interaction between plants significantly affect the recruitment, growth, death, and resource utilization of species [36], and might depend on biological characteristics of plants and environmental heterogeneity [35,39]. Those two features can be used for dynamic analysis of vegetation [40].
Shrub species, which are dominant in most desert grasslands in arid and semiarid regions [41,42], are the crucial element in rehabilitating a degraded grassland ecosystem. They can provide shelter for other species [43,44] and significantly improve seedling survival for its “nursing” feature in arid and semiarid regions [45,46], eventually facilitating more complex vegetation restoration. Artemisia ordosica, a native subshrub species, is the main dominant species widely distributed in the Mu Us desert [15]. Therefore, the study on the A. ordosica shrub species is important for understanding the progress of vegetation restoration in this study area. Recent research on A. ordosica population dynamics mostly concentrated on natural recruitment [15,47], effects of habitat changes [15,48,49], spatial pattern [50,51], plant morphology and structural characteristics [52], soil moisture characteristics [53,54], and soil physicochemical properties [55]. However, the dynamics of the population structure and spatial pattern of the A. ordosica population in a restored desert grassland after removing livestock have not received enough attention.
To better understand and clarify the effects of enclosure on the A. ordosica population in the southeast Mu Us desert, different periods of enclosure and one comparable, continuously grazed grassland were selected to analyze spatial pattern and population structure. The objectives of this study were to: (1) assess how population structure and spatial pattern of A. ordosica change after enclosure and (2) reveal the recruitment mechanism of the A. ordosica population across the succession time. The main findings of this study aim to provide the scientific basis and theoretical support for grassland restoration in Mu Us desert and to provide policy guidance for the restoration and management of degraded grasslands in other arid and semiarid regions around the world.

2. Materials and Methods

2.1. Study Site

The study was conducted in an artificially fenced area (37°50′ N, 107°23′ E) located in the Mu Us desert, northwest China. The elevation is approximately 1395 m. The local climate is classified as temperate continental, characterized by frequent drought, low precipitation, strong evaporation, short and hot summers, and long and cold winters. The annual precipitation mainly ranges from 250 mm to 350 mm. More than 80% of the precipitation occurs in the growing season of plants (from June to September). The mean annual evapotranspiration is 2024 mm. The annual mean air temperature is 8.1 °C, with the highest (16.2 °C) and lowest (−14.4 °C) monthly mean temperatures in July and January, respectively. The average frost-free period is about 165 days. The annual mean wind speed is 2.8 m/s, but bouts of wind over 6 m/s are common and dominated by west and northwest winds. The dominant species in the study area are some shrubs, semishrubs, or perennial herbs, such as Artemisia ordosica, Heteropappus altaicus, Caragana Korshinskii, Artemisia scoparia, and Salix psammophila. A. ordosica is the main focus of study area.
A. ordosica is a long-lived (about 10 years), deciduous, dwarf perennial shrub with linearly lobate leaves and is widely distributed in Mu Us desert, China [15,56,57]. It is easier to be grazed by livestock in the winter and spring seasons when forage grass is absent, although the palatability is poor during the growth season of plants. The flowering season mainly occurs around August and lasts approximately 50 days. The peak growing season is around July. A. ordosica, a low-growing (usually less than 100 cm) shrub species, has some paratactic branches near the ground [50,58]. Its root system is mainly distributed in the 30 cm of upper sand [47,48]. A. ordosica is generally recruited through seeds [59], which are tiny, light, and susceptible to wind-dispersion of more than several miles [60].

2.2. Study Design

Owing to excessive grazing by local Ning-Xia Tan sheep, local vegetation in the study area has been degraded seriously since late 20th century. Enclosure is the main practice implemented to restore the ecological system and the livelihood of the environment. We encircled the grassland with a chain link fence to remove herbivores such as local Ning-Xia Tan sheep, which are generally grazed by herdsman. However, there were still many small wild herbivores in enclosed grassland, such as rabbit and mice. The space-for-time method was used to monitor the effects of enclosure on the A. ordosica population. In this study, four successional ages, including enclosed 5 years (F5), 10 years (F10), 15 years (F15), and 25 years (F25), contrasting with adjacent continuously grazed grassland (CG), were selected for the plant survey (Table 1). Before enclosure establishment, the study area was grazed under similar intensity ensuring uniform natural conditions over the whole relatively flat study site. The five selected plots were contiguous. All plots were in the same continuous flat area and have the same slope, altitude, and soil type, insuring that their edaphic conditions were comparable before enclosure—although these conditions later changed after enclosure (Table 1). The vegetation coverage of the enclosed grassland was higher than that of continuously grazed grassland. All vegetation surveys were conducted in mid-August 2017, the peak growing season of plants.
Data were collected on five 0.25 ha plots (50 × 50 m), each of which was established randomly for vegetation survey at F5, F10, F15, F25, and CG, respectively. A. ordosica was the dominant species in these five sample plots. Each sample plot was divided into 100 contiguous 5 × 5 m subplots. In each subplot, relative location coordinates, height (H), wide (WC) and narrow (NC) crown diameters, and the health status of each A. ordosica was recorded. The crown diameter of the shrub was calculated as WC/2 + NC/2 [61,62], like many other studies [15,49,63,64]. The height of the shrub was measured at the center of the crown where it is maximum [63,64]. All individuals, including seedling, adult, aging, and dead groups were crown-mapped (Figure 1). Further, in each selected plot, 10 random quadrats (1 × 1 m) were selected for an herb layer survey. In each quadrat, the density of each plant species was recorded.

2.3. Data Analysis

2.3.1. Population Structure

A. ordosica, a subshrub species, does not produce notable annual rings [51]. Therefore, the morphological characteristics were used to distinguish age structures. In this study, crown classes were used to explain the age structures of A. ordosica population [65]. All A. ordosica individuals were classified into different classes (Class 1: from 0.1 cm to 10 cm; 10 cm were added gradually from Class 2 to 4 age groups (seedling group: average heights and crown were less than 30 cm; adult group: average heights and crown were more than 30 cm, and the dead branches proportion was less than two-thirds; aging group: the dead branches proportion was more than two-thirds; and dead group: individual which was lifeless) which were consistent with previous research [15,51,65,66]. The relationship between crown and height was fitted by a linear model. A Gaussian model was used for the curve fitting of the crown distribution structure to analyze the population dynamics.

2.3.2. Spatial Pattern Analysis

The univariate spatial autocorrelation and intraspecific spatial correlation of different age groups were analyzed at multiple scales by the pair-correlation function g(r) instead of the commonly used Ripley’s K-function in this study. The g-function, which is related to the derivative of the K-function, is more sensitive to small-scale effects in comparison with K-function [67,68,69]. The formula of function g(r) is as follows [39,70]:
g ( r ) = 1 2 π t A 2 n 2 ∑ i = 1 n ∑ i = 1 j ≠ 1 n w i j − 1 k h ( t − | x i − x j | )
where A represents the plot area, n is the total plants number, and wij means a weighting factor correcting for edge effects. kh, a kernel function, is used for applying maximum weight to point pairs within a distance t.
The g-function, based on point-to-point distances analysis, describes aggregating, random, or regularity spatial distribution within a given radius r by using a standardized density [71]. At a given distance r, g(r) > 1 indicates an aggregation trend, g(r) = 1 indicates complete spatial randomness (CSR), and g(r) < 1 indicates a regularity trend. To test significance departures of CSR, the 5th lowest and 199th highest values from the Monte Carlo simulations approach was used to generate approximately two-sided 95% simulation envelopes after calculating for each distance r of a point process.

2.3.3. Species Diversity Index

The Shannon-Wiener (H), Richness (R), and Evenness indices (E) were used to illustrate species diversity [25]. These species diversity indices were calculated as:
Richness index (R):
R = S
Shannon-Wiener index (H):
H = − ∑ i = 1 S ( P i l n P i )
Evenness index (E):
E = H/lnS
where S is the total species number found in each quadrat and Pi is the density proportion of ith species.
All descriptive statistical parameters were carried out using EXCEL 2013. The significance test (p < 0.05) was calculated by SPSS 20.0 software (SPSS for Windows, Chicago, IL, USA). Spatial analysis was completed by Programita Febrero 2014. All figures were plotted by OriginPro 2015 (OriginLab Corporation, Northampton, MA, USA).

3. Results

3.1. Population Structure of A. ordosica Populations

The data showed that enclosure could significantly influence the individuals of different groups, coverage, crown, and height of the A. ordosica population (Table 2, Figure 1). Enclosed plots (F5, F10, F15, and F25) had more individuals of different age classes (seedling, adult, aging, and dead group), and greater population coverage, with significantly lower crown and height (p < 0.05) than the continuously grazed plot (CG). With the increase of enclosure time, the number of adult and total A. ordosica increased initially, and then decreased. Both parameters were highest in the F10 plot (3523 adult individuals, 4643 total individuals). In contrast, the proportion of aging (34.1%) and dead (13.4%) groups in the F10 plot, and the proportion of the seedling group (19.1%) in the F5 plot was higher than the other four plots. The proportion of adults was minimal in the F25 plot (45.0%), while maximal in the CG plot (87.0%). The F5 plot had the lowest population coverage (19.7%) and height (0.289 ± 0.102) of the four enclosed plots. The crown and height of the F10 plot were significantly higher (p < 0.05) than the other three enclosed plots (F5, F15, and F25). The F15 plot had minimum crown (0.556 ± 0.221) (Table 2, Figure 1).
There was a significant linear positive correlation (p < 0.001) between crown and height of the A. ordosica population in all plots. The regression slope (0.346) of linear equation in continuously grazed plot was higher than that in the enclosure plots (F5: 0.263, F10: 0.234, F15: 0.251, and F25: 0.310). The slope of the linear equation initially declined in the first 10 years after enclosure and subsequently slowly increased (Figure 2).
According to the polynomial fitting results (Figure 3), the crown distribution structure of the A. ordosica population approximated a Gaussian distribution model in all survey plots (R2 > 0.9, p < 0.001). Probably due to differences in the site environmental conditions after long-term enclosure, A. ordosica grew smaller in the F15 and F25 plots (up to crown = 140 cm) than in the F5 and F10 plots. The main distribution of crown is within 30–120 cm. Seedling and medium-sized (crown within 40–100 cm) A. ordosica thrived in the F5, F10, and F15 plots, which suggests a stable, renewable population. The dead and aging individuals of medium-sized A. ordosica were larger than seedlings and big-sized adults (crown more than 100 cm). At the F25 plot, the alive rate of individuals was lower than other survey plots (Figure 3).

3.2. Spatial Autocorrelation of A. ordosica Populations

The univariate spatial autocorrelation and intraspecific spatial correlation of the entire population, S (seedling groups), AD (adult groups), AG (aging groups), and D (dead groups) were unlikely in five survey plots (Figure 4, Figure 5 and Figure 6). Entire population, seedling, and adult groups exhibited aggregated distribution at small scales, especially at scales of 0–5 m (Figure 4 and Figure 5). The aging and dead A. ordosica groups showed random distribution at almost all scales in the F5, F10, F25, and CG plots, while those two groups showed aggregated distribution within 0–11 m and 12–15 m in the F15 plot (Figure 5). The aggregation degree of A. ordosica adult groups was greater than that of entire and seedling populations in the F10, F15, and CG plots, but was less in the F5 and F25 plots (Figure 4 and Figure 5). Further, the aggregation degree of A. ordosica seedling groups declined with the distances increasing in all plots (Figure 5). The seedlings of A. ordosica showed positive correlations with adults at small scales in all plots except for F10, especially at scales of 0–5 m, although they showed independent correlations at almost all scales with aging and dead A. ordosica (Figure 6).

3.3. Species Diversity of Community

In this study, enclosure grasslands (F5, F10, F15, and F25) had significantly higher (p < 0.05) Shannon-Wiener (H) and Evenness (E) indices but had a significantly lower (p < 0.05) Richness index (R) than continuously grazed grassland (CG). Further, with the increase of enclosure time, H and E indices were all increased initially, and then decreased, which were highest in the F10 plot and the F15 plot, respectively. However, the R index decreased initially and then increased with the increase of enclosure time, which was lowest in the F10 plot (Figure 7).

4. Discussion

4.1. Population Structure of A. ordosica Population

Population structure is the synthetic action of biological characteristics, environmental factors, and intraspecific interactions [40,72] which can be used to reflect the dynamics and development trends of a population [37]. Some essential characteristics of plants, such as density, height, crown, and diameter at breast height (DBH), can reflect the growth status and ecological functions of plants and affect the ecological process of a population, community, or even an ecosystem [73,74,75]. With the gradual restoration of an ecological environment after enclosure, the vegetation community gradually recovered [23,24,26,27], characterized by a significant increase (p < 0.05) of Shannon-Wiener and Evenness indices and rapidly growing development of the A. ordosica population (Table 2, Figure 1 and Figure 7). In the early stage of restoring grassland, natural conditions are favorable for the recruitment of A. ordosica because of adequate resource supply and lack of competition from other species [42,71]. However, the massive reproduction of A. ordosica restricted the growing development of plant species with weak competitiveness. As shown in this study, for individuals and coverage of the A. ordosica population, Shannon-Wiener and Evenness indices increased quickly, but the Richness index greatly reduced in the first 10 years after enclosure. At this time, a great deal of seeds germinated and grew into seedlings or adults, which caused a high reduction of shrubs’ height and crown (Table 2). However, elevated plant density and size lead to some individuals dying or withering because of a lack of ability to face intense intraspecific competition for light, water, and nutrients [76,77,78]. Further, the appearance and development of soil crusts, which is water-absorbing and hinders water infiltration, would prevent seeds from germinating [79,80,81]. Therefore, for the individual seedlings and adults, the population coverage, height, crown diameter, Shannon-Wiener index, and Evenness index decreased, while the proportion of aging and dead individuals of different-sized A. ordosica increased after being enclosed for 25 years (Table 2, Figure 3 and Figure 7). This reduction of dominant species reduced competitive pressure on other species [25], so the Richness index increased after being enclosed for 25 years (Figure 7).
The allometric model of plant characteristics can reflect the plastic responses of plants to resource availability and the natural selection or ecological sorting of species in the competitive environment of local resources [82,83]. This is the long-term response of the interaction between plants and environment [84]. As shown in Figure 2, the regression model indicated how plant height changed with crown. The data showed that the linear equation slope in enclosure plots was lower than the continuously grazed plot (Figure 2), which was probably caused by the removal of livestock activities such as ingestion and trampling [85]. The significant positive linear correlation between height and crown of A. ordosica indicated that it was feasible to analyze age structure by crown distribution.
As a dominant shrub species in Mu Us desert, A. ordosica adapts well to the local environment [15]. Our study showed that the crown distribution structure of the A. ordosica population approximated a Gaussian distribution model (Figure 3). In addition, the abundance of seedling and medium-sized A. ordosica demonstrated a stable crown-diameter distribution structure in enclosed grassland (Figure 3) [40]. The exclusion of livestock had little effect on crown-diameter distribution structure across the succession time, but it could influence the recruitment and mortality of the population (Figure 3).

4.2. Spatial Pattern of A. ordosica Population

The spatial pattern of plants is the performance of comprehensive interactions between plants and environmental factors over a long term, which can be used to reveal the succession law of vegetation [35,36,37,40,86,87]. The spatial-patterns structure analysis of different growth stages can obtain some essential information of a population under natural or human disturbance, such as recruitment, wilting or mortality of individuals, and intraspecific competition [51,70,88,89]. In this study, different age A. ordosica groups showed unpredictable spatial patterns in enclosed and grazing grasslands. The entire A. ordosica population, seedling, and adult groups showed significantly aggregated trends at small scales but random or even a uniform trend at large scales in all survey plots (Figure 4 and Figure 5), which was consistent with previous studies in this area [15,51]. This is mainly determined by plant–plant interactions [90] and patches resources [91]. Moreover, because of the limited resources and spaces for reproduction and expansion, aggregated distribution at small patches is beneficial for individuals, especially for seedlings, to resist adverse nature conditions [15,42].
The spatial distribution pattern of plants is closely related to human or animal disturbance [88,89,92], such as grazing, which is widely considered to be one of the most direct and essential factors for vegetation succession [25,93,94]. Removal of livelihood usually led to the reduction of environmental spatial heterogeneity [95]. Therefore, the aggregation degree of the A. ordosica population in the F10 plot was lower than that in the F5 and CG plots. Further, the random mortality pattern of plants [96] results in a random distribution of aging and dead groups (Figure 5). With the ecological environment gradually improved after the grassland was enclosed, the A. ordosica population rapidly developed (Table 2) due to the reduction of environmental stress (i.e., sand burial, wind erosion, water erosion) [26,27,51]. The elevated shrub density caused nonrandom intraspecific competition [97,98], characterized by an aggregation of dead and aging groups in the F15 plot (Figure 5). Therefore, with the continuation of enclosure, the surviving A. ordosica showed an even more aggregated spatial autocorrelation pattern in the F25 plot (Figure 5 and Figure 6).
It is likely that adult groups were relatively important for A. ordosica population recruitment in this study (Figure 6). The seed of A. ordosica is tiny and light and has some special structures to allow airborne dispersal [60]. The seeds are more likely to be caught by adult shrubs, and then generate and survive under or nearby the canopy because of the amelioration of the microenvironment [99,100,101]. Moreover, large shrub patches can provide a shelter for seedlings under canopy and protect the seedlings from water–wind erosion and sand burial [102,103,104]. However, when shrubs are wilted or dead, the nursing effects will diminish [105]. Spatial pattern analysis could well support this theory. From the results of this study, A. ordosica seedlings had an aggregated distribution within 0–5 m and had a significant positive correlation with adults at small scales (0–5 m) at the beginning of grassland enclosure (Figure 5 and Figure 6). The seedling group showed independent correlation with aging and dead groups (Figure 6). Soil nutrients were gradually enriched under shrub canopy after enclosure because they intercepted atmospheric dust and litter decomposition [106]. Therefore, the soil crust cover and thickness under shrub canopy are larger than in open spaces nearby [107], which limit the germination of seeds [100]. So, after enclosure for 10 years, the spatial correlation between seedlings and adults became independent (Figure 6). At this time, limited resources result in density-dependent self-thinning, although spatial heterogeneity (i.e., microtopography, soil, and water) causes nonrandom mortality of plants [40,71]. Therefore, some individuals gradually died because of increasing competition, leading to a reduction of population density (Table 2, Figure 3 and Figure 4). Finally, after 15 years of enclosure, the seedlings were positively associated with adults at small scales (spatially) and changed to an independent correlation at large scales (Figure 6).

5. Conclusions

Enclosure is an effective practice to recover the A. ordosica population in Mu Us desert. It can increase Shannon-Wiener and Evenness indices while decreasing the Richness index. Moderate enclosure had positive effects on the A. ordosica population in degraded grazing grassland, including an increase in the seedling ratio, population individuals, population coverage, and reduction in the height–crown ratio. However, the aging and death rate of the population increased after a long period of enclosure. The entire population, seedling, and adult A. ordosica groups showed an aggregated distribution at small scales, although aging and dead groups were randomly distributed at almost all scales in different plots. The seedling A. ordosica group showed a significant positive correlation with adults, while independent correlation with aging and dead groups. Enclosure is an important factor affecting the spatial autocorrelation and intraspecific spatial correlation of alive A. ordosica. With the continuation of enclosure, the surviving A. ordosica showed a more aggregated spatial autocorrelation pattern. Further, the spatial correlation between seedlings and adults changed from positive at small scales in the F5 plot to independent at all scales in the F10 plot and eventually became more positive at small scales. Those results indicated that continuous long-term enclosure might cause re-degradation of the A. ordosica population. Therefore, other artificial measures, such as fertilization, mowing, interval grazing, and seasonal grazing, should be implemented to maintain the stability and sustainability of shrub-dominated grasslands after restoration.

Author Contributions

J.L. and K.Z. conceived, designed and conducted the experiments. J.L. analyzed the data and wrote the manuscript. K.Z. revised the manuscript.

Funding

This research was funded by the National Key Technology R&D Program (2016YFC0500908), the National Natural Science Foundation of China (30771764), and the National Forestry Administration Desertification Positioning Monitoring Project (660550).

Acknowledgments

We thank Li Ma, Shuqin Liu, Hongda Li, Rui Ma, Yan Zhao, Mengwan Wu, and Haiying Li who collected the data of this study. Many thanks to Wanxue You, Cunlin Bai, Zhilin Huang, and other staff of Habahu Nature Reserve of China for support field survey. We are grateful to the editor and anonymous reviewers for their helpful comments on improving previous manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. O’Mara, F.P. The role of grasslands in food security and climate change. Ann. Bot. 2012, 110, 1263–1270. [Google Scholar] [CrossRef] [PubMed]
  2. Henzell, E.F. Perspectives in Grassland Ecology. Agro-Ecosystems 1981, 7, 169–170. [Google Scholar] [CrossRef]
  3. Zheng, Y.; Xie, Z.; Jiang, L.; Wu, Y.; Shimizu, H. Model simulation and comparison of the ecological characteristics of three degraded grassland types in China. Belg. J. Bot. 2005, 138, 109–118. [Google Scholar]
  4. Wang, Z.; Deng, X.; Song, W.; Li, Z.; Chen, J. What is the main cause of grassland degradation? A case study of grassland ecosystem service in the middle-south Inner Mongolia. Catena 2017, 150, 100–107. [Google Scholar] [CrossRef]
  5. Sivakumar, M.V.K. Interactions between climate and desertification. Agric. For. Meteorol. 2007, 142, 143–155. [Google Scholar] [CrossRef]
  6. John, R.; Chen, J.; Ou-Yang, Z.T.; Becker, R.; Samanta, A.; Batkhishig, O. Vegetation response to extreme climate events on the Mongolian Plateau from 2000 to 2010. Environ. Res. Lett. 2013, 8, 035033. [Google Scholar] [CrossRef]
  7. Gang, C.; Zhou, W.; Chen, Y.; Wang, Z.; Sun, Z.; Li, J.; Odeh, I. Quantitative assessment of the contributions of climate change and human activities on global grassland degradation. Environ. Earth Sci. 2014, 72, 4273–4282. [Google Scholar] [CrossRef]
  8. Qi, J.; Chen, J.; Wan, S.; Ai, L. Understanding the coupled natural and human systems in Dryland East Asia. Environ. Res. Lett. 2012, 7, 015202. [Google Scholar] [CrossRef]
  9. Wang, H.; Li, Z.; Gong, Y.; Wang, Z.; Huang, D. Single irrigation can achieve relatively high production and water use efficiency of Siberian wildrye grass in the semiarid agropastoral ecotone of North China. Agron. J. 2009, 101, 996–1002. [Google Scholar] [CrossRef]
  10. Andrade, B.O.; Koch, C.; Boldrini, I.I.; Vélezmartin, E.; Hasenack, H.; Hermann, J.M.; Overbeck, G.E. Grassland degradation and restoration: A conceptual framework of stages and thresholds illustrated by southern Brazilian grasslands. Nat. Conserv. 2015, 13, 95–104. [Google Scholar] [CrossRef]
  11. Yang, Y.; Wang, Z.; Li, J.; Gang, C.; Zhang, Y.; Zhang, Y.; Qi, J. Comparative assessment of grassland degradation dynamics in response to climate variation and human activities in China, Mongolia, Pakistan and Uzbekistan from 2000 to 2013. J. Arid Environ. 2016, 135, 164–172. [Google Scholar] [CrossRef]
  12. Jiao, F.; Wen, Z.M.; An, S.S. Changes in soil properties across a chronosequence of vegetation restoration on the Loess Plateau of China. Catena 2011, 86, 110–116. [Google Scholar] [CrossRef]
  13. Handlová, V.; Münzbergová, Z. Seed banks of managed and degraded grasslands in the Krkonoše Mts., Czech Republic. Folia Geobotanica 2006, 41, 275–288. [Google Scholar] [CrossRef]
  14. Henning, J.A.G.; Kellner, K. Degradation of a soil (aridosol) and vegetation in the semi-arid grasslands of southern Africa. Bot. Bull. Acad. Sin. 1994, 35, 195–199. [Google Scholar]
  15. Wu, B.; Yang, H. Spatial patterns and natural recruitment of native shrubs in a semi-arid sandy land. PLoS ONE 2013, 8, e58331. [Google Scholar] [CrossRef] [PubMed]
  16. Wu, B.; Ci, L.J. Landscape change and desertification development in the Mu Us Sandland, Northern China. J. Arid Environ. 2002, 50, 429–444. [Google Scholar] [CrossRef]
  17. Wang, X.; Chen, F.H.; Dong, Z.; Xia, D. Evolution of the southern Mu Us Desert in north China over the past 50 years: An analysis using proxies of human activity and climate parameters. Land Degrad. Dev. 2005, 16, 351–366. [Google Scholar] [CrossRef]
  18. He, Y.H.; Zhao, H.L.; Liu, X.P.; Zhao, X.Y.; LI, Y.Q. Soil physical and chemical characteristics of sandy meadow in natural restoration process. J. Soil Water Conserv. 2008, 22, 159–161. [Google Scholar] [CrossRef]
  19. Zeng, Q.; Liu, Y.; Xiao, L.; Huang, Y. How Fencing Affects the Soil Quality and Plant Biomass in the Grassland of the Loess Plateau. Int. J. Environ. Res. Public Health 2017, 14, 1117. [Google Scholar] [CrossRef] [PubMed]
  20. Zhou, H.; Zhao, X.; Tang, Y.; Gu, S.; Zhou, L. Alpine grassland degradation and its control in the source region of the Yangtze and Yellow Rivers, China. Grassl. Sci. 2005, 51, 191–203. [Google Scholar] [CrossRef]
  21. Batoyun, T.; Shinoda, M.; Cheng, Y.; Purevdorj, Y. Effects of grazing and precipitation variability on vegetation dynamics in a Mongolian dry steppe. J. Plant Ecol. 2016, 9, 508–519. [Google Scholar] [CrossRef]
  22. Zhu, G.Y.; Deng, L.; Zhang, X.B.; Shangguan, Z.P. Effects of grazing exclusion on plant community and soil physicochemical properties in a desert steppe on the Loess Plateau, China. Ecol. Eng. 2016, 90, 372–381. [Google Scholar] [CrossRef]
  23. Li, W.; Cao, W.; Wang, J.; Li, X.; Xu, C.; Shi, S. Effects of grazing regime on vegetation structure, productivity, soil quality, carbon and nitrogen storage of alpine meadow on the Qinghai-Tibetan Plateau. Ecol. Eng. 2017, 98, 123–133. [Google Scholar] [CrossRef]
  24. Deng, L.; Zhang, Z.; Shangguan, Z. Long-term fencing effects on plant diversity and soil properties in China. Soil Tillage Res. 2014, 137, 7–15. [Google Scholar] [CrossRef]
  25. Wu, G.L.; Du, G.Z.; Liu, Z.H.; Thirgood, S. Effect of fencing and grazing on a Kobresia-dominated meadow in the Qinghai-Tibetan Plateau. Plant Soil 2009, 319, 115–126. [Google Scholar] [CrossRef]
  26. Shang, Z.H.; Ma, Y.S.; Long, R.J.; Thirgood, S. Effect of fencing, artificial seeding and abandonment on vegetation composition and dynamics of ‘black soil land’ in the headwaters of the Yangtze and the Yellow Rivers of the Qinghai-Tibetan Plateau. Land Degrad. Dev. 2008, 19, 554–563. [Google Scholar] [CrossRef]
  27. Gao, Q.; Li, Y.; Wan, Y.; Jiangcun, W.Z.; Qin, X.B.; Wang, B.S. Significant achievements in protection and restoration of alpine grassland ecosystem in Northern Tibet, China. Restor. Ecol. 2009, 17, 320–323. [Google Scholar] [CrossRef]
  28. Gómezaparicio, L. The role of plant interactions in the restoration of degraded ecosystems: A meta-analysis across life-forms and ecosystems. J. Ecol. 2009, 97, 1202–1214. [Google Scholar] [CrossRef]
  29. Brooker, R.W.; Maestre, F.T.; Callaway, R.M.; Lortie, C.L.; Cavieres, L.A.; Kunstler, G.; Armas, C. Facilitation in plant communities: The past, the present, and the future. J. Ecol. 2008, 96, 18–34. [Google Scholar] [CrossRef]
  30. Bezemer, T.; Lawson, C.S.; Hedlund, K.; Edwards, A.R.; Brook, A.J.; Igual, J.M.; Van Der Putten, W.H. Plant species and functional group effects on abiotic and microbial soil properties and plant–soil feedback responses in two grasslands. J. Ecol. 2006, 94, 893–904. [Google Scholar] [CrossRef]
  31. Zuo, X.; Zhao, X.; Zhao, H.; Zhang, T.; Guo, Y.; Li, Y.; Huang, Y. Spatial heterogeneity of soil properties and vegetation–soil relationships following vegetation restoration of mobile dunes in Horqin Sandy Land, Northern China. Plant Soil 2009, 318, 153–167. [Google Scholar] [CrossRef]
  32. Forup, M.L.; Henson, K.S.E.; Craze, P.G.; Memmott, J. The restoration of ecological interactions: Plant–pollinator networks on ancient and restored heathlands. J. Appl. Ecol. 2008, 45, 742–752. [Google Scholar] [CrossRef]
  33. Tielbörger, K.; Kadmon, R. Temporal environmental variation tips the balance between facilitation and interference in desert plants. Ecology 2000, 81, 1544–1553. [Google Scholar] [CrossRef]
  34. An, S.; Zheng, F.; Zhang, F.; Van Pelt, S.; Hamer, U.; Makeschin, F. Soil quality degradation processes along a deforestation chronosequence in the Ziwuling area, China. Catena 2008, 75, 248–256. [Google Scholar] [CrossRef]
  35. Guo, Y.X.; Kang, B.; Li, G.; Wang, D.X.; Yang, G.H.; Wang, D.W. Species composition and point pattern analysis of standing trees in secondary Betula albosinensis forest in Xiaolongshan of west Qinling Mountains. J. Appl. Ecol. 2011, 22, 2574–2580. [Google Scholar] [CrossRef]
  36. Druckenbrod, D.L.; Shugart, H.H.; Davies, I. Spatial pattern and process in forest stands within the Virginia piedmont. J. Veg. Sci. 2005, 16, 37–48. [Google Scholar] [CrossRef]
  37. Mwavu, E.N.; Witkowski, E.T.F. Population structure and regeneration of multiple-use tree species in a semi-deciduous African tropical rainforest: Implications for primate conservation. For. Ecol. Manag. 2009, 258, 840–849. [Google Scholar] [CrossRef]
  38. Zhang, D.; Yang, L.; Zhang, Z. Numeric Dynamics of Artemisia ordosica Population under Different Habitats in Southeastern Edge of the Tengger Desert. J. Desert Res. 2017, 5, 893–901. [Google Scholar] [CrossRef]
  39. Condit, R.; Ashton, P.S.; Baker, P.; Bunyavejchewin, S.; Gunatilleke, S.; Gunatilleke, N.; Lee, H.S. Spatial Patterns in the Distribution of Tropical Tree Species. Science 2000, 288, 1414–1418. [Google Scholar] [CrossRef] [PubMed]
  40. Kang, D.; Guo, Y.; Ren, C.; Zhao, F.; Feng, Y.; Han, X.; Yang, G. Population structure and spatial pattern of main tree species in secondary Betula platyphylla forest in Ziwuling Mountains, China. Sci. Rep. 2014, 4, 6873. [Google Scholar] [CrossRef] [PubMed]
  41. Li, X.R. Study on shrub community diversity of Ordos Plateau, Inner Mongolia, Northern China. J. Arid Environ. 2001, 47, 271–279. [Google Scholar] [CrossRef]
  42. Wang, Y.; Yang, X.; Shi, Z. The formation of the patterns of desert shrub communities on the Western Ordos Plateau, China: The roles of seed dispersal and sand burial. PLoS ONE 2013, 8, e69970. [Google Scholar] [CrossRef] [PubMed]
  43. King, E.G.; Stanton, M.L. Facilitative effects of Aloe shrubs on grass establishment, growth, and reproduction in degraded Kenyan rangelands: Implications for restoration. Restor. Ecol. 2008, 16, 464–474. [Google Scholar] [CrossRef]
  44. Li, Y.; Chen, J.; Cui, J.; Zhao, X.; Zhang, T. Nutrient resorption in Caragana microphylla along a chronosequence of plantations: Implications for desertified land restoration in North China. Ecol. Eng. 2013, 53, 299–305. [Google Scholar] [CrossRef]
  45. Yang, X.; Baskin, C.C.; Baskin, J.M.; Liu, G.; Huang, Z. Seed mucilage improves seedling emergence of a sand desert shrub. PLoS ONE 2012, 7, e34597. [Google Scholar] [CrossRef] [PubMed]
  46. Padilla, F.M.; Pugnaire, F.I. The role of nurse plants in the restoration of degraded environments. Front. Ecol. Environ. 2006, 4, 196–202. [Google Scholar] [CrossRef]
  47. Li, S.L.; Werger, M.J.A.; Zuidema, P.A.; Yu, F.H.; Dong, M. Seedlings of the semi-shrub Artemisia ordosica are resistant to moderate wind denudation and sand burial in Mu Us sandland, China. Trees 2010, 24, 515–521. [Google Scholar] [CrossRef]
  48. Li, S.L.; Zuidema, P.A.; Yu, F.H.; Werger, M.J.; Dong, M. Effects of denudation and burial on growth and reproduction of Artemisia ordosica in Mu Us sandland. Ecol. Res. 2010, 25, 655–661. [Google Scholar] [CrossRef]
  49. She, W.; Zhang, Y.; Qin, S.; Wu, B.; Liu, Z.; Liu, J.; Zhang, W. Habitat effect on allometry of a xeric shrub (Artemisia ordosica Krasch) in the Mu Us Desert of northern China. Forests 2015, 6, 4529–4539. [Google Scholar] [CrossRef]
  50. Yang, H.; Zhang, J.; Li, Z.; Wu, B.; Zhang, Z.; Wang, Y. Comparative study on spatial patterns of the Artemisia ordosica population in the Mu Us sandy land. Acta Ecol. Sin. 2008, 28, 1901–1910. [Google Scholar] [CrossRef]
  51. Zhang, J.; Zhang, Y.; Fan, D.; Qin, S.; Jia, X.; Wu, B.; Zhu, L. Effects of sand dune stabilization on the spatial pattern of Artemisia ordosica population in Mu Us Desert, Northwest China. PLoS ONE 2015, 10, e0129728. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, J.H.; Wu, B.; Lei, Y.K.; Li, X.M.; He, J.; Gao, D.X. Analysis of Artemisia ordosica plant morphology and structure characteristics in Mu Us Sandland. J. Southwest For. Univ. 2011, 31, 6–9. [Google Scholar] [CrossRef]
  53. Xue, P.; He, X. Physiological water strategy of Artemisia ordosica around soil threshold of drought. Front. Biol. China 2008, 3, 32–39. [Google Scholar] [CrossRef]
  54. Zhang, J.H.; Han, H.Y.; Lei, Y.K.; Yang, W.B.; Li, Y.H.; Yang, D.F.; Zhao, X.B. Correlations between distribution characteristics of Artemisia ordosica root system and soil moisture under different fixation stage of sand dunes. J. Southwest For. Univ. 2012, 32, 1–5. [Google Scholar] [CrossRef]
  55. Zhao, X.L.; He, X.D.; Xue, P.P.; Zhang, N.; Wu, W.; Li, R.; Zhao, H. Effects of soil stoichiometry of the CaCO3/available phosphorus ratio on plant density in Artemisia ordosica communities. Chin. Sci. Bull. 2012, 57, 492–499. [Google Scholar] [CrossRef]
  56. She, W.; Bai, Y.; Zhang, Y.; Qin, S.; Liu, Z.; Wu, B. Plasticity in Meristem Allocation as an Adaptive Strategy of a Desert Shrub under Contrasting Environments. Front. Plant Sci. 2017, 8, 1933. [Google Scholar] [CrossRef]
  57. Kobayashi, T.; Liao, R.T.; Li, S.Q. Ecophysiological behavior of Artemisia ordosica on the process of sand dune fixation. Ecol. Res. 1995, 10, 339–349. [Google Scholar] [CrossRef]
  58. Wei, B.; Ding, G.D.; Wu, B.; Zhang, Y.Q. Study on Windbreak Effect of Artemisia ordosica Community over Different Vegetation Coverage. Arid Zone Res. 2012, 29, 213–217. [Google Scholar] [CrossRef]
  59. Li, S.L.; Yu, F.H.; Werger, M.J.A.; Dong, M.; Zuidema, P.A. Habitat-specific demography across dune fixation stages in a semi-arid sandland: Understanding the expansion, stabilization and decline of a dominant shrub. J. Ecol. 2011, 99, 610–620. [Google Scholar] [CrossRef]
  60. Ma, J.; Li, J. A comparative study on microstructure features of seed from Artemisia plants in desert of China. J. Desert Res. 2002, 22, 586–590. [Google Scholar]
  61. Martinezmeza, E.; Whitford, W.G. Stemflow, throughfall and channelization of stemflow by roots in three Chihuahuan desert shrubs. J. Arid Environ. 1996, 32, 271–287. [Google Scholar] [CrossRef]
  62. Pugnaire, F.I.; Haase, P.; Puigdefábregas, J.; Cueto, M.; Clark, S.C.; Incoll, L.D. Facilitation and succession under the canopy of a leguminous shrub, Retama sphaerocarpa, in a semi-arid environment in south-east Spain. Oikos 1996, 76, 455–464. [Google Scholar] [CrossRef]
  63. Zhang, Z.S.; Zhao, Y.; Li, X.R.; Huang, L.; Tan, H.J. Gross rainfall amount and maximum rainfall intensity in 60-minute influence on interception loss of shrubs: A 10-year observation in the Tengger Desert. Sci. Rep. 2016, 6, 26030. [Google Scholar] [CrossRef] [PubMed]
  64. Zhang, Y.; Wang, X.; Hu, R.; Pan, Y.X. Stemflow volume per unit rainfall as a good variable to determine the relationship between stemflow amount and morphological metrics of shrubs. J. Arid Environ. 2017, 141, 1–6. [Google Scholar] [CrossRef]
  65. Wang, J.H.; Jin, H.J.; Ma, Q.L.; Zhang, D.K.; Liu, Y.J.; Xu, B.Y. Structure and Distribution Pattern of Artermisia Ordosica Population in Arid Region. J. Desert Res. 2010, 30, 534–538. [Google Scholar]
  66. Zhao, Y.; Zhang, Z.; Chen, Y.; Xu, B.X.; Chen, D. Effect of Degradation Stage of Artemisia ordosica on Morphology of Nebkhas. J. Desert Res. 2015, 35, 1136–1140. [Google Scholar] [CrossRef]
  67. Illian, J.; Penttinen, A.; Stoyan, H.; Stoyan, D. Spatial Analysis and Modelling of Spatial Point Patterns. Technometrics 2008, 47, 516–517. [Google Scholar] [CrossRef]
  68. Jiang, L.; Wang, X.; Li, L.; Shi, Z.; Yang, X. Spatial Association of Shrubs and Their Interrelation to Burrowing Site Preference of Subterranean Rodents on Dune Slope in the Otindag Sandy Land, China. Sustainability 2017, 9, 1729. [Google Scholar] [CrossRef]
  69. Wiegand, T.; A Moloney, K. Rings, circles, and null-models for point pattern analysis in ecology. Oikos 2004, 104, 209–229. [Google Scholar] [CrossRef]
  70. Wang, X.; Zhang, B.; Zhang, K.; Zhou, J.; Ahmad, B. The Spatial Pattern and Interactions of Woody Plants on the Temperate Savanna of Inner Mongolia, China: The Effects of Alternating Seasonal Grazing-Mowing Regimes. PLoS ONE 2015, 10, e0133277. [Google Scholar] [CrossRef] [PubMed]
  71. Getzin, S.; Dean, C.; He, F.; A Trofymow, J.; Wiegand, K.; Wiegand, T. Spatial patterns and competition of tree species in a Douglas-fir chronosequence on Vancouver Island. Ecography 2006, 29, 671–682. [Google Scholar] [CrossRef]
  72. Abeli, T.; Cauzzi, P.; Rossi, G.; Adorni, M.; Vagge, I.; Parolo, G.; Orsenigo, S. Restoring population structure and dynamics in translocated species: Learning from wild populations. Plant Ecol. 2016, 217, 183–192. [Google Scholar] [CrossRef]
  73. Mori, A.; Takeda, H. Effects of undisturbed canopy structure on population structure and species coexistence in an old-growth subalpine forest in central Japan. For. Ecol. Manag. 2004, 200, 89–100. [Google Scholar] [CrossRef]
  74. Vidal, S. Plant biodiversity and vegetation structure in traditional cocoa forest gardens in southern Cameroon under different management. Biodivers. Conserv. 2008, 17, 1821–1835. [Google Scholar] [CrossRef]
  75. Wadt, L.H.O.; Kainer, K.A.; Gomes-Silva, D.A.P. Population structure and nut yield of a Bertholletia excelsa stand in Southwestern Amazonia. For. Ecol. Manag. 2005, 211, 371–384. [Google Scholar] [CrossRef]
  76. Hagiwara, Y.; Kachi, N.; Suzuki, J. Combined effects between temporal heterogeneity of water supply, nutrient level, and population density on biomass of four broadly distributed herbaceous species. J. Plant Res. 2012, 125, 77–83. [Google Scholar] [CrossRef] [PubMed]
  77. Craine, J.M.; Dybzinski, R. Mechanisms of plant competition for nutrients, water and light. Funct. Ecol. 2013, 27, 833–840. [Google Scholar] [CrossRef]
  78. Kraaij, T.; Ward, D. Effects of rain, nitrogen, fire and grazing on tree recruitment and early survival in bush-encroached savanna, South Africa. Plant Ecol. 2006, 186, 235–246. [Google Scholar] [CrossRef]
  79. Gómez, D.A.; Aranibar, J.N.; Tabeni, S.; Villagra, P.E.; Garibotti, I.A.; Atencio, A. Biological soil crust recovery after long-term grazing exclusion in the Monte Desert (Argentina). Changes in coverage, spatial distribution, and soil nitrogen. Acta Oecol. 2012, 38, 33–40. [Google Scholar] [CrossRef]
  80. Read, C.F.; Duncan, D.H.; Vesk, P.A.; Elith, J. Surprisingly fast recovery of biological soil crusts following livestock removal in southern Australia. J. Veg. Sci. 2011, 22, 905–916. [Google Scholar] [CrossRef]
  81. Oba, G.; Vetaas, O.R.; Stenseth, N.C. Relationships between biomass and plant species richness in arid-zone grazing lands. J. Appl. Ecol. 2001, 38, 836–845. [Google Scholar] [CrossRef]
  82. King, D.A.; Maindonald, J.H. Tree architecture in relation to leaf dimensions and tree stature in temperate and tropical rain forests. J. Ecol. 1999, 87, 1012–1024. [Google Scholar] [CrossRef]
  83. Kitajima, K.; Mulkey, S.S.; Wright, S.J. Variation in crown light utilization characteristics among tropical canopy trees. Ann. Bot. 2004, 95, 535–547. [Google Scholar] [CrossRef] [PubMed]
  84. Grellier, S.; Barot, S.; Janeau, J.L.; Ward, D. Grass competition is more important than seed ingestion by livestock for Acacia recruitment in South Africa. Plant Ecol. 2012, 213, 899–908. [Google Scholar] [CrossRef]
  85. Kimball, S.; Schiffman, P.M. Differing effects of cattle grazing on native and alien plants. Conserv. Biol. 2003, 17, 1681–1693. [Google Scholar] [CrossRef]
  86. Lin, Y.C.; Chang, L.W.; Yang, K.C.; Wang, H.H.; Sun, I.F. Point patterns of tree distribution determined by habitat heterogeneity and dispersal limitation. Oecologia 2011, 165, 175–184. [Google Scholar] [CrossRef] [PubMed]
  87. Coop, J.D.; Massatti, R.T.; Schoettle, A.W. Subalpine vegetation pattern three decades after stand-replacing fire: Effects of landscape context and topography on plant community composition, tree regeneration, and diversity. J. Veg. Sci. 2010, 21, 472. [Google Scholar] [CrossRef]
  88. Lin, Y.; Hong, M.; Han, G.; Zhao, M.; Bai, Y.; Chang, S.X. Grazing intensity affected spatial patterns of vegetation and soil fertility in a desert steppe. Agric. Ecosyst. Environ. 2010, 138, 282–292. [Google Scholar] [CrossRef]
  89. Alados, C.L.; ElAich, A.; Papanastasis, V.P.; Ozbek, H.; Navarro, T.; Freitas, H.; Cabezudo, B. Change in plant spatial patterns and diversity along the successional gradient of Mediterranean grazing ecosystems. Ecol. Model. 2004, 180, 523–535. [Google Scholar] [CrossRef]
  90. Wiegand, T.; Gunatilleke, S.; Gunatilleke, N. Species associations in a heterogeneous Sri Lankan dipterocarp forest. Am. Nat. 2007, 170, 77–95. [Google Scholar] [CrossRef] [PubMed]
  91. Gutterman, Y. Regeneration of Plants in Arid Ecosystems Resulting from Patch Disturbance; Springer Science & Business Media: Berlin, Germany, 2001. [Google Scholar]
  92. Adler, P.; Raff, D.; Lauenroth, W. The effect of grazing on the spatial heterogeneity of vegetation. Oecologia 2001, 128, 465–479. [Google Scholar] [CrossRef] [PubMed]
  93. Akiyama, T.; Kawamura, K. Grassland degradation in China: Methods of monitoring, management and restoration. Grassl. Sci. 2007, 53, 1–17. [Google Scholar] [CrossRef]
  94. Nedessa, B.; Ali, J.; Nyborg, I. Exploring ecological and socio-economic issues for the improvement of area enclosure management: A case study from Ethiopia. Drylands Coord. Group Rep. 2005, 12, 95–107. [Google Scholar]
  95. Zuo, X.; Zhao, H.; Zhao, X.; Zhang, T.; Guo, Y.; Wang, S.; Drake, S. Spatial pattern and heterogeneity of soil properties in sand dunes under grazing and restoration in Horqin Sandy Land, Northern China. Soil Tillage Res. 2008, 99, 202–212. [Google Scholar] [CrossRef]
  96. He, F.; Duncan, R.P. Density-dependent effects on tree survival in an old-growth Douglas fir forest. J. Ecol. 2000, 88, 676–688. [Google Scholar] [CrossRef]
  97. Getzin, S.; Wiegand, T.; Wiegand, K.; He, F. Heterogeneity influences spatial patterns and demographics in forest stands. J. Ecol. 2008, 96, 807–820. [Google Scholar] [CrossRef]
  98. Coomes, D.A.; Rees, M.; Turnbull, L. Identifying aggregation and association in fully mapped spatial data. Ecology 1999, 80, 554–565. [Google Scholar] [CrossRef]
  99. Defalco, L.A.; Esque, T.C.; Kane, J.M.; Nicklas, M.B. Seed banks in a degraded desert shrubland: Influence of soil surface condition and harvester ant activity on seed abundance. J. Arid Environ. 2009, 73, 885–893. [Google Scholar] [CrossRef]
  100. Giladi, I.; Segoli, M.; Ungar, E.D. Shrubs and herbaceous seed flow in a semi-arid landscape: Dual functioning of shrubs as trap and barrier. J. Ecol. 2013, 101, 97–106. [Google Scholar] [CrossRef]
  101. Aguiar, M.R.; Soriano, A.; Sala, O.E. Competition and facilitation in the recruitment of seedlings in Patagonian steppe. Funct. Ecol. 1992, 66–70. [Google Scholar] [CrossRef]
  102. Hao, H.M.; Lu, R.; Liu, Y.; Fang, N.F.; Wu, G.L.; Shi, Z.H. Effects of shrub patch size succession on plant diversity and soil water content in the water-wind erosion crisscross region on the Loess Plateau. Catena 2016, 144, 177–183. [Google Scholar] [CrossRef]
  103. Brantley, S.T.; Young, D.R. Contribution of sunflecks is minimal in expanding shrub thickets compared to temperate forest. Ecology 2009, 90, 1021–1029. [Google Scholar] [CrossRef] [PubMed]
  104. Aarssen, L.W.; Schamp, B.S.; Wight, S. Big plants—Do they affect neighbourhood species richness and composition in herbaceous vegetation? Acta Oecol. 2014, 55, 36–42. [Google Scholar] [CrossRef]
  105. Pérez, F.L. Geoecological alteration of surface soils by the Hawaiiansilversword (Argyroxiphium sandwicense DC.) in Haleakala’s crater, Maui. Plant Ecol. 2001, 157, 215–233. [Google Scholar] [CrossRef]
  106. Shumway, S.W. Facilitative effects of a sand dune shrub on species growing beneath the shrub canopy. Oecologia 2000, 124, 138–148. [Google Scholar] [CrossRef] [PubMed]
  107. Zheng, J.; He, M.; Li, X.; Chen, Y.; Liu, L. Effects of Salsola passerina shrub patches on the microscale heterogeneity of soil in a montane grassland, China. J. Arid Environ. 2008, 72, 150–161. [Google Scholar] [CrossRef]
Ijerph 15 00946 g001 550
Figure 1. Maps of seeding, adult, aging, and dead A. ordosica groups in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots. Symbol sizes represent the crown of individuals.
Figure 1. Maps of seeding, adult, aging, and dead A. ordosica groups in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots. Symbol sizes represent the crown of individuals.
Ijerph 15 00946 g001
Ijerph 15 00946 g002 550
Figure 2. The crown and height allometric model of A. ordosica population in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots.
Figure 2. The crown and height allometric model of A. ordosica population in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots.
Ijerph 15 00946 g002
Ijerph 15 00946 g003 550
Figure 3. Crown distribution structure of alive, aging, and dead A. ordosica in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Figure 3. Crown distribution structure of alive, aging, and dead A. ordosica in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Ijerph 15 00946 g003
Ijerph 15 00946 g004 550
Figure 4. Spatial autocorrelation of entire A. ordosica population on F5, F10, F15, F25 and CG plots. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Figure 4. Spatial autocorrelation of entire A. ordosica population on F5, F10, F15, F25 and CG plots. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Ijerph 15 00946 g004
Ijerph 15 00946 g005 550
Figure 5. Univariate spatial autocorrelation of four age classes of the A. ordosica population on the F5, F10, F15, F25, and CG plots. S: seedling group, AD: adult group, AG: aging group, and D: dead group. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Figure 5. Univariate spatial autocorrelation of four age classes of the A. ordosica population on the F5, F10, F15, F25, and CG plots. S: seedling group, AD: adult group, AG: aging group, and D: dead group. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Ijerph 15 00946 g005
Ijerph 15 00946 g006 550
Figure 6. Intraspecific spatial correlations among different age classes of the A. ordosica population on F5, F10, F15, F25, and CG sites. S: seedling group, AD: adult group, AG: aging group, D: dead group. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Figure 6. Intraspecific spatial correlations among different age classes of the A. ordosica population on F5, F10, F15, F25, and CG sites. S: seedling group, AD: adult group, AG: aging group, D: dead group. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Ijerph 15 00946 g006
Ijerph 15 00946 g007 550
Figure 7. The Shannon-Wiener index (H), Richness index (R), and Evenness index (E) of communities in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots.
Figure 7. The Shannon-Wiener index (H), Richness index (R), and Evenness index (E) of communities in F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed) plots.
Ijerph 15 00946 g007
Table
Table 1. Basic characteristics information in different plots. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Table 1. Basic characteristics information in different plots. F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Plot NumberCoverage (%)SlopeAltitudeSoil TypeGeographic Coordinates
Latitude (N)Longitude (E)
F5371–3°1393 mSierozem37°50′54.4”107°23′12.7”
F10751–2°1393 mSierozem37°50′49.9”107°23′49.0”
F15401–3°1394 mSierozem37°50′46.3″107°23′48.3″
F26481–3°1396 mSierozem37°50′46.1″107°24′07.8″
CG252–4°1395 mSierozem37°50′40.5″107°24′31.3″
Table
Table 2. Category and basic parameters of A. ordosica population in five plots, F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Table 2. Category and basic parameters of A. ordosica population in five plots, F5 (enclosed 5 years), F10 (enclosed 10 years), F15 (enclosed 15 years), F25 (enclosed 25 years), and CG (continuously grazed).
Plot NumberPopulation Individuals of Different GroupsPopulation Coverage (%)Crown (m)Height (m)
Seedling (%)Adult (%)Aging (%)Dead (%)Total
F5279 (19.1)1089 (74.7)46 (3.2)44 (3.0)145819.70.601 ± 0.262 b0.289 ± 0.102 a
F10257 (5.5)3523 (75.9)536 (11.5)327 (7.0)464371.70.661 ± 0.236 c0.355 ± 0.115 c
F15387 (13.3)2201 (75.9)229 (7.9)83 (2.9)290032.60.556 ± 0.221 a0.341 ± 0.111 b
F25216 (7.5)1294 (45.0)980 (34.1)385 (13.4)287536.70.596 ± 0.226 b0.352 ± 0.126 c
CG53 (6.8)678 (87.0)40 (5.1)8 (1.0)77918.30.814 ± 0.294 d0.387 ± 0.138 d
Data are presented as the mean ±SD. Different small letters in same column indicated significant differences (p < 0.05).

Share and Cite

MDPI and ACS Style

Liu, J.; Zhang, K. Spatial Pattern and Population Structure of Artemisia ordosica Shrub in a Desert Grassland under Enclosure, Northwest China. Int. J. Environ. Res. Public Health 2018, 15, 946. https://doi.org/10.3390/ijerph15050946

AMA Style

Liu J, Zhang K. Spatial Pattern and Population Structure of Artemisia ordosica Shrub in a Desert Grassland under Enclosure, Northwest China. International Journal of Environmental Research and Public Health. 2018; 15(5):946. https://doi.org/10.3390/ijerph15050946

Chicago/Turabian Style

Liu, Jiankang, and Kebin Zhang. 2018. "Spatial Pattern and Population Structure of Artemisia ordosica Shrub in a Desert Grassland under Enclosure, Northwest China" International Journal of Environmental Research and Public Health 15, no. 5: 946. https://doi.org/10.3390/ijerph15050946

APA Style

Liu, J., & Zhang, K. (2018). Spatial Pattern and Population Structure of Artemisia ordosica Shrub in a Desert Grassland under Enclosure, Northwest China. International Journal of Environmental Research and Public Health, 15(5), 946. https://doi.org/10.3390/ijerph15050946

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop