1 Introduction
Grain for green (GFG) refers to projects that promote the gradual cessation of farming on cultivated land with serious problems (e.g., soil erosion, desertification, and high salinity) and the introduction of plants that are suitable for the local environment to improve its ecology (
Li et al., 2020). Previous studies on GFG have mainly focused on vegetation, soil properties, and community patterns of soil bacteria, fungi, and invertebrates (
Ouyang, 2010;
Yang et al., 2010;
Sabais et al., 2011;
Chang et al., 2012; Piché and Kelting, 2015;
Feng et al., 2019;
Hao et al., 2020). Protists are valuable tools for biomonitoring studies in terrestrial ecosystems. They are sensitive to farming systems and contribute significantly to mineralization processes and key ecosystem functions in soil (
Heger et al., 2012;
Domonell et al., 2013). Changes in soil ciliate community structure have also been suggested as bioindicators of environmental stress, improvements in soil quality, CO
2 flux increments in soil, and peat-bog restoration (
Lüftenegger et al.,1985;
Yeates et al., 1991;
Foissner, 1999; Mayzlish and Steiberger, 2004;
Mieczan et al., 2017;
Gabilondo et al., 2018;
Abraham et al., 2019). Several previous studies on the diversity of soil ciliate communities at GFG sites in north-western China have been carried out by
Ning et al. (2011b,
2016,
2017a,
2018,
2019).
Because of their abundance, bacteria, fungi, and protists play important roles in nutrient cycling and energy flow in soils (
Zhou et al., 2012;
Schulz-Bohm et al., 2016;
Xiong et al., 2018; Fiore-Donno et al., 2019;
Liu et al., 2019;
Shen, 1999;
Song et al., 1999;
Chen et al., 2020;
Liu et al., 2020). Traditionally studies of soil ciliates focused mainly on their morphology, assemblage, ecology, and use in toxicology, but in recent years increasing emphasis is placed on molecular phylogeny, relationships with other microorganisms, and the effects of environmental pollution on soil ciliate community structure. (
Liu et al., 2008; Buonanno and Ortenzi, 2010;
Ning et al., 2011a;
Shao et al., 2014a, 2014b;
Geisen et al., 2016;
Lin et al., 2017;
Luo et al., 2017;
Kosiba et al., 2019;
Wanner et al., 2020).
In this study, spatial variations in community patterns of soil ciliates and their potential utility as a bioindicator of environmental restoration resulting from a GFG project at a site in north-western China were studied during 1 year. The main objectives of the study were: (1) to document the community patterns of the soil ciliates; (2) to reveal differences in ciliate community patterns along with the ecological restoration quality of GFG; and (3) to summarize the utilities of the ecological features based on the community structure of the soil ciliates for assessing the ecological restoration in GFG.
2 Materials and methods
2.1 Sampling sites
The four sampling sites were located in Gulang County, north-western China (Fig. 1), east of the Hexi Corridor (37°09′–37°54′ N, 102°38′–103°54′ E). The climate was temperate continental arid with an average annual precipitation of 300 mm and an annual average temperature of 4.9°C. The soil type was Castanozems. The vegetation in the sample area was dominated by Stipa krylovii, Populus davidiana, Picea asperata, and Rosa xanthina. The vegetation coverage and other characteristic parameters are shown in Table 1.
Fig.1 Four sampling sites were investigated in Gulang County, northwestern China. A, site A: GFG for 15 years; B, site B: GFG for 13 years; both A and B are in mixed forests planted with Armeniaca sibirica and Hippophae rhamnoides; C, site C: layland; D, site D: cultivated land. |
Full size|PPT slide
Tab.1 Parameters of plant community at sites A and B. |
Parameter of plant community | A | B |
Species number | 33.000 | 22.000 |
Density | 98.240 | 69.280 |
Coverage | 53.000 | 46.000 |
Shannon-wiener index | 2.278 | 1.984 |
Pielou index | 0.652 | 0.642 |
Simpson index | 0.159 | 0.186 |
Margalef index | 4.111 | 2.816 |
2.2 Experiment design and data collection
Four sampling sites were selected. Sites A and B had 15 and 13 years of GFG, respectively, and dominant plants in the two sites were
Armeniaca sibirica and
Hippophae rhamnoides. Site C was an area of layland transitioning to becoming forested after being formerly cultivated, where the dominant plants were
Capsella bursa-pastoris and
Potentilla anserina. Site D was in an area of cultivated land planted with
Brassica napus. We placed a 10 m
´10 m quadrat at sites A and B to investigate the types and numbers of trees and shrubs and a 1 m × 1 m quadrat to investigate the number of herbaceous plants. According to Flora Reipublicae Popularis Sinicae (1993) and Flora Tsinlingensis (1976), all plant species were identified. After removing branches and fallen leaves from the soil surface, we used the “plum blossom five-point sampling method” for collecting soil samples from the upper 5 cm layer. Additional soil samples were collected using a cutting ring with a volume of 100 cm
3 for determining soil bulk density and porosity. Other data recorded were the geographical coordinates of the sampling points, soil temperature, altitude, vegetation type, and soil texture. A part of the soil sample was brought back to the laboratory to determine soil physical-chemical variables (pH, water content, and electrical conductivity). The remaining part was naturally air-dried indoors for subsequent qualitative and quantitative experiments. During the drying process, the soil samples were covered with straw paper to prevent contamination by air-born ciliate cysts (
Ning et al., 2016).
The non-flooded Petri dish method was used for the qualitative study of ciliates (
Foissner, 1992). Soil (50 g) was placed in a Petri dish with a diameter of 12 cm, and soil extract was added to make the soil samples thoroughly wet but not submerged. Two or three rice grains were added to enrich the growth of bacteria as food for the ciliates. The Petri dishes were incubated at 25°C. Microscopic examination of the ciliates began after 24 h of cultivation. Ciliates were observed
in vivo and by protargol staining. Species were identified using published guides and keys (
Shen et al., 1990;
Foissner, 1993,
1998; Ning and Shen, 1998a, 1998b;
Berger, 1999,
2011;
Adl et al., 2019).
According to
Ning et al. (2019), the directing-culture counting method was used for the quantitative study of ciliates. Briefly, 30 g soil samples were added to a Petri dish with a diameter of 5 cm. Soil extract (30 mL) was added, and the level of the liquid in the Petri dish was marked. The Petri dish was placed in a 25°C incubator, and soil extract was added every day to maintain the liquid level in the dish. All the soil suspension was removed from the Petri dish after 5, 7, 9 and 11 days and the volume of liquid was determined with a measuring tube. Five drops of suspension were placed on a glass slide using a pipette and observed under a microscope. The ciliates in each drop were enumerated, and the number of individuals per mL of suspension (based on 22 drops having a volume of about 1 mL). From these data, the total number of ciliate cells in 30 g soil was calculated.
Soil pH was measured with a potentiometric method, used extraction method to determining soil electrical conductivity. Soil temperature was measured directly using a curved tube thermometer. Water content was determined following the method of
Davidson et al. (1998). Soil porosity was determined by the soil cutting ring method (Vogel and Roth, 2001). Soil organic matter content was determined by the potassium dichromate titrimetric method. Soil total nitrogen, total potassium, and total phosphorus were determined by the Kjeldahl method, flame spectrophotometry, and sodium hydroxide alkali fusion-molybdenum-antimony colorimetric method, respectively (
Bremner., 1960;
Jackson., 1973;
Xue et al., 2004).
2.3 Data analysis
The number of species at the order level was defined as dominant species, and the number of secondary species was defined as inferior dominant species. Species that were observed only once were defined as endemic species; ubiquitous species are defined as those present at each sampling site (Ning and Shen, 1998a, 1998b). One-way ANOVA and the least significant difference (LSD) were used to analyze the correlations between physical-chemical variables and ciliates. The characteristic community parameters of soil ciliates and vegetation were analyzed by cluster and multidimensional scale (MDS) analyses. Data was analyzed by SPSS 19.0, Primer 6.0 (for cluster), MDS and PCA. Histograms was made by Origin 8.0. Arcview Gis 10.6 was used for sample plotting.
3 Results
3.1 Vegetation community characteristics
A total of 4134 plants representing 42 species were observed at sites A and B. This comprised four species (127 individuals) of woody plants and 38 species (4007 individuals) of herbaceous plants. At site A, 33 plant species were recorded, and the dominant species was Poa sphondylodes. At site B, 22 plant species were recorded, and the dominant species was Roegneria kamoji. It can be seen that with the increase of the years of GFG, species number, individual abundance, coverage, Shannon-Wiener index, Pielou index, and Margalef index of vegetation community were all higher at site A than at site B, whereas the Simpson index was higher at site B than at site A (Table 1). That is, the vegetation structure of site A was more complex than site B.
3.2 Soil physical-chemical variables
The soil physical-chemical variables of each sample are shown in Fig. 2. Values for soil porosity, organic matter, water content, and total nitrogen were all higher in the soils at sites A and B than at sites C and D (P<0.05), and were in the following order: A>B>C>D. Values of pH, conductivity, total phosphorus, and total nitrogen were significantly higher in soil from cultivated land (site D) than at sites where GFG had been implemented for 13 years or more (sites A and B) (P<0.05), and were in the following order: D>C>B>A. Soil temperature did not differ significantly among the four sites (P>0.05).
Fig.2 The sol physical-chemical variables (a: Temperature; b: Water content; c: pH value; d: Conductivity; e: Soil porosity; f: Organic matter; g: Total nitrogen; h: Total phosphorus; i: Total potassium) in each season (spring, summer, autumn, winter) and each sampling site (A, B, C, D). |
Full size|PPT slide
3.3 Composition of soil ciliate community
The dominant group was Sporadotrichida, with 35 species, accounting for 30.7% of the soil ciliate species richness. The secondary dominant group was Haptorida, with 21 species accounting for 18.42% of soil ciliate species richness (Table 2).
Tab.2 Community composition of soil ciliates in various sampling sites. |
Phylum | Class | Order | Family | Genus | Species |
| Heterotrichea | Heterotrichida | 1 | 1 | 1 |
Spirotrichea | Euplotida | 1 | 1 | 1 |
| Stichotrichida | 3 | 4 | 6 |
| Sporadotrichida | 2 | 9 | 35 |
| | Urostylida | 1 | 5 | 14 |
| Armophorea | Armophorida | 1 | 1 | 3 |
Ciliophora | Litostomatea | Haptorida | 4 | 6 | 21 |
| Nassophorea | Nassulida | 1 | 1 | 4 |
| | Microthoracida | 1 | 2 | 3 |
| Colpodea | Colpodida | 1 | 1 | 14 |
| | Cyrtolophosidida | 1 | 1 | 1 |
| Prostomatea | Prorodontida | 2 | 2 | 4 |
| Oligohymenophorea | Pleuronematida | 2 | 2 | 4 |
| Phyllopharyngea | Chlamydodontida | 1 | 1 | 3 |
Total | 9 | 14 | 22 | 37 | 114 |
The dominant groups of soil ciliates in various samples underwent obvious succession. Species of Sporadotrichida were dominant at sites A and B, whereas species of Colpodida were dominant at site D (Fig. 3). It can be seen that the number of soil ciliate species increased, the community structure tended to be more complex.
Fig.3 Community composition of soil ciliates at the four sampling sites. |
Full size|PPT slide
3.4 Distribution of soil ciliates
A total of 114 species of soil ciliates were identified at the four sampling sites. At site A, 91 species were identified, representing 33 genera, 20 families, 13 orders, and nine classes, accounting for 79.82% of the total species number. 87 species were identified at site B, representing 31 genera, 19 families, 14 orders, and nine classes, accounting for 76.32% of the total species number. At site C, 74 species were identified, representing 28 genera, 17 families, 11 orders, and eight classes, accounting for 64.91% of the total species number. At site D, 55 species were identified, representing 23 genera, 15 families,10 orders, and eight classes, accounting for 48.25% of the total species number.
The results of the analyses are shown in Fig. 4. In the cluster analysis, sites A and B clustered together, and then clustered with C, whereas site D clustered with sites A, B and C sample at the furthest distance. The results of MDS analysis are consistent with the clustering analysis; that is, sites A, B, and C are close to each other, and site D is the most isolated. Through clustering and MDS analyses, it is shown that sites A and B have a closer similarity with site C than with site D.
Fig.4 Cluster analysis (a) and MDS (b) analysis on species distribution of soil ciliates at the four sampling sites. |
Full size|PPT slide
3.5 Ubiquitous and endemic species
Ubiquitous species are defined as those present at each sampling site. Among the total of 114 ciliate species recorded, 30 were ubiquitous, accounting for 26.32% of total ciliate species richness. These were primarily composed of colpodids, stichotrichids, sporadotrichids, and urostylids. Endemic species were defined as those that appeared only at one sample site. There were 10 endemic species at both sites A and B, accounting for 8.77% of ciliate species. It includes omnivorous ciliates such as sporadotrichids, haptorids, and stichotrichids. There were two endemic species at site C, accounting for 1.75% of the total species, one urostylid and one prorostomatid. There was one endemic species at site D, accounting for 0.88% of the total species. It can be seen that the dominant of ciliate species at the four sites was in the order: D>C>B>A.
3.6 Community parameters of soil ciliates
The species number, individual abundance, and C/P index of soil ciliates at each sampling site are shown in Table 3. In any given season, the number of ciliate species at the four sites was in the order: A>B>C>D, and the individual abundance of ciliates at sites A and B was significantly larger than at sites C and D (P<0.05). The individual abundance of ciliates in any given season was in the order: A>B>C>D, and the total ciliate abundance among the four seasons was in the order autumn>summer>spring>winter. The C/P index of soil ciliates was the reverse of species number and abundance, i.e., the C/P index was in the order D>C>B>A. The number of soil ciliate species and their individual abundance in the soils after GFG (sites A and B) were higher than at sites C and D, whereas the C/P index was the opposite.
Tab.3 Species number, individual abundance and C/P index of soil ciliates in each season (spring, summer, autumn, winter) and at each sampling site (A, B, C, D). |
Sampling sites | Seasons | Species number | Individual abundance | C/P index |
A | Spring | 41 | 72.51±0.26 c | 0.32 |
B | 37 | 57.43±0.26 c | 0.35 |
C | 29 | 48.52±0.28 c | 0.44 |
D | 24 | 38.91±0.36 c | 0.73 |
A | Summer | 46 | 80.28±0.45 b | 0.31 |
B | 33 | 68.37±0.71 b | 0.35 |
C | 32 | 59.26±0.69 b | 0.37 |
D | 21 | 48.56±0.51 b | 0.7 |
A | Autumn | 50 | 110.61±0.84 a | 0.32 |
B | 41 | 96.89±1.06 a | 0.4 |
C | 39 | 77.96±0.54 a | 0.37 |
D | 28 | 60.01±0.57 a | 0.69 |
A | Winter | 36 | 57.55±0.88 d | 0.43 |
B | 32 | 48.24±0.65 d | 0.45 |
C | 28 | 38.76±0.39 d | 0.5 |
D | 20 | 28.59±0.43 d | 0.8 |
The soil ciliate community indices at each of the four sites are shown in Table 4. For each season, the values of the species diversity and richness indices at sites A and B were significantly higher than those at sites C and D (P<0.05). These indices were in the order: A>B>C>D.
Tab.4 Shannon-Wiener diversity index, Margalef richness index, Pielou evenness index and Simpson dominance index of soil ciliates in each season (spring, summer, autumn, winter) and at each sampling soie (A, B, C, D). |
Sampling sites | Seasons | Shannon-Wiener diversity index | Margalef richness index | Pielou evenness index | Simpson dominance index |
A | Spring | 2.290±0.011 a | 2.849±0.039 a | 0.889±0.005 b | 0.281±0.014 a |
B | 2.052±0.016 b | 2.071±0.024 b | 0.916±0.006 ab | 0.159±0.005 b |
C | 1.807±0.012 c | 1.479±0.023 c | 0.944±0.003 a | 0.196±0.004 a |
D | 1.695±0.007 d | 1.363±0.003 c | 0.946±0.004 a | 0.198±0.007 b |
A | Summer | 2.250±0.018 a | 2.965±0.035 a | 0.853±0.004 c | 0.140±0.005 a |
B | 2.084±0.011 b | 2.178±0.038 b | 0.899±0.004 b | 0.147±0.003 b |
C | 1.879±0.011 c | 1.567±0.023 c | 0.940±0.004 a | 0.168±0.003 b |
D | 1.804±0.017 c | 1.499±0.019 c | 0.942±0.003 a | 0.172±0.010 c |
A | Autumn | 2.673±0.016 a | 4.041±0.026 a | 0.892±0.004 a | 0.093±0.010 a |
B | 2.480±0.012 b | 3.321±0.037 b | 0.891±0.005 a | 0.110±0.005 c |
C | 2.222±0.025 c | 2.617±0.050 c | 0.884±0.006 a | 0.127±0.004 c |
D | 1.918±0.011 d | 1.809±0.047 d | 0.904±0.006 a | 0.170±0.006 c |
A | Winter | 1.908±0.023 a | 1.879±0.070 a | 0.891±0.002 b | 0.178±0.005 a |
B | 1.749±0.018 b | 1.549±0.038 a | 0.901±0.008 b | 0.197±0.004 a |
C | 1.628±0.014 c | 1.259±0.033 b | 0.948±0.004 a | 0.211±0.003 a |
D | 1.464±0.021 d | 1.185±0.032 b | 0.910±0.014 b | 0.242±0.007 a |
3.7 Relationship between soil ciliates and soil physical-chemical variables
As shown in Table 5, for temperature, the most significant negative correlation was with the dominance index (P<0.01). For water content, the most significant positive correlation was with species diversity (P<0.01). For soil porosity, the most significant positive correlations were with species number and diversity (P<0.01). For total nitrogen, the most significant positive correlations were with species number, individual abundance, diversity, and richness (P<0.01), and the most significant negative correlation was with species evenness (P<0.01). For organic matter, the most significant positive correlations were with species number and richness (P<0.01). Therefore, water content, soil porosity, total nitrogen, and organic matter were the key soil physical and chemical variables affecting the soil ciliate communities.
Tab.5 Relationships between the community patterns of soil ciliates and the soil physical-chemical variables (Pearson’s correlation coefficient). |
Physical-chemical variables | Species number | Individual abundance | Diversity index | Evenness index | Dominance index | Richness index |
Temperature (°C) | 0.168 | 0.472 | 0.379 | -0.235 | -0.631** | 0.340 |
Water content (%) | 0.584* | 0.558* | 0.651** | -0.264 | -0.267 | 0.560* |
Soil porosity (%) | 0.628** | 0.608* | 0.697** | -0.335 | -0.310 | 0.605* |
T-nitrogen (g Kg–1) | 0.807** | 0.671** | 0.676** | -0.670** | -0.261 | 0.731** |
T-potassium (g Kg–1) | -0.205 | 0.102 | 0.011 | 0.428 | -0.369 | 0.018 |
T-phosphorus (g Kg–1) | -0.524* | -0.336 | -0.331 | 0.246 | 0.063 | -0.315 |
Organic matter (g Kg–1) | 0.770** | 0.564* | 0.621* | -0.567* | -0.033 | 0.665** |
pH value | -0.265 | -0.082 | -0.008 | 0.418 | 0.219 | -0.056 |
Electric conductivity (μS) | -0.604* | -0.394 | -0.412 | 0.343 | 0.122 | -0.376 |
The results of principal component analysis of soil ciliate community parameters and soil physical-chemical variables are shown in Fig. 5. Nine physical-chemical variables were measured: temperature (T), water content (WC), soil porosity (SP), total nitrogen (TN), total potassium (TK), total phosphorus (TP), organic matter (OM), electrical conductivity (SC), and pH. Six ciliate community parameters were included in the analysis: species number (SN), individual abundance (D), diversity index (H′), evenness index (E), dominance index (C) and richness index (M). As can be seen from the figure, the line segments of vectors for WC, SP, OM, TN, SN, D, H' and M were relatively long and the angles between them were relatively small, indicating significant correlations among these factors. The angles between the vectors of other factors were greater than 90 degrees, indicating no significant correlations among them. Therefore, water content, soil porosity, organic matter, and total nitrogen were the main soil physical-chemical variables that affected the soil ciliate communities, which was consistent with Pearson correlation analysis.
Fig.5 The PCA analysis between the community characteristic parameters of soil ciliates and soil environmental variables. |
Full size|PPT slide
3.8 Relationship between soil ciliates and vegetation parameters
The results of cluster and MDS analysis of the community patterns of soil ciliates and vegetation community characteristic parameters at the four sampling sites are shown in Fig. 6. In the cluster analysis, the species number and individual abundance of soil ciliates were highly correlated with vegetation coverage, species number of vegetation, and individual abundance of vegetation. Soil ciliate species richness, diversity, and evenness had high correlations with vegetation richness, diversity, and evenness. Soil ciliate dominance index and vegetation community dominance index had poor correlations with other indices. The results of the MDS analysis were consistent with the clustering analysis.
Fig.6 The cluster (a) and MDS (b) analysis on soil ciliates and plant community at the four sampling sites. |
Full size|PPT slide
4 Discussion
4.1 Effects of different patterns of GFG on vegetation communities
GFG is one of the most important measures to restore vegetation and improve ecology. The indexes used by domestic scholars to study vegetation diversity include those of species richness, diversity, evenness, and dominance (
Wan et al., 2014). It is essential to analyze plant community species diversity and reveal community succession and stability (Kai and Jaanus, 1999; Burton and Mark, 2001). In this study, the measures of coverage, number of species and individual abundance, diversity, richness, and evenness of the plant community were higher for site A than for site B, while the dominance index was the opposite. These findings reflect the greater number of years of GFG at site A (15 years) than at site B (13 years). The effects of GFG on plant species abundance and community composition are mainly caused by enhanced deciduous vegetation layers, root systems and nitrogen fixation. With increasing stocks of leaf litter in forests, the storage of C and N in soil increases, and soil’s physical properties are improved (
Zhang et al., 2011;
Shang et al., 2015;
Liu et al., 2020).
4.2 Effects of GFG on soil physical and chemical variables under different return modes
With the conversion of habitats by GFG, the community structure of above-ground vegetation tends to be stable, and the soil condition is improved. In this study, soil physical and chemical variables, such as soil water content, soil porosity, organic matter, and total nitrogen, were significantly higher where GFG has been implemented (sites A and B) compared to sites C and site D (
P<0.05). The soil water content of site A was 70.92% higher than site D, mainly because with the progress of GFG and the increase of aboveground vegetation, there were increases in the surface litter, humus, and soil water storage capacity, and a decrease in transpiration (
Zhang et al., 2018;
Mo et al., 2020). Furthermore, the soil porosity of site A was 4.09% higher than site D, at the sites where GFG had been implemented, probably due to the increased development of plant root systems. The results showed that the soil porosity of the GFG sites was significantly higher than that of the cultivated land, which is consistent with the present findings (
Ning et al., 2019). Soil organic matter of site A was 62.41% higher than site D, probably as a result of higher microbial biomass, which decomposes vegetation litter and humus on the soil surface resulting in more soil organic matter (
Cui et al., 2010;
Lin et al., 2011;
Zhang et al., 2018). The total nitrogen content of site A was 58.59% higher than site D, probably as a result of nitrogen from the leaf litter entering the soil with rainwater following decomposition by microorganisms. The total nitrogen content in soil was consistent with that predicted by the vertical variation law of soil organic matter content (Pang and Bao, 2011).
4.3 Effects on soil ciliate community of GFG under different return regimes
The diversity of soil animal communities changes with different vegetation types; that is, the diversity of surface vegetation communities affects the diversity of soil animal species (
Chauvat et al., 2011). Furthermore, the population dynamics of soil animals respond to the changes in vegetation and environmental restoration, so it is susceptible to external environmental factors. Therefore, the dynamic changes of soil animal communities can serve as an important bioindicator of vegetation restoration and environmental improvement (
Sabais et al., 2011;
Zhang et al., 2012). The findings of the present study showed that with the progress of GFG, the species number, individual abundance, diversity index, richness index, and evenness index of soil ciliate communities at each site were in the order A>B>C>D, while the dominance index and C/P index showed the opposite. Due to its longer return time, site A had the most complete vegetation community and the most stable and complex ciliate community, indicating that the soil conditions were the most suitable for the survival of soil ciliates. The C/P index at site D was higher than that at sites A and B, indicating that the soil environment of farmland is not as suitable for the survival of soil ciliates, which is consistent with previous findings (
Ning et al., 2017b). Correlation analysis showed that water content, soil porosity, total nitrogen, organic matter, and soil ciliate community parameters were the most significantly positively correlated, also consistent with the findings of
Zou et al. (2009). Colpodids were the dominant ciliate species at site D (cultivated land), whereas sporadtrichids were dominant in the soil at sites A and B (where GFG had been implemented).
5 Conclusions
This study shows that with the implementation of GFG, the above-ground vegetation community structure is complicated. When the physical and chemical variables of the surface soil are all significantly changed, the species number, individual abundance, diversity index, evenness index, richness index of soil ciliates (increases), and the dominance index (decreases) change. However, the community structure of soil ciliates remains complicated. The correlation analysis showed that water content, soil porosity, total nitrogen, and organic matter had the most influence on the characteristics of the soil ciliate community. Therefore, the community parameters of soil ciliates can be used as an evaluation index of the ecological restoration effect of GFG.
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}