- 1School of Water and Environment, Chang’an University, Xi’an, China
- 2Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang’an University, Xi’an, China
- 3Key Laboratory of Eco-hydrology and Water Security in Arid and Semi-arid Region of Ministry of Water Resources, Chang’an University, Xi’an, China
- 4Key Laboratory of State Forest Administration on Soil and Water Conservation & Ecological Restoration of Loess Plateau, Shaanxi Academy of Forestry, Xi’an, China
Introduction: The impact of groundwater table depth (GTD) on bacterial communities and soil nutrition in revegetated areas remains unclear.
Methods: We investigated the impacts of plant growth and soil physicochemical factors on rhizosphere bacterial communities under different GTD.
Results: The four plant growth indices (Pielou, Margalef, Simpson, and Shannon-Wiener indices) and soil water content (SWC) at the Artem and Salix sites all showed a decreasing trend with increasing GTD. Salix had a higher nutrient content than Artem. The response of plant rhizosphere bacterial communities to GTD changes were as follows. Rhizosphere bacteria at the Artem and Salix sites exhibited higher relative abundance and alpha diversity in SW (GTD < 5 m) compared than in DW (GTD > 5 m). Functional microbial predictions indicated that the rhizosphere bacterial communities of Artem and Salix promoted carbon metabolism in the SW. In contrast, Artem facilitated nitrogen cycling, whereas Salix enhanced both nitrogen cycling and phototrophic metabolism in the DW.
Discussion: Mantel test analysis revealed that in the SW of Artem sites, SWC primarily governed the diversity of rhizosphere and functional bacteria involved in the nitrogen cycle by affecting plant growth. In DW, functional bacteria increase soil organic carbon (SOC) to meet nutrient demands. However, higher carbon and nitrogen availability in the rhizosphere soil was observed in the SW of the Salix sites, whereas in DW, carbon nutrient availability correlated with keystone bacteria, and changes in nitrogen content could be attributed to nitrogen mineralization. This indicates that fluctuations in the groundwater table play a role in regulating microbes and the distribution of soil carbon and nitrogen nutrients in arid environments.
1 Introduction
In recent years, irrational human activities, including water diversion for irrigation and overexploitation, coupled with the impacts of climate change, have resulted in significant groundwater depletion. This phenomenon poses a severe threat to vegetation survival in arid and semi-arid regions worldwide (Castle et al., 2014; Barbeta et al., 2015). Groundwater is a crucial water resource in these drylands and is tightly linked to vegetation productivity and the ecosystem carrying capacity (Karthe, 2018; Ma et al., 2022). Alterations in the groundwater table, whether upward or downward, exert a substantial influence on the soil water content, consequently affecting regional vegetation growth (Chen et al., 2016). Vegetation growth predominantly relies on water within the capillary zone above the groundwater surface (Barron-Gafford et al., 2021). Numerous studies have demonstrated that drought disrupts vital ecosystem processes, pushing the sustainable application of water resources close to vegetation growth limitations, particularly in revegetated areas (Han et al., 2020; Li et al., 2021). Although vegetation restoration plays a pivotal role in rehabilitating degraded ecosystems, unregulated large-scale expansion of vegetation can disrupt the overall water balance in a region, potentially leading to soil fertility loss and detrimental alterations to the soil microbial structure. This can affect the progress of ecological restoration (Wall et al., 2004; Schlaepfer et al., 2017; Jia et al., 2020).
Numerous studies have demonstrated that a decrease in the groundwater table induces alterations in plant communities, modifies soil physicochemical attributes, and affects soil nutrient cycling, such as carbon and nitrogen dynamics (Liu et al., 2018; Zeng J. et al., 2020; Shen et al., 2022). Xi et al. (2018) established a robust association between the groundwater table and plant growth, particularly root development and distribution, emphasizing that GTD is a key determinant of rooting depth under relatively stable conditions. Mao et al. (2018) indicated the positive impact of groundwater table on both aboveground and belowground plant traits, along with a concurrent reduction in soil moisture content with decreasing groundwater table in the scope of shallow groundwater table depth. Fluctuations in the groundwater table exert diverse effects on soil physicochemical properties and environmental factors such as oxygen content (O2) (Meyerhof et al., 2016), temperature (Lai et al., 2021), humidity, atmospheric pressure and chemical composition (Shen et al., 2022). Vegetation and inter-root microorganisms are sensitive to environmental changes, allowing for rapid adaptation (Bissett et al., 2010; Hill et al., 2015). Therefore, variations in the groundwater table will inevitably affect the composition of inter-root microbial communities, resulting in distinct assemblages adapted to shifting groundwater conditions.
Microbial communities in the rhizosphere exhibited distinct characteristics among the different vegetation types. To some extent, changes in rhizosphere soil physicochemical properties are influenced by root regeneration and rhizosphere metabolism (Galland et al., 2019). To rehabilitate soil in the Mu Us Sandy Land, shrubs such as Artem and Salix have been extensively introduced. Artemisia ordosica has been reported to release allelochemicals through rain leaching, inhibiting the growth of other plants and potentially establishing a dominant community (Wang et al., 2022). In addition, Salix litter significantly increased soil organic carbon (SOC), total nitrogen (TN), and soil water content (SWC) (Li T. et al., 2019). The involvement of microorganisms is vital for litter decomposition and influences nutrient cycling within the ecosystem (Osburn et al., 2022). Rhizosphere microbes play a critical role in nutrient cycles, particularly in soil carbon and nitrogen dynamics, as well as in the conversion of both organic and inorganic matter, especially rhizosphere functional microorganisms (Zhao et al., 2020; Kunkun et al., 2022). Functional microorganisms, including photoheterotrophic, chemotrophic, and those involved in the degradation of aromatic compounds and methane generation, govern the decomposition and oxidative synthesis of organic carbon under varying conditions (Liang et al., 2019; Huang et al., 2020). Nitrogen fixation, nitrification, and nitrate reduction are the primary microorganisms responsible for regulating the nitrogen cycle (Chen et al., 2020; Xiang et al., 2020). Plant root secretions, including organic acids, affect microbial species and activity (Zhao et al., 2020). Nonetheless, the effects of the application of Artem and Salix in sandy soil restoration on properties such as SOC, TN, SWC, and TP, as well as soil microbial composition, remain unclear. Addressing this knowledge gap requires a systematic assessment of how associations between plant species and soil properties, including microbial diversity, evolve during the revegetation process.
The Mu Us Sandy Land, located in the arid and semi-arid regions of northwestern China, has encountered substantial ecological and environmental challenges to vegetation restoration due to drought and soil impoverishment (Veldkornet et al., 2016). Effective vegetation restoration measures have emerged as a pivotal strategy for conserving soil and water resources (Wu and Ci, 2002). However, certain areas have witnessed deterioration in plantation forest health (Bi et al., 2021). Since the implementation of the vegetation restoration program in 1999, rapid vegetation growth has induced notable fluctuations in the local groundwater cycle, potentially affecting the regional groundwater table (Jia et al., 2020). Furthermore, these alterations in the groundwater table can probably affect vegetation growth and nutrient cycling (Han et al., 2020; Jiakun et al., 2023), posing a serious challenge in the battle against secondary desertification in Mu Us Sandy Land.
This study utilized two common sand-fixing vegetation types, Salix and Artem sites, to examine changes in plant community composition, soil physicochemical factors, and rhizosphere bacterial communities across varying GTD. The hypotheses were as follows. In both the SW and DW sites, (1) it was postulated that significant variations would be evident in the characteristics of plants, soil physicochemical factors, composition of rhizosphere bacterial communities, functional bacteria, and co-occurrence networks of sand-fixing vegetation. (2) an extension of the hypothesis suggested that the impact of plant diversity and soil physicochemical factors on soil bacterial communities and functional bacteria would display noteworthy distinctions between the two sites. (3) it was conjectured that the regulation of soil nutrient levels might be contingent on differing levels of vegetation bacterial diversity and the presence of specific functional bacteria within these ecosystems.
2 Materials and methods
2.1 Study sites
The study area, spanning approximately 259 km2, is situated on the southeastern fringe of the Mu Us Sandy Land (109.13°–109.72°E, 38.32°–38.45°N) at the transitional boundary between the Mu Us Sandy Land and the Loess Plateau. This region has a unique geographical context and hosts an exceptionally sensitive and delicate ecological environment (Figure 1). The average annual temperature ranges from 8.15°C to 10°C, with peak temperature typically occurring in July. The average annual precipitation ranges from 137 mm to 578 mm, gradually decreasing toward the northwest. Precipitation exhibits substantial annual variability, with 60–75% concentrated from July to September. Evaporation is intense, averaging 1,606–3,240 mm annually (Zhao et al., 2021). The total annual radiation ranges from 5,140 MJ/m2 to 6,281 MJ/m2, ensuring ample sunshine availability.
The study area exhibits complex and undulating terrain with a gradual descent from northwest to southeast. The groundwater table displays significant spatial variation, ranging from 0 to 60 m (Rao et al., 2022). In this region, the limit of evaporation depth for the groundwater table generally falls within the range of 4–6 m. For our delineation between SW and DW, we adopted a groundwater depth of 5 m based on distinct differences in soil water content and vegetation coverage, as established in our previous investigation (He et al., 2023). Common plant species in sandy land areas include Artemisia desertorum, Salix cheilophila, Corethrodendron fruticosum, Caragana korshinskii, and Amorpha fruticosa. Additionally, mudflats and river valleys feature aquatic and halophytic vegetation including Phragmites australis, Thermopsis lanceolata, Aster tataricus, Sueda salsa, and Neotrinia splendens.
2.2 Experimental design and sample collection
According to different groundwater depth conditions (1–10 m), 20 sample sites were established to investigate the responses of plant diversity, soil properties, and bacterial communities to GTD changes in the Artem and Salix coverage areas. Groundwater table data were sourced from the Yulin Water Resources Bureau of Shaanxi Province. First, within each sample site, 5 m × 5 m quadrats were set up to assess plant diversity, including parameters such as plant coverage, plant height, and associated plant species. Subsequently, rhizosphere soil samples were obtained from a depth of 0–20 cm below the ground surface using a small sterile shovel. The ‘five point method’ was employed within the same quadrat, preceded by the removal of surface litter. Finally, three soil samples were collected from each quadrat after on-site mixing. All soil samples were collected from the rhizosphere. One section was allocated for the analysis of soil chemical and physical properties, while the remainder was preserved at −80°C to facilitate future soil DNA extraction. Soil sampling was conducted in July 2021, and subsequent testing was conducted within a two-week timeframe from the date of collection.
2.3 Measurement of soil properties
Soil water content (SWC, %) was measured using the drying method (constant temperature of 105°C, drying for 8 h). The soil pH was assessed using a soil extract with a soil–water ratio of 1:2.5. Soil organic carbon content (SOC, g/kg) was determined using the potassium dichromate oxidation external heating method (Nelson and Sommers, 1982). Total nitrogen (TN, g/kg) was quantified using the semi-micro-Kjeldahl method (Bremner and Mulvaney, 1982). Available phosphorus (AP, mg/kg) was determined using the NaHCO3 extraction-molybdenum antimony anti-colorimetric method (Olsen and Sommers, 1982).
2.4 High-throughput sequencing of microbial communities
Total soil DNA was extracted from the soil samples (0.5 g) using the OMEGA Soil DNA Kit (M5635-02) (Omega Bio-Tek, Norcross, GA, USA). The extracted DNA was qualitatively assessed via agar gel electrophoresis and its concentration and purity were determined using a NanoDrop NC2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Subsequently, it was stored at −20°C to prevent evaporation. Amplification of the 16S rRNA gene was performed using primers 338F (5′-GCCAGCMGCCGCGGTAA-3′) and 806R (5′-CCGGACTA CHVGGGTWTCTAAT-3′), following the approach outlined. The PCR conditions consisted of 25 cycles, including pre-denaturation at 98°C for 2 min, denaturation at 98°C for 30 s, annealing at 55°C for 30 s, and elongation at 72°C for 5 min. After purification of the PCR products by gel electrophoresis, dsDNA was accurately quantified using Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). The samples were homogenized based on the quantified values. The integrity of the DNA fragments was assessed using Illumina measurements from Shanghai Paisenor Biotechnology Co. LTD (Shanghai, China). The resulting DNA stock solution was stored at −20°C and diluted tenfold with ddH2O before being stored at 4°C. Each sample was analyzed in four replicates.
2.5 Data analysis
Plant community diversity and evenness indices were obtained from a study by Christian et al. (2022). Variations in the relative abundance of rhizosphere bacteria, alpha diversity, plant diversity, and edaphic variables in SW and DW were assessed using one-way analysis of variance (ANOVA). To investigate the impact of GTD on the composition of rhizospheric bacterial taxa, both alpha and beta diversities were calculated. Four specific indices were employed to comprehensively evaluate the alpha diversity of bacterial communities. Chao1 (Chao, 1984) represented richness, whereas Margalef, Shannon and Simpson indices denoted diversity. Pielou evenness index (Pielou, 1966) was used to characterize the evenness. Beta diversity was determined through Principal Coordinate Analysis (PCoA) based on unweighted UniFrac dissimilarity matrices, exploring differences in soil bacterial community composition between the two sandy plant species under SW and DW conditions.
The network structure was defined using igraph packages incorporating different topological parameters, such as modularity, average degree, and inter-community interactions. For the statistical analysis of rhizosphere soil bacterial community networks in both sandy plants across various GTD, igraph packages were used. Gephi software was used to visualize the co-occurrence networks (Bastian et al., 2009). Rhizosphere bacterial abundance, diversity, and functional bacteria were assessed under diverse GTD conditions through high-throughput 16S rRNA gene sequencing using the Functional Annotation of Prokaryotic Taxa (FAPROTAX) approach. Mantel tests were used to detect the connections between plant diversity, soil edaphic variables, microbial community structure, and functional microorganisms, utilizing the Personalbio Gene Cloud Platform.1
3 Results
3.1 Characteristics of plant communities and edaphic factors with GTD
All four diversity indices for both Artem and Salix displayed a decreasing trend as GTD increased (Figure 2). The Pielou index, reflecting evenness in species distribution, exhibited relatively low values, ranging from 0.04 to 0.26 (Artem) and 0.01 to 0.14 (Salix), indicating that the overall distribution of plant species in the sandy soil was not uniform (Figure 2C). The Margalef index, indicative of species richness, ranged from 0.18 to 0.62 (Artem) and 0.14 to 0.49 (Salix), indicating non-uniform species distribution in both Artem and Salix plots, with varying degrees of alteration (Figure 2A). Simpson’s index ranged from 0.21 to 0.69 (Artem) and 0.05 to 0.74 (Salix), indicating a relatively simplified community or habitat species structure and composition, coupled with high interspecific heterogeneity (Figure 2D). Moreover, the Shannon-Wiener index, which provides a comprehensive measure of species richness and homogeneity, exhibited a broad range, varying from 0.79 to 1.62 (Artem) and 0.6 to 2.12 (Salix) (Figure 2B). This wide variability demonstrates the significant impact of groundwater depth on plant growth in the Mu Us sandy soil.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g002.jpg)
Figure 2. Variation of plant community with increasing groundwater table depth. Panels (A-D) represent Margalef, Shannon-Wiener, Pielou evenness, and Simpson diversity index of Artem and Salix sites, respectively.
With increasing GTD, the SWC exhibited a downward trend at both plant sites. Notably, the SWC of Salix consistently exceeded that of Artem, and displayed an evident declining pattern (Figure 3A). This suggests that the rhizosphere soil of Salix has a higher water availability than that of Artem. Changes in soil nutrients (SOC, TN, and AP) and pH in response to GTD followed quadratic polynomial trends, indicating initial increments followed by subsequent declines. The nutrient content of Salix consistently surpassed that of Artem, suggesting a higher soil nutrient content of Salix (Figures 3B–D). In addition, the soil pH of both plants exhibited significant alterations within the GTD range of 4–6 m (Supplementary Figure S1), indicating the strong influence of GTD within this interval.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g003.jpg)
Figure 3. Changes in soil physicochemical properties of both plants with increasing groundwater table depth. Panels (A-D) represent soil water content (SWC), soil organic carbon (SOC), total nitrogen (TN), and Available phosphorus (AP) of Artem and Salix sites, respectively.
3.2 Characteristics of soil bacterial communities with GTD
Characteristic of bacterial community composition.
The soil bacterial communities exhibited significant differences between the sand-fixing plants under various GTD. Specifically, Artem and Salix ASVs were notably more abundant at a GTD < 5 m than GTD > 5 m (Supplementary Figure S2). Hence, it is verified that the groundwater depth of 5 m can be used as the dividing line between SW and DW. Although the abundance of Artem’s rhizospheric bacteria at the order, family, and genus levels were higher at GTD < 5 m than GTD > 5 m, the difference was not statistically significant. In contrast, the rhizosphere bacterial community of Salix at GTD < 5 m exhibited significantly greater species diversity at the phylum, class, order, and family levels (p < 0.05) than at GTD > 5 m (Supplementary Table S1). In summary, GTD of 5 m influenced both the quantity and species composition of rhizospheric microbial communities for both plants.
The relative abundance of major bacterial phyla did not differ significantly between the rhizosphere bacterial communities of Artem and Salix under varying GTDs (Figure 4). The dominant phyla in the rhizosphere microbial community (with a relative abundance >1%) were Actinobacteria, Proteobacteria, Chloroflexi, and Acidobacteria, collectively constituting 86.49–90.47% of the total relative abundance. Specifically, the relative abundance of Chloroflexi in the rhizosphere bacterial community was significantly higher at DW Salix sites than at Artem sites, whereas the relative abundance of Acidobacteria was notably higher at SW Artem sites than at Salix sites (p < 0.05). Actinobacteria, Thermoleophilia, Alphaproteobacteria, and Gammaproteobacteria were the dominant bacterial classes (with a relative abundance >1%), collectively representing 44.06–49.72% of the total relative abundance (Supplementary Table S2). Significantly divergent trends were also observed in some major bacterial classes within the rhizosphere communities of Artem and Salix under different GTD, including Subgroup_6 (p < 0.05), Chloroflexia (p < 0.05), Deltaproteobacteria (p < 0.05), Anaerolineae (p < 0.05) and NC10 (p < 0.05) (Supplementary Table S3). These findings indicate that variations at the major bacterial classes were more pronounced than those at the phylum level in the rhizosphere soil of both plants in response to changing GTD. Consequently, the relative abundance of sandy soil bacteria was influenced by both the GTD and plant species.
3.2.1 Changes in bacterial community diversity
The alpha diversity (Chao1, Shannon-Wiener, Simpson diversity, and Pielou evenness index) of rhizosphere bacteria in both Artem and Salix were notably higher at GTD < 5 m (Figure 5). However, no significant differences in bacterial alpha diversity were observed between the SW and DW groups. The beta diversity of the soil bacterial community was further analyzed through principal coordinate analysis (PCoA) (Figure 6), complemented by statistical assessment via multiple analysis of variance (PERMANOVA) (Table 1). The first and second axes of the principal coordinate analysis explained 18 and 15.3% of the unconstrained ranking variance for rhizosphere bacterial communities across different GTD conditions in both Artem and Salix, respectively. Figure 6; Table 1 illustrate that using GTD = 5 m as the dividing line, the rhizosphere bacteria of Artem and Salix formed distinct and stable clusters in the SW and DW regions, respectively. Notably, the bacterial communities in the rhizosphere of the Artem (DW) sites were significantly different from those in the other groups (PERMANOVA: F = 2.00, p = 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g005.jpg)
Figure 5. Alpha diversity of soil bacterial communities of two plants across different groundwater areas. Panels (A-D) represent Chao1, Shannon-Wiener, Simpson diversity and Pielou evenness index of rhizosphere bacteria community in the shallow groundwater table depth range (SW) and deep groundwater table range (DW) of Artem and Salix sites, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g006.jpg)
Figure 6. Ordination plot derived from principle coordinate analyses of weighted UniFrac distances of soil rhizosphere bacterial communities of two plants in varying groundwater conditions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-t001.jpg)
Table 1. Differences in soil bacterial communities between the two plants across different groundwater areas (PERMANOVA).
3.2.2 Changes in bacterial community structure
Co-occurrence network analysis was used to investigate the ecological interactions among the microbial species in diverse environmental contexts. In this study, we constructed a co-occurrence network based on the microbial taxa at the ASV level of two plants subjected to varying groundwater depths (Figure 7). The results showed that the soil bacterial community of Artem in SW exhibited a more extensive and stable network structure, whereas the bacterial network of Salix in both SW and DW displayed a similar structure and was largely unaffected by changes in GTD. In the co-occurrence network, Actinobacteria, Proteobacteria, Chloroflexi, and Acidobacteria exhibited robust connectivity with the surrounding nodes, constituting 88.79–92.02% of all nodes (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g007.jpg)
Figure 7. Panels (A-D) exhibit co-occurrence networks of soil rhizosphere bacteria community in the shallow groundwater table depth range (SW) and deep groundwater table range (DW) of Artem and Salix sites, respectively. Connections indicate significant correlations (Screening conditions: Relative abundance of bacterial phyla>0.004). The nodes are colored by phylum and represent an operational taxonomic unit (97% sequence identify threshold, ASV). The size of each node is proportional to the number of connections (degrees). The thickness of each connection between two nodes (edge) is proportional to the values of Spearman’s correlation coefficient.
The topological properties of the bacterial co-occurrence network are presented in Supplementary Table S4, which displays notable discrepancies below the threshold of 0.7. Specifically, at GTD < 5 m, Artem exhibited a greater number of network nodes and edges (230 nodes and 9,326 edges) than GTD > 5 m (132 nodes and 1961 edges). In the case of Salix, the number of network nodes and edges at GTD < 5 m (223 nodes and 5,518 edges) was nearly equivalent to that at GTD > 5 m (238 nodes and 6,095 edges) (Supplementary Table S4). Positive interactions accounted for 52.27 to 55.6%, while negative interactions ranged from 44.38 to 47.73% across different GTD levels. This suggests that the rhizosphere bacterial communities of both sand plants adapted to the arid soil conditions in the Mu Us region by engaging in mutual facilitation, coordination, predation, and antagonism, thereby forming a tight network. In the Artem co-occurrence network at a GTD < 5 m, the total number of links, mean node degree, complexity index, and mean path length were higher, indicating shorter links between nodes, a higher degree of linkage, greater complexity, and a more stable rhizosphere bacterial community. However, Salix exhibited distinct topological characteristics. In both the SW and DW, the total number of connections, average node degree, average path length, and complexity index in the co-occurrence network were similar. This further underscores the robust adaptability of the rhizosphere microorganisms in Salix to fluctuating groundwater table.
Four nodes meeting the criteria (Zi ≥ 2.5, Pi ≤0.62) were selected within the module hubs through an assessment of intra-module and inter-module connectivity within the bacterial community co-occurrence network. These nodes, specifically ASV17731, ASV187664, and ASV233705, were identified within the Artem rhizosphere bacterial co-occurrence network in the DW (Supplementary Figure S3). Additionally, ASV192677 was identified in the Salix rhizosphere bacterial co-occurrence network in the SW. Notably, these four keynote species were from Actinobacteria, Firmicutes, and Acidobacteria. ASV17731 and ASV192677 belong to the Subgroup_17 class of Actinobacteria, representing significant species in this phylum. They are associated with soil fertility enhancement due to their capacity for organic matter decomposition. ASV187664, a member of the phylum Firmicutes, participates in soil respiration metabolism and chemoheterotrophic processes, demonstrating resilience under harsh environmental conditions. ASV233705, which originates from Catenulispora within Actinobacteria, serves as a nutrient provider to plant roots, ensuring optimal plant growth.
3.2.3 Changes in functional bacterial communities
A chord diagram depicting the relative abundances of the top 15 functional bacteria was generated (Figure 8). The functions exhibited by the rhizosphere bacterial communities of Artem and Salix under varying GTD primarily encompassed chemoheterotrophy (33.55–37.77%), aerobic chemoheterotrophy (33.12–37.24%), aromatic compound degradation (3.15–5.31%), hydrocarbon degradation (2.84–4.93%), nitrate reduction (2.22–3.95%), nitrification (1.74–3.92%), aerobic ammonia oxidation (1.37–3.39%), and nitrogen fixation (1.06–1.79%). Eight functional microbial categories with relative abundance exceeding 1% were identified.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g008.jpg)
Figure 8. Prediction and analysis of soil bacterial community functions at various groundwater depths.
Subsequently, the functional characteristics of the rhizosphere bacterial communities in both plant species exhibited significant variations in response to GTD (Figure 9). In the rhizospheric soil of Artem, there was a notable increase in the relative abundance of bacteria engaged in carbon metabolism, aliphatic non-methane hydrocarbon degradation, ligninolysis, aromatic hydrocarbon degradation, and other functional categories under SW conditions. Conversely, functional bacteria associated with the nitrogen cycle, including ureolysis, aerobic ammonia oxidation, and nitrification, showed a marked elevation in their relative abundance under DW conditions.
In the Salix rhizospheric soil, there was a noteworthy increase in the relative abundance of functional bacteria involved in the carbon cycle, including chitinolysis and aromatic compound degradation, under SW conditions. Conversely, under DW conditions, a significant increase in relative abundance was observed for functional bacteria associated with phototrophy and nitrogen metabolism, including phototrophy, photoheterotrophy, photoautotrophy, anoxygenic photoautotrophy, aerobic chemoheterotrophy, chemoheterotrophy, nitrogen respiration, nitrate respiration, nitrite respiration, denitrification, nitrous oxide respiration, nitrate denitrification and nitrite denitrification.
The results revealed that both Artem and Salix exhibited a higher presence of functional bacteria related to carbon metabolism in SW. However, in DW, Artem demonstrated a greater abundance of functional microorganisms linked to nitrogen cycle, whereas Salix displayed a higher proportion of functional microorganisms involved in the nitrogen cycle and phototrophic processes.
3.3 Response of rhizosphere bacterial communities to plant–soil factors across various GTD
The Mantel test correlation analysis (Figure 10) revealed a significant positive correlation between plant and soil factors and the microbial community. Specifically, the SWC exhibited a remarkable association with the Shannon-Wiener diversity index of the Artem, which may exerted a prominent influence on bacterial alpha diversity (Figure 10A). This indicated that plant growth in SW was highly dependent on SWC. Moreover, enhanced plant diversity promoted greater involvement of functional microorganisms in facilitating plant nutrient uptake, thereby leading to an increase in alpha diversity within the rhizosphere bacterial community. Afterwards, plant diversity indices (Margalef, Simpson Wiener diversity, and Pielou evenness) demonstrated a substantial positive impact on functional bacteria that participate in the soil nitrogen cycle. Hence, SWC affected both rhizosphere bacterial diversity and the presence of functional microorganisms integral to the nitrogen cycle, all through its effects on plant growth under SW conditions. In contrast, in DW, except for the notable effect of SOC on microbial beta diversity, there was no correlation between plant and soil factors. Furthermore, these factors did not significantly affect the functional microorganisms involved in carbon and nitrogen cycles (Figure 10B). This implies that SOC exerts a more pronounced effect on the rhizospheric bacterial diversity of the Artem in DW conditions. And this observation result aligns with the identification of critical species, particularly the representative Actinobacteria, specifically ASV17731, from the Subgroup_17 class, known for their association with soil fertility and their role as decomposers of soil organic matter.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1285922/fmicb-14-1285922-HTML/image_m/fmicb-14-1285922-g010.jpg)
Figure 10. (A-D) show the shallow relationship between plant, soil factors and major functional microbes in the shallow groundwater table depth range (SW) and deep groundwater table range (DW) of Artem and Salix sites, respectively.
In Salix sites subjected to SW conditions, although plant and edaphic factors did not significantly affect microbial diversity and functional composition, a noteworthy correlation was detected between almost all soil–plant factors (except the Pielou index) and functional microorganisms engaged in the C cycle compared to the N cycle (Figure 10C). As far as Salix sites under DW conditions was concerned, not only plant diversity (Shannon diversity, Pielou evenness index) was strongly correlated with the C cycle functional microorganisms, but both plant diversity (Margalef, Simpson diversity index) and soil nutrients (pH, SOC, TN, AP and SWC) also had a significant effect on the N cycle (Figure 10D). These findings were similar with earlier microbial function predictions, affirming that rhizosphere microbes in Salix played a facilitative role in the carbon cycle under SW conditions. This could be attributed to representative species of Actinomycetes, specifically ASV192677 from the Subgroup_17 order, which contributed to soil fertility enhancement and acted as decomposers of soil organic matter. This, in turn, contributes to soil carbon cycling. Notably, these microbes made substantial contributions to both carbon and nitrogen cycling under DW conditions.
4 Discussion
4.1 Responses of plant diversity and soil physicochemical factors to GTD
The microtopographic features of the aquifer structure in Mu Us Sandy Land contributed to the variation in the GTD. As the GTD increased, the overall diversity of plant communities exhibited a decrease trend (Figure 2). This aligns with the study by Chen et al. (2006) in the Tarim River Basin, where increasing GTD led to a decrease in species diversity and richness, along with an overall low Pielou evenness index. In the study area, dominant vegetation primarily comprised shrubs, including Salix cheilophila, Artemisia desertorum, and Corethrodendron fruticosum, along with herbs, such as Phragmites australis, Thermopsis lanceolata, and Aster tataricus in the SW. In contrast, the DW area was characterized by shrubs such as Salix and Artem, accompanied by species such as Lespedeza bicolor and saplings of Pinus sylvestris var. mongolica. This pattern indicated fewer plant species and a simpler community structure in DW, reflecting soil water stress as the principal driver of changes in vegetation species diversity. This was attributed to the linear relationship between the GTD and SWC (Figure 3A), which is consistent with the findings of Zhang et al. (2018). In general, when the GTD is shallower, the root zone benefits from recharging via rising capillary water, thereby maintaining a higher soil water content. Conversely, with a deeper GTD, the soil receives less recharge from rising capillary water, resulting in lower soil water content. This indicates that groundwater and soil moisture dynamics are pivotal factors that affect vegetation growth (Huang et al., 2019; Barron-Gafford et al., 2021; Ma et al., 2022).
Previous studies have demonstrated that water efficiency in arid regions has notable implications for soil nutrient sustainability (Zeng J. et al., 2020; Condon et al., 2021; Reed et al., 2022). In both Artem and Salix sites, soil nutrients increased first and then decreased with the increase of GTD (Figures 3B–D). On the one hand, when the groundwater table is shallow, higher soil moisture levels reduce oxygen penetration into the soil. This reduces the activity of aerobic microorganisms involved in decomposition, resulting in a lower rate of organic matter breakdown and consequent soil nutrient accumulation (Sahrawat, 2003). On the other hand, under SW conditions, an adequate supply of vegetative nutrients is facilitated by hydraulic uplift within the root system (Prieto et al., 2012). This may explain why soil nutrient levels increased as GTD increased in SW. Conversely, as the groundwater table continued to decrease, reduced soil water content and heightened aeration in the vadose zone accelerated the decomposition of soil organic matter, leading to a decrease in SOC. This aligns with the trend reported by Zhao et al. (2023) for the Horqin sand plain, where nutrients exhibited a decline with increasing GTD. In this study, SOC, TN, and AP exhibited an initial increase, followed by a decrease in the GTD (Figure 3C). There are two potential explanations. First, TN and AP contents followed a similar trend to SOC with decreasing GTD, likely because 95% of TN and 40 to 60% of AP were derived from organic matter, thereby linking the decrease in SOC content to TN and AP levels (Wu et al., 2018). Second, this phenomenon may also be associated with the trend in vegetation diversity, which demonstrated a significant decrease trend (Figure 2), indirectly affecting the nutrient change. This implies that alterations in environmental factors, particularly groundwater table depth, in drylands exert varying influence on functional traits at the plant community level, potentially modifying their nutrient utilization and storage strategies (Reich, 2014). This also indicates that soil nutrient variations are contingent on GTD, plant cover, and their interactions (Li et al., 2021). In addition, SWC, SOC, TN, and AP levels in Salix surpassed those in Artem (Figure 3). Since Salix owns more extensive and robust root system, which confers greater efficiency in water and nutrient uptake than Artem (Chen et al., 2022). In arid sandy soils, shrubs typically exhibit well-developed root systems that mitigate water limitations (Xiao-Dong et al., 2022). This advantage facilitates nutrient enrichment around the plant root system through the ‘fertile island effect’, yielding pronounced spatial heterogeneity in soil physicochemical properties, particularly in nutrient content (Wang et al., 2019; Gao et al., 2022). Consequently, the groundwater table depth and vegetation patterns tend to influence nutrient accumulation in unsaturated soils (Zhang et al., 2022), thereby modifying the microenvironment, including factors such as pH, through alterations in soil redox processes and cation exchange rates (Supplementary Figure S1).
4.2 Response of rhizosphere bacterial community structure and co-occurrence networks to GTD
As depicted in Figure 4, the predominant phyla (with a relative abundance exceeding 1%) within the soil bacterial communities under varying GTDs for both plants included Actinobacteria, Proteobacteria, Chloroflexi, and Acidobacteria (Supplementary Table S2). The principal classes (relative abundance >1%) were Actinobacteria, Thermoleophilia, Alphaproteobacteria, and Gammaproteobacteria (Supplementary Table S3). Although the proportions of the major taxa varied, the major dominant taxa were basically the same as in other studies (Bi et al., 2021; Li et al., 2021). Similarly, the co-occurrence network within the rhizosphere bacterial communities of Artem and Salix revealed distinct compositions in the SW and DW. Specifically, in SW, the network was primarily composed of Proteobacteria and Acidobacteria, whereas Actinobacteria played a prominent role in DW (Figure 7). The observed impact of the GTD on the dominant bacterial phylum was notable. Specifically, alterations in the relative abundance of bacterial communities revealed a notable prevalence of the phylum Proteobacteria and its associated classes within the rhizosphere soil bacterial community (Figure 7). This phylum consists of four major groups (α-Proteobacteria, β-Proteobacteria, δ-Proteobacteria, and γ-Proteobacteria), all of which are Gram-negative, indicating their propensity to accumulate in the rhizosphere and facilitate organic matter decomposition. The presence of Acidobacteria, which produce extracellular polysaccharides that facilitate prokaryotic adhesion to root surfaces, demonstrates their positive interaction with plants, ultimately promoting plant growth (Kielak et al., 2016). In DW, the Artem co-occurrence network identified ASV187664, ASV17731, and ASV233705 from Bacillus of Firmicutes, Subgroup_17, and Catenulispora of Actinobacteria, respectively (Supplementary Figure S3). Gram-positive bacteria play pivotal roles in nutrient cycling processes, including organic matter decomposition and chemoenergetic heterotrophy. This suggests that the heightened GTD places Artem in an inhospitable environment, compelling key functional microorganisms involved in nutrient cycling to extract nutrients from soil. Zeng L. et al. (2020) corroborated that Actinobacteria and Acidobacteria, which function as saprophytic bacteria, actively participate in organic carbon degradation and contribute to carbon cycling, with Acidobacteria demonstrating a preference for promoting plant growth in nutrient-scarce soil environments (Fierer, 2007). In contrast, Salix exhibited selection for only one critical species from Actinobacteria, ASV192677, associated with soil fertility in SW (Figure 8), suggesting that Salix may display greater resilience to GTD variations than Artem in SW, consequently exerting a more pronounced influence on limiting GTD.
Under SW conditions, the relative abundance of Deltaproteobacteria, Anaerolineae, and NC10 class exhibited significant elevation compared to DW, indicating that specific families within Deltaproteobacteria, such as Comamonadaceae, are adept at utilizing nitrate-derived nitrogen for denitrification, thereby contributing to nitrogen cycling in SW (Khan et al., 2002). Furthermore, the relative abundance of Chloroflexia in the rhizosphere bacteria of Artem was notably higher in SW than in DW, whereas in Salix, the trend was reversed. Both Anaerolineae and Chloroflexia belong to the phylum Chloroflexi, characterized as anaerobic organisms that do not engage in oxygenic photosynthesis or nitrogen fixation. Instead, they fix carbon dioxide through the 3-hydroxypropionate pathway and are well suited for higher temperatures (Sprague et al., 1981). Therefore, the prevalence of Chloroflexia in SW of Artem may be attributed to its adaptation to elevated temperatures, whereas in DW of Salix, it may contribute to carbon fixation.
In the SW of Artem, there is relatively high microbial diversity and activity owing to the ample availability of essential nutrients for plant growth. Consequently, the soil bacterial community exhibited a more intricate network structure (Figure 7; Supplementary Table S4) and displayed heightened resilience to environmental fluctuations (de Vries et al., 2018; Tao et al., 2018; Luo et al., 2023). Furthermore, microbial networks characterized by increased co-occurrence and complexity tend to positively affect microbial behavior in nutrient cycling and subsequent ecosystem functioning (Wagg et al., 2019). This is attributed to their greater efficiency in resource utilization and information exchange (Faust, 2021; Ma et al., 2021). In contrast, for Salix, regardless of whether it was in SW or DW conditions, the complexity of the co-occurrence network did not exhibit significant differences. This implies that rhizosphere microorganisms in Salix demonstrated higher adaptability to changes in water table levels than those in Artem. Rhizosphere effect refers to the phenomenon that the exudates of plant roots and tissues provide abundant nutrients and energy for rhizosphere microorganisms. Therefore, this adaptability was likely attributed to the well-developed root system of Salix, which amplified the ‘rhizosphere effect’ in regulating microbial structure and influencing nutrient cycling (López et al., 2023). Furthermore, the co-occurrence network displayed a prevalence of positive connections over negative ones, indicating that microbial cooperation outweighs competition in adapting to the challenging conditions of arid desert ecosystems.
4.3 Regulation of soil nutrients by microbial community changes in Mu Us Sandy Land
In the SW of Artem sites (Figure 10A), further analysis can be conducted on the effects of plant diversity and soil factors on the diversity of rhizosphere microorganisms, as well as functional microorganisms engaged in carbon and nitrogen cycling. This indicates a potentially higher nitrogen content in the rhizosphere soils of Artem under SW conditions. In contrast, in DW conditions, there was minimal correlation between plant and soil factors, with only organic carbon significantly affecting microbial beta diversity. The substantial increase in GTD significantly affects various soil properties, including pH, SWC, and oxygen availability, thereby exerting selective pressure on the microbial community (Baselga, 2010). Consequently, microbial-mediated partitioning of carbon and nitrogen in the rhizosphere community becomes susceptible to change through alterations in critical species and functions crucial for sustaining material and energy production under adverse conditions (Yerbury et al., 2005; Bhattarai et al., 2015; Ghezzehei et al., 2019). For example, in response to drought stress, bacteria undergo increased synthesis of amino compounds to maintain intracellular osmotic pressure and prevent dehydration, thereby consuming substantial amounts of C and N resources (Schimel et al., 2007). Notably, the rate of carbon decomposition tends to increase with shallower groundwater tables, leading to heightened CO2 fluxes (Burkett and Kusler, 2000). This could potentially shift the status of the sandy soil from a carbon sink to a carbon source. This hypothesis is supported by the key species identified within the co-occurrence network, particularly the representative Actinomyces species ASV17731 and ASV233705. These species are associated with soil fertility and function as decomposers of soil organic matter. In addition, ASV187664 from Bacillus of phylum Firmicutes exhibits resilience in extreme environments and plays a role in respiratory metabolism and chemoenergetic heterotrophy.
In the case of Salix, nearly all plant–soil factors demonstrated robust correlations with functional bacteria involved in the C cycle at GTD < 5 m (Figure 10). This suggests that the rhizosphere soil of Salix in SW may offer a more conducive environment for the survival of functional microbes associated with the C cycle. Specifically, the representative species ASV192677 from the Actinomycete Subgroup_17 class played a role in enhancing soil fertility through its involvement in organic matter decomposition. This is in line with the findings of Bi et al. (2021) in the Yulin Botany Garden, where pH, a crucial regulatory factor governing microbial community distribution, significantly affected the dynamic processes of the soil bacterial community. Furthermore, SWC, pH, SOC, nutrients, and litter have been reported to be significantly correlated with microbial communities in the arid and semi-arid regions of northern China (Li X. Q. et al., 2019; Na et al., 2019; Yang et al., 2019). Notably, under DW conditions, plant diversity not only displayed a robust correlation with the C cycle, but both plant diversity and soil nutrients also exerted a notable effect on the N cycle. This was consistent with the results predicted by the microbial function (Figure 9), suggesting that GTD affected the dynamics of carbon and nitrogen cycling. Carbon nutrient availability may be related to functional microorganisms, whereas alterations in the nitrogen content may involve more intricate processes. Chen et al. (2012) observed that the lowering of the groundwater table induced vigorous nitrogen mineralization, leading to an accumulation of ammonium in the soil. This can provide more effective nitrogen for the soil and stimulate the participation of functional microorganisms in the nitrogen cycle. This further indicates an effect on nitrogen recovery in the root zone in the deeper groundwater tables. In addition, certain studies have indicated that limited oxygen supply in the soil can markedly affect gaseous nitrogen losses (specifically N2O and N2 emissions). Under highly anoxic conditions facilitated by a shallow groundwater table, there is a propensity for enhanced N2O consumption, consequently reducing the total N2O emissions (Hansen et al., 2014). Given that optimal conditions for complete denitrification of N2 are typically identified when the oxygen supply is restricted, the presence of Soil Water (providing an anaerobic environment) has been suggested to potentially diminish soil N2O emissions while augmenting N2 emissions (Flessa et al., 1998).
5 Conclusion
We conducted a study investigating the correlation between rhizosphere bacterial communities and soil nutrition in the Mu Us Sandy Land, focusing on two representative sandy plants subjected to varying GTD. Our findings revealed distinct characteristics between Salix and Artem sites. Salix sites exhibited higher SWC levels and soil nutrient concentrations than Artem sites, suggesting greater adaptability of Salix to fluctuations in groundwater table. Notably, the co-occurring network of functional bacteria in the rhizosphere of the Artem sites was more intricate in SW, displaying heightened resilience to the environment. This complexity exerts a substantial influence on regulating the microbial structure and affecting nutrient cycling. In contrast, Salix demonstrated a lower responsiveness to changes in GTD, attributed to its well-developed root system, consequently amplifying the ‘rhizosphere effect’ in governing microbial structure and influencing nutrient cycling. In terms of functional rhizosphere bacteria, both Artem and Salix positively affected carbon metabolism in SW. Conversely, in DW, the rhizosphere bacterial community of Artem was found to be more conducive to nitrogen cycling, whereas Salix exhibited a dual promotion of both nitrogen cycling and phototrophic functions. Importantly, we found that the rhizosphere soil of SW had a higher carbon content at Artem sites, as GTD increases, the rate of organic carbon decomposition accelerates, potentially causing a shift from carbon sink to carbon source. Salix sites in SW contain higher carbon levels and consumption of greenhouse gas N2O, as GTD increases, the soil undergoes significant nitrogen mineralization, increasing soil nitrogen effectiveness. These findings indicate the crucial role of microorganisms in mediating carbon and nitrogen redistribution in the soil to mitigate the adverse conditions caused by fluctuating groundwater table depths. It is imperative to acknowledge that groundwater tables have a notable impact on greenhouse gas emissions from sandy rhizosphere soils. Further research is warranted to broaden our understanding of the various vegetation types.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
LH: Data curation, Formal analysis, Methodology, Visualization, Writing – original draft. XL: Funding acquisition, Project administration, Supervision, Writing – review & editing. RJ: Data curation, Methodology, Visualization, Writing – review & editing. YM: Data curation, Funding acquisition, Investigation, Writing – review & editing. PW: Investigation, Resources, Writing – review & editing. QC: Investigation, Resources, Writing – review & editing. YX: Investigation, Resources, Writing – review & editing.
Funding
The author (s) declare financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by grants from the National Natural Science Foundation of China (42372288, 41877179, and 42001033), the Central University Special Fund (300102292904), the Natural Science Basic Research Program of Shaanxi (2023-JC-YB-280), the Shaanxi Forest Science and Technology Innovation Team (SXLK2023-0301), and the Central Fiscal Forestry Science and Technology Extension Project (SLTG202303). The Forestry Science and Technology Innovation Project of Shaanxi Province (SXLK2022-06-3).
Acknowledgments
The authors thanks to the National Forestry and Grassland Administration and Yulin Desert Ecological Station for providing long-term experimental sites and Yulin Municipal Water Resources Bureau for providing relevant groundwater depth data.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1285922/full#supplementary-material
Abbreviations
Artem Salix, Artemisia desertorum Salix cheilophila; GTD, Groundwater table depth; SW, shallow groundwater table depth range (GTD < 5 m); DW, deep groundwater table depth range (GTD > 5 m); SWC, Soil water content (%); SOC, Soil organic carbon
Footnotes
References
Barbeta, A., Mejia-Chang, M., Ogaya, R., Voltas, J., Dawson, T. E., and Penuelas, J. (2015). The combined effects of a long-term experimental drought and an extreme drought on the use of plant-water sources in a Mediterranean forest. Glob. Chang. Biol. 21, 1213–1225. doi: 10.1111/gcb.12785
Barron-Gafford, G. A., Knowles, J. F., Sanchez-Canete, E. P., Minor, R. L., Lee, E., Sutter, L., et al. (2021). Hydraulic redistribution buffers climate variability and regulates grass-tree interactions in a semiarid riparian savanna. Ecohydrology 14:e2271. doi: 10.1002/eco.2271
Baselga, A. (2010). Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 19, 134–143. doi: 10.1111/j.1466-8238.2009.00490.x
Bastian, M., Heymann, S., and Jacomy, M. (2009). "Gephi: an open source software for exploring and manipulating networks", in: Proceedings of the international AAAI conference on web and social media, San Jose
Bhattarai, A., Bhattarai, B., and Pandey, S. (2015). Variation of soil microbial population in different soil horizons. J. Microbiol. Exp. 2:00044. doi: 10.15406/jmen.2015.02.00044
Bi, B. Y., Zhang, H., Yuan, Y., Wu, Z. H., Wang, Y., and Han, F. P. (2021). Dynamic changes of soil microbial community in Pinus sylvestris var. mongolica plantations in the mu us Sandy land. J. Environ. Manag. 287:112306. doi: 10.1016/j.jenvman.2021.112306
Bissett, A., Richardson, A. E., Baker, G., Wakelin, S., and Thrall, P. H. (2010). Life history determines biogeographical patterns of soil bacterial communities over multiple spatial scales. Molecular Ecology. doi: 10.1111/j.1365-294x.2010.04804.x
Bremner, J. M., and Mulvaney, C. S. (1982). “Nitrogen-total” in Methods of soil analysis. Part 2. Chemical and microbiological properties. eds. A. L. Page, R. H. Miller, and D. R. Keeney (Madison, Wisconsin: American Society of Agronomy, Soil Science Society of America)
Burkett, V., and Kusler, J. (2000). Climate change: potential impacts and interactions in wetlands of the United States 1. JAWRA J. Am. Water Resour. Associat. 36, 313–320. doi: 10.1111/j.1752-1688.2000.tb04270.x
Castle, S. L., Thomas, B. F., Reager, J. T., Rodell, M., Swenson, S. C., and Famiglietti, J. S. (2014). Groundwater depletion during drought threatens future water security of the Colorado River Basin. Geophys. Res. Lett. 41, 5904–5911. doi: 10.1002/2014GL061055
Chao, A. (1984). “Nonparametric estimation of the number of classes in a population” in Scand. J. Statist, 11.
Chen, Y. T., Borken, W., Stange, C. F., and Matzner, E. (2012). Dynamics of nitrogen and carbon mineralization in a fen soil following water table fluctuations. Wetlands 32, 579–587. doi: 10.1007/s13157-012-0295-7
Chen, Y., Chen, Y., Xu, C., and Li, W. (2016). The effects of groundwater depth on water uptake of Populus euphratica and Tamarix ramosissima in the hyperarid region of northwestern China. Environ. Sci. Pollut. Res. 23, 17404–17412. doi: 10.1007/s11356-016-6914-8
Chen, Y. F., He, J. Q., He, Y., Gao, W. D., Zheng, C., and Liu, X. H. (2022). Seasonal hydrological traits in Salix psammophila and its responses to soil moisture and meteorological factors in desert areas. Ecol. Indic. 136:108626. doi: 10.1016/j.ecolind.2022.108626
Chen, X., Lu, J., Zhu, J., and Liu, C. (2020). Characteristics of denitrifying bacteria in different habitats of the Yongding River wetland, China. J. Environ. Manag. 275:111273. doi: 10.1016/j.jenvman.2020.111273
Chen, Y.-N., Zilliacus, H., Li, W.-H., Zhang, H.-F., and Chen, Y.-P. (2006). Ground-water level affects plant species diversity along the lower reaches of the Tarim river, Western China. J. Arid Environ. 66, 231–246. doi: 10.1016/j.jaridenv.2005.11.009
Christian, R., Alexander, W., Marcel, C., Yvonne, P., Christiane, R., Fredd, V., et al. (2022). Plant diversity effects on herbivory are related to soil biodiversity and plant chemistry. J. Ecol. 111, 412–427. doi: 10.1111/1365-2745.14032
Condon, L. E., Kollet, S., Bierkens, M. F. P., Fogg, G. E., Maxwell, R. M., Hill, M. C., et al. (2021). Global groundwater modeling and monitoring: opportunities and challenges. Water Resour. Res. 57:e2020WR029500. doi: 10.1029/2020WR029500
de Vries, F. T., Griffiths, R. I., Bailey, M., Craig, H., Girlanda, M., Gweon, H. S., et al. (2018). Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9:3033. doi: 10.1038/s41467-018-05516-7
Faust, K. (2021). Open challenges for microbial network construction and analysis. ISME J. 15, 3111–3118. doi: 10.1038/s41396-021-01027-4
Fierer, N. (2007). Finding a place for microorganisms in the field of ecology. Hoboken Wiley Online Library
Flessa, Wild, and Klemisch, Pfadenhauer. (1998). Nitrous oxide and methane fluxes from organic soils under agriculture. EUR J SOIL SCI.
Galland, W., Piola, F., Burlet, A., Mathieu, C., Nardy, M., Poussineau, S., et al. (2019). Biological denitrification inhibition (BDI) in the field: a strategy to improve plant nutrition and growth. Soil Biol. Biochem. 136:107513. doi: 10.1016/j.soilbio.2019.06.009
Gao, Y. J., Tariq, A., Zeng, F. J., Sardans, J., Zhang, Z. H., Islam, W., et al. (2022). “Fertile islands” beneath three desert vegetation on soil phosphorus fractions, enzymatic activities, and microbial biomass in the desert-oasis transition zone. Catena 212:106090. doi: 10.1016/j.catena.2022.106090
Ghezzehei, T. A., Sulman, B., Arnold, C. L., Bogie, N. A., and Berhe, A. A. (2019). On the role of soil water retention characteristic on aerobic microbial respiration. Biogeosciences 16, 1187–1209. doi: 10.5194/bg-16-1187-2019
Han, Z., Huang, S., Huang, Q., Bai, Q., Leng, G., Wang, H., et al. (2020). Effects of vegetation restoration on groundwater drought in the loess plateau, China. J. Hydrol. 591:125566. doi: 10.1016/j.jhydrol.2020.125566
Hansen, M., Clough, T. J., and Elberling, B. (2014). Flooding-induced N2O emission bursts controlled by pH and nitrate in agricultural soils. Soil Biol. Biochem. 69, 17–24. doi: 10.1016/j.soilbio.2013.10.031
Hill, R., Saetnan, E. R., Scullion, J., Gwynn-Jones, D., Ostle, N., and Edwards, A. (2015). Temporal and spatial influences incur reconfiguration of Arctic heathland soil bacterial community structure. Environ. Microbiol. doi: 10.1111/1462-2920.13017
He, J., Zhao, T., Chen, Y., Shang, X., Liu, S., Li, A., et al. (2023). Analysis on the influence of groundwater on the spatial pattern of vegetation coverage in mu us sand area. J. Soil Water Conserv. 37, 90–99. doi: 10.13870/j.cnki.stbcxb.2023.02.012, in Chinese
Huang, J., Yang, J., Jiang, H., Wu, G., Liu, W., Wang, B., et al. (2020). Microbial responses to simulated salinization and desalinization in the sediments of the Qinghai-Tibetan Lakes. Front. Microbiol. 11:1772. doi: 10.3389/fmicb.2020.01772
Huang, F., Zhang, D. R., and Chen, X. (2019). Vegetation response to groundwater variation in arid environments: visualization of research evolution, synthesis of response types, and estimation of groundwater threshold. Int. J. Environ. Res. Public Health 16:1849. doi: 10.3390/ijerph16101849
Jia, X. X., Shao, M. A., Wei, X. R., Zhu, Y. J., Wang, Y. Q., and Hu, W. (2020). Policy development for sustainable soil water use on China's Loess Plateau. Sci. Bull. 65, 2053–2056. doi: 10.1016/j.scib.2020.09.006
Jiakun, Y., Li, L., Wenhui, B., Suiqi, Z., and Ningning, Z. (2023). Phosphorus deficiency is the main limiting factor for re-vegetation and soil microorganisms in Mu Us Sandy Land, Northwest China. Science of the Total Environment
Karthe, D. (2018). Environmental changes in central and east Asian drylands and their effects on Major River-Lake systems. Quat. Int. 475, 91–100. doi: 10.1016/j.quaint.2017.01.041
Khan, S. T., Horiba, Y., Yamamoto, M., and Hiraishi, A. (2002). Members of the family Comamonadaceae as primary poly (3-hydroxybutyrate-co-3-hydroxyvalerate)-degrading denitrifiers in activated sludge as revealed by a polyphasic approach. Appl. Environ. Microbiol. doi: 10.1128/aem.68.7.3206-3214.2002
Kielak, A. M., Scheublin, T. R., and Mendes, L. W. (2016). Bacterial Community Succession in Pine-Wood Decomposition. Front Microbiol. 2016:231.
Kunkun, F., Hannah, H.-M., Corinne, W., Xisheng, G., Daozhong, W., Yang, B., et al. (2022). Identification of the rhizosphere microbes that actively consume plant-derived carbon. Soil Biol. Biochem. 166:108577. doi: 10.1016/j.soilbio.2022.108577
Lai, T. V., Ryder, M. H., Rathjen, J. R., Bolan, N. S., Croxford, A. E., and Denton, M. D. (2021). Dissimilatory nitrate reduction to ammonium increased with rising temperature. Biol. Fertil. Soils 57, 363–372. doi: 10.1007/s00374-020-01529-x
Li, Y. T., He, J. S., Wang, H., Zhou, J. Z., Yang, Y. F., and Chu, H. Y. (2021). Lowered water table causes species substitution while nitrogen amendment causes species loss in alpine wetland microbial communities. Pedosphere 31, 912–922. doi: 10.1016/S1002-0160(21)60023-1
Li, T., Jia, Y., Shao, M. A., and Shen, N. (2019). Camponotus japonicus ant mounds indirectly accelerate leaf litter decomposition by altering soil temperature and moisture. Soil Tillage Res. 194:104312. doi: 10.1016/j.still.2019.104312
Li, X. Q., Xia, J. B., Zhao, X. M., and Chen, Y. P. (2019). Effects of planting Tamarix chinensis on shallow soil water and salt content under different groundwater depths in the Yellow River Delta. Geoderma 335, 104–111. doi: 10.1016/j.geoderma.2018.08.017
Liang, Z., Olesen, J. E., Jensen, J. L., and Elsgaard, L. (2019). Nutrient availability affects carbon turnover and microbial physiology differently in topsoil and subsoil under a temperate grassland. Geoderma 336, 22–30. doi: 10.1016/j.geoderma.2018.08.021
Liu, L., Chen, H., Jiang, L., Hu, J., Zhan, W., He, Y., et al. (2018). Water table drawdown reshapes soil physicochemical characteristics in Zoige peatlands. Catena. doi: 10.1016/j.catena.2018.05.025
López, J. L., Fourie, A., Poppeliers, S. W. M., Pappas, N., Sánchez-Gil, J. J., de Jonge, R., et al. (2023). Growth rate is a dominant factor predicting the rhizosphere effect. ISME J. 17, 1396–1405. doi: 10.1038/s41396-023-01453-6
Luo, S., Png, G. K., Ostle, N. J., Zhou, H. K., Hou, X. Y., Luo, C. L., et al. (2023). Grassland degradation-induced declines in soil fungal complexity reduce fungal community stability and ecosystem multifunctionality. Soil Biol. Biochem. 176:108865. doi: 10.1016/j.soilbio.2022.108865
Ma, B., Stirling, E., Liu, Y. H., Zhao, K. K., Zhou, J. Z., Singh, B. K., et al. (2021). Soil biogeochemical cycle couplings inferred from a function-taxon network. Research 2021:7102769. doi: 10.34133/2021/7102769
Ma, W. J., Wei, F. L., Zhang, J. P., Karthe, D., and Opp, C. (2022). Green water appropriation of the cropland ecosystem in China. Sci. Total Environ. 806:150597. doi: 10.1016/j.scitotenv.2021.150597
Mao, W., Felton, A. J., Ma, Y., Zhang, T., Sun, Z., Zhao, X., et al. (2018). Relationships between aboveground and belowground trait responses of a dominant plant species to alterations in watertable depth. Land Degrad \u0026amp; Dev. doi: 10.1002/ldr.3159
Meyerhof, M.S., Wilson, J.M., Dawson, M.N., and Michael Beman, J. (2016). Microbial community diversity, structure and assembly across oxygen gradients in meromictic marine lakes, Palau. Environmental Microbiology
Na, X., Yu, H., Wang, P., Zhu, W., Niu, Y., and Huang, J. (2019). Vegetation biomass and soil moisture coregulate bacterial community succession under altered precipitation regimes in a desert steppe in northwestern China. Soil Biology and Biochemistry 136:107520. doi: 10.1016/j.soilbio.2019.107520
Nelson, D.W., and Sommers, L.E. (1982). Total carbon, organic carbon and organic matter, in: Methods of soil analysis part 2. Chemical and Microbial Properties. ASA-SSSA, Madison
Olsen, S. R., and Sommers, L. E. (1982). “Phosphorous, methods of soil analysis part 2, chemical and microbial properties” in Agronomy monograph. eds. A. L. Page, R. H. Miller, and D. R. Keeney, vol. 9 (Madison, Wisconsin: Agronomy Society of America), 403–430.
Osburn, E. D., Hoch, P. J., Lucas, J. M., McBride, S. G., and Strickland, M. S. (2022). Evaluating the roles of microbial functional breadth and home-field advantage in leaf litter decomposition. Funct. Ecol. 36, 1258–1267. doi: 10.1111/1365-2435.14026
Pielou, E. C. J. (1966). The measurement of diversity in different types of biological collections. J. Theor. Biol. 13, 131–144. doi: 10.1016/0022-5193(66)90013-0
Prieto, I., Armas, C., and Pugnaire, F. I. (2012). Hydraulic lift promotes selective root foraging in nutrient-rich soil patches. Funct. Plant Biol. 39, 804–812. doi: 10.1071/FP12070
Rao, P., Wang, Y., Liu, Y., Wang, X., Hou, Y., Pan, S., et al. (2022). A comparison of multiple methods for mapping groundwater levels in the Mu Us Sandy Land, China. J. Hydrol. 43:101189. doi: 10.1016/j.ejrh.2022.101189
Reed, C. C., Berhe, A. A., Moreland, K. C., Wilcox, J., and Sullivan, B. W. (2022). Restoring function: positive responses of carbon and nitrogen to 20 years of hydrologic restoration in montane meadows. Ecol. Appl. 32:e2677. doi: 10.1002/eap.2677
Reich, P. B. (2014). The world‐wide ‘fast–slow’ plant economics spectrum: a traits manifesto. J. Ecol. doi: 10.1111/1365-2745.12211
Schimel, J., Balser, T. C., and Wallenstein, M. (2007). Microbial stress-response physiology and its implications for ecosystem function. Ecology 88, 1386–1394. doi: 10.1890/06-0219
Schlaepfer, D. R., Bradford, J. B., Lauenroth, W. K., Munson, S. M., Tietjen, B., Hall, S. A., et al. (2017). Climate change reduces extent of temperate drylands and intensifies drought in deep soils. Nat. Commun. 8:14196. doi: 10.1038/ncomms14196
Shen, J., Liu, H., Zhou, H., and Chen, R. (2022). Specific characteristics of the microbial community in the groundwater fluctuation zone. Environ. Sci. Pollut. Res. 29, 76066–76077. doi: 10.1007/s11356-022-21166-1
Sprague, S. G., Staehelin, L. A., DiBartolomeis, M. J., and Fuller, R. C. (1981). Isolation and development of chlorosomes in the green bacterium Chloroflexus aurantiacus. J. Bacteriol. doi: 10.1128/jb.147.3.1021-1031.1981
Tao, J. M., Meng, D. L., Qin, C., Liu, X. D., Liang, Y. L., Xiao, Y. H., et al. (2018). Integrated network analysis reveals the importance of microbial interactions for maize growth. Appl. Microbiol. Biotechnol. 102, 3805–3818. doi: 10.1007/s00253-018-8837-4
Veldkornet, D. A., Potts, A. J., and Adams, J. B. (2016). The distribution of salt marsh macrophyte species in relation to physicochemical variables. S. Afr. J. Bot. doi: 10.1016/j.sajb.2016.08.008
Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E. E., and van der Heijden, M. G. A. (2019). Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat. Commun. 10:4841. doi: 10.1038/s41467-019-12798-y
Wall, D.H., Bardgett, R.D., Covich, A.R., and Snelgrove, P.V.R. (2004). "The need for understanding how biodiversity and ecosystem functioning affect ecosystem services in soils and sediments", D. H. Wall (Ed.), Sustaining biodiversity and ecosystem services in soils and sediments (SCOPE; Vol. 64). Washington D.C Island Press.
Wang, H. F., Cai, Y., Yang, Q., Gong, Y. M., and Lv, G. H. (2019). Factors that alter the relative importance of abiotic and biotic drivers on the fertile island in a desert-oasis ecotone. Sci. Total Environ. 697:134096. doi: 10.1016/j.scitotenv.2019.134096
Wang, Q., Zhang, H., Yang, Q., Wang, T., Zhang, Z., Liu, J., et al. (2022). The impact of grazing intensity on the allelopathic effect of Artemisia frigida in a temperate grassland in northern China. Flora 288:152005. doi: 10.1016/j.flora.2022.152005
Wu, B., and Ci, L. J. (2002). Landscape change and desertification development in the mu us Sandland, northern China. J. Arid Environ. 50, 429–444. doi: 10.1006/jare.2001.0847
Wu, J.-Q., Ma, W.-W., Li, G., Yan, L.-J., Zhu, H.-D., Luo, Y.-Z., et al. (2018). Comparative analysis of physicochemical property of soil with different groundwater level in Gahai Swamp Meadow Wetland. Acta Agrestia Sinica 26:341. doi: 10.11733/j.issn.1007-0435.2018.02.010
Xi, B., Di, N., Liu, J., Zhang, R., and Cao, Z. (2018). Hydrologic regulation of plant rooting depth: Pay attention to the widespread scenario with intense seasonal groundwater table fluctuation [Physical Sciences]. Proc. Natl. Acad. Sci. U.S.A. doi: 10.1073/pnas.1803987115
Xiang, Y., Shao, Z., Chai, H., Ji, F., and He, Q. (2020). Functional microorganisms and enzymes related nitrogen cycle in the biofilm performing simultaneous nitrification and denitrification. Bioresour. Technol. 314:123697. doi: 10.1016/j.biortech.2020.123697
Xiao-Dong, Y., Elhamjan, A., Jie, Z., Dong, H., Yong-Chao, G., Guang-Hui, L., et al. (2022). Higher association and integration among functional traits in small tree than shrub in resisting drought stress in an arid desert. Environ Exp Bot. doi: 10.1016/j.envexpbot.2022.104993
Yang, Y., Zhang, L., Wei, X., Chen, Y., Yang, W., Tan, B., et al. (2019). Litter removal reduced soil nitrogen mineralization in repeated freeze-thaw cycles. Scientific Reports. doi: 10.1038/s41598-018-38431-4
Yerbury, J. J., Stewart, E. M., Wyatt, A. R., and Wilson, M. R. (2005). Quality control of protein folding in extracellular space. EMBO Rep. 6, 1131–1136. doi: 10.1038/sj.embor.7400586
Zeng, J., Chen, H., Bai, Y., Dong, F., Peng, C., Yan, F., et al. (2020). Water table drawdown increases plant biodiversity and soil polyphenol in the Zoige Plateau. Ecol. Indic. 121:107118. doi: 10.1016/j.ecolind.2020.107118
Zeng, L., Li, J., Qin, K., Liu, J., Zhou, Z., and Zhang, Y. (2020). The total suitability of water yield and carbon sequestration under multi-scenario simulations in the Weihe watershed, China. Environ. Sci. Pollut. Res. Int. 27, 22461–22475. doi: 10.1007/s11356-020-08205-5
Zhang, X. L., Guan, T. Y., Zhou, J. H., Cai, W. T., Gao, N. N., Du, H., et al. (2018). Groundwater depth and soil properties are associated with variation in vegetation of a desert riparian ecosystem in an arid area of China. Forests 9:34. doi: 10.3390/f9010034
Zhang, B., Zeng, F. J., Gao, X. P., Shareef, M., Zhang, Z. H., Yu, Q., et al. (2022). Groundwater depth alters soil nutrient concentrations in different environments in an arid desert. Front. Environ. Sci. 10:939382. doi: 10.3389/fenvs.2022.939382
Zhao, M., Geruo, A., Zhang, J. E., Velicogna, I., Liang, C. Z., and Li, Z. Y. (2021). Ecological restoration impact on total terrestrial water storage. Nat. Sustain. 4, 56–62. doi: 10.1038/s41893-020-00600-7
Zhao, S. T., Zhao, X. Y., Li, Y. L., Chen, X. P., Li, C. Y., Fang, H., et al. (2023). Impact of deeper groundwater depth on vegetation and soil in semi-arid region of eastern China. Front. Plant Sci. 14:1186406. doi: 10.3389/fpls.2023.1186406
Keywords: groundwater table depth, plant diversity, rhizosphere functional bacteria, soil nutrient, sand-fixing plant
Citation: Hao L, Liu X, Ji R, Ma Y, Wu P, Cao Q and Xin Y (2023) Indirect regulation of topsoil nutrient cycling by groundwater depth: impacts on sand-fixing vegetation and rhizosphere bacterial communities. Front. Microbiol. 14:1285922. doi: 10.3389/fmicb.2023.1285922
Edited by:
Xuebin Qi, Chinese Academy of Agricultural Sciences, ChinaReviewed by:
Lei Deng, Northwest A&F University, ChinaYanhua Wang, Shaanxi Normal University, China
Xiaolin Liao, Nanjing Forestry University, China
Weiming Yan, Northwest A&F University, China
Copyright © 2023 Hao, Liu, Ji, Ma, Wu, Cao and Xin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiuhua Liu, liuxh68@chd.edu.cn