- 1Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Naiman Desertification Research Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Tongliao, China
- 4Tongliao Hydrology and Water Resources Sub-center, Tongliao, China
Understanding the effects of groundwater depth on soil microbiota and multiple soil functions is essential for ecological restoration and the implementation of groundwater conservation. The current impact of increased groundwater levels induced by drought on soil microbiota and multifunctionality remains ambiguous, which impedes our understanding of the sustainability of water-scarce ecosystems that heavily rely on groundwater resources. This study investigated the impacts of altered groundwater depths on soil microbiota and multifunctionality in a semi-arid region. Three groundwater depth levels were studied, with different soil quality and soil moisture at each level. The deep groundwater treatment had negative impacts on diversity, network complexity of microbiota, and the relationships among microbial phylum unites. Increasing groundwater depth also changed composition of soil microbiota, reducing the relative abundance of dominant phyla including Proteobacteria and Ascomycota. Increasing groundwater depth led to changes in microbial community characteristics, which are strongly related to alterations in soil multifunctionality. Overall, our results suggest that groundwater depth had a strongly effect on soil microbiota and functionality.
1. Introduction
The sufficiency of water resources is crucial for arid and semi-arid regions, as the availability of water is closely linked to vegetation productivity and energy cycling (Karthe, 2018; Ma et al., 2022). Direct and indirect accessibility of water are considered as potential mechanisms explaining plant drought tolerance (Barron-Gafford et al., 2021; Ma et al., 2022). For example, The ability of Prosopis tamarugo Phil. to survive in the extremely arid ecosystem of the Atacama Desert is attributed to its ability to extend roots to access groundwater (Garrido et al., 2016). Alterations of soil water resources availability is the dominant driver factor of restoration, sustainability, and interactions of soil nutrients and vegetation community (Chai et al., 2019; Condon et al., 2021; Reed et al., 2022). Drought can lead to soil fertility and nutrient loss, which reduces the support for microbial habitats and soil ecosystem services, ultimately resulting in decreased productivity (Wall et al., 2004; Schlaepfer et al., 2017; Jia et al., 2020). The simultaneous occurrence of multiple ecological functions in an ecosystem should not be considered in isolation. Therefore, the use of integrated multifunctional indicators, such as soil multifunctionality that includes soil nutrients, soil water content, soil physical properties, soil temperature, soil pH, etc., will enhance our comprehension and prediction of the functions supplied by soils and habitats, as well as the responsiveness to biodiversity loss and ecological changes (Fanin et al., 2018; Garland et al., 2021).
Soil microorganisms play a vital role in driving nutrition utilization and biogeochemical cycling, supporting restoration of ecosystem (Bardgett and van der Putten, 2014; Delgado-Baquerizo et al., 2016; Crowther et al., 2019). Soil microbial diversity, composition, interactions, and functions are responsive to soil ecological alterations (de Vries et al., 2018; Guo et al., 2018). Desiccation induced by increasing groundwater depth has been shown to strongly modify the diversity of soil microbiota and functionality in different ecosystems (Hu et al., 2019; Zhang et al., 2020; Wu H. et al., 2021). Meanwhile, the soil microbial community is highly complex, comprising tens of thousands of species per gram of soil, and it should be noted that microbial diversity alone is insufficient to fully explain microbial functions (Van Der Heijen, 2008). Analyzing the relationships between individual microbial communities and their functional groups in co-occurring networks can reveal the interdependence arising from soil environmental heterogeneity (Banerjee et al., 2019). Network analysis can suggest whether microbial groups are more critical for maintaining network stability (Freilich et al., 2018; Zhang M. et al., 2023). The use of relevant metrics like betweenness and closeness centrality in assessing the complexity of microbial networks allows for the elucidation of the interconnectivity among operational taxonomic units (OTUs), which can be attributed to the heterogeneity of soil environmental conditions (Banerjee et al., 2019; Wagg et al., 2019; Martinez, 2023). This analysis provides insights into the influence of environmental factors on soil microorganisms and sheds light on the dynamics of microbial complexity, as well as how microbial interactions affect ecosystem function. Understanding the correlations between soil microbiota, soil multifunctionality, and drought (caused by groundwater depth increase) is crucial for the sustainability of ecosystems, as soil microbes are sensitive to soil quality.
Horqin Sandy Land, an ecologically fragile region, is one of the largest four deserts in northern China, with an area of approximately 139,300 km2, of which the desertification area is as high as 71,884 km2 (Wang, 2016; Luo et al., 2017; Zuo et al., 2020). The geomorphological features of the region are characterized by alternating sand dunes and slightly undulating lowlands (Luo et al., 2020). In recent years, expedited population growth and the resulting boosted requirement for food, housing, employment, and land have led to more ecological problems (Kooch et al., 2022), especially the decline of the ground water level (Maihemuti et al., 2021). As the temperature increases in the future, the potential enhancement in evapotranspiration will be more pronounced than that of precipitation, leading to more water shortage and exacerbating desertification or desertification in this region (Dai and Zhao, 2017; Su et al., 2018; Zhang Y. et al., 2023). The future changes in climate conditions may enhance the vulnerability of sandy ecosystems and strongly impact the biotic and abiotic processes in the soil environment. Groundwater depth may be a determinant regulating the changes in soil multifunctionality, soil microbial diversity and composition. Consequently, the objectives of our study were to systematically analyze the soil microbiota and soil properties in the agro-pasture crisscross region of the Horqin Sandy Land, considering variations in groundwater depth. We aimed to assess the responses of these soil factors to different groundwater depths and provide support for ecological revegetation and the implementation of groundwater conservation. Our research sought to explore the following objectives:
1. investigate the effect of altered groundwater depth on soil multifunctionality, soil microbial diversity, and composition;
2. analyze the co-occurrence networks and microbial functionality in different groundwater depth; and
3. evaluate the relative contribution of different groundwater depth on soil microbiota characteristics, soil moisture, and multiple soil functions and their coupling relationship among the four.
2. Materials and methods
2.1. Study area
The study area is in the Horqin Sandy Land of eastern Inner Mongolia, China (42°55′N, 120°42′E). The elevation in our study site is 327 m. This region belongs to a typical temperate continental monsoon climate, with an annual average temperature of 5.8–6.4°C and an average annual precipitation of 351.7 mm. The spatial and temporal distribution of rainfall is uneven, with approximately 80% of the precipitation occurring from June to September (Huang et al., 2022).
The region of study is characterized by low soil nitrogen content, with levels ranging from 0.057 to 0.199 g/kg. The bulk density of the topsoil layer (0–30 cm) ranges from 1.29 to 1.59 g/cm3 (Mao et al., 2012). The dominant vegetation species of this region are Pennisetum centrasiaticum, Setaria viridis, Artemisia halodendron, and Caragana microphylla. This study selected Pennisetum centrasiaticum and Artemisia halodendron because they are capable of surviving in various habitats, especially in areas with severe water shortage, and can help to maintain the stability of soil structure and prevent land degradation.
2.2. Experimental design
This experiment began on June, 2017 at the Naiman Desertification Research Station (Su et al., 2021). As shown in Figure 1, three groundwater depth levels and two main species were used in the experimental manipulations (Pennisetum flaccidum Griseb. and Artemisia halodendron Turcz.et Bess.), with six replicates per each treatment. In 1 m2 concrete pools, an experimental operation was conducted, which included three treatment groups with groundwater depths of 50, 100, and 200 cm (Su et al., 2021). PVC shelters were set up to protect the cement pits from the rainfall. Establishment of a controlled experimental gradient was carried out to assess the effect of increasing groundwater depth on microbial diversity, network complexity, and multifunctionality of soil. Cement pits were buried to 1.0, 1.5, and 2.5 m according to each treatment gradient (Su et al., 2021). Stones with a diameter greater than 10 mm and a thickness of 10 cm were placed in cement pits, followed by sand to build the groundwater depth level. The bottom of each cement pit contained a 50 cm layer of stone–sand–water mixture. Adult individuals of Pennisetum flaccidum Griseb. and Artemisia halodendron Turcz.et Bess. were selected from the same sandy dune environment. The root system preserving the bud was excavated, and a root stem of equivalent diameter and length, possessing the bud, was selected. The selected rhizomes were then planted into the experimental pits.
Figure 1. Alpha diversity of soil bacterial and fungal as altered by groundwater depth and the correlations of soil multifunctionality to alpha diversity. Diversity of bacterial and fungal communities (Observed species, Shannon index and chao1) in soils from 50, 100, and 200 cm groundwater depth plots, and the correlations of soil multifunctionality to diversity index. Box plots exhibit the first (25%) and third (75%) quartiles, and median, and the maximum and minimum observed values. Asterisks denote significant differences based on the Wilcoxon test. *p < 0.05, **p < 0.01, and ***p < 0.001.
2.3. Soil sampling and analyses
Sampling was conducted in August 2022. This study measured 11 soil variables, including soil moisture (SM), total nitrogen (TN), nitrate nitrogen (NN), ammonium nitrogen (AN), soil organic carbon (SOC), total phosphorus (TP), and pH in the 0–50 cm soil layer. The soil characteristics, such as nutrient availability, biogeochemical cycling, and microbial productivity, were assessed using these variables, and standard methods were employed for the measurements (Perez-Harguindeguy et al., 2013). Soil moisture (SM) was determined by oven-drying soil samples at 105°C to a constant weight. Soil organic carbon (SOC) and total nitrogen (TN) were quantified using the Walkley-Black and Kjeldahl method, respectively. total phosphorus (TP) was quantified through colorimetric analysis following digestion with sulfuric acid and perchloric acid, employing standard procedures. nitrate–nitrogen (NN) and ammonium–nitrogen (AN) were analyzed using a continuous flow analyzer (AutoAnalyzer-AA3, Seal Analytical, Norderstedt, Germany) following extraction with 2 mol/L KCl. Soil pH was measured in a soil:water extract (1:2.5) with a pH meter (Mettler Toledo, Germany). The C:N ratio was calculated as the mass ratio of SOC and TN, while the C:P ratio and N:P ratio were calculated as the mass ratios of SOC and TP, and TN and TP, respectively. Soil multifunctionality (MF) was evaluated based on seven soil variables (TN, NN, AN, SOC, TP, SM, and pH) according to the method described by Fanin et al. (2018). This method has been widely employed in various studies (Maestre et al., 2012; Wang et al., 2019). Before analysis, the normality of the data was tested using the Shapiro–Wilk test. Non-normally distributed data was transformed by taking the logarithm or square root to approximate normality. For variables with negative values, the variable was transformed to positive values by subtracting the minimum value of the entire dataset. Variables were then standardized to a scale of 0 to 1 and the mean of these transformed values was taken as the functional value for each pit.
2.4. DNA extraction, illumina sequencing, and data processing
Soil microbes’ DNA was extracted using FastDNA Spin Kits (Boaojingdian BiotECH Co., Beijing, China) and the primers 341F (5-CTAYGGGRBGCASCAG-3) and 806R (5-GGACTACNNGGGTATCTAAT-3) were used to amplify the V3-V4 regions of the 16S rRNA gene (Yu et al., 2005). The fungal community was targeted using the broad-spectrum primers ITS1 (5-CTTGGTCATTTAGAGGAAGTAA-3) and ITS2 (5-GCTGCGTTCTTCATCGATGC-3; Li et al., 2020a). The PCR reaction was performed using Phusion® High-Fidelity PCR Master Mix (New England Biolabs) with a total volume of 15 μL, 0.2 μM of each forward and reverse primer, and approximately 10 ng of template DNA. The PCR cycling conditions included an initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, extension at 72°C for 30 s, and a final extension step at 72°C for 5 min (Liang et al., 2020). The amplicons obtained from triplicate PCR reactions were combined for each sample and analyzed by electrophoresis on a 2% (w/v) agarose gel. The PCR products were purified using a Qiagen Gel Extraction Kit (Qiagen, Germany). Libraries for sequencing, with sample index tags, were generated using the Illumina TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States). The quality of the library was evaluated using the Qubit® 2.0 Fluorometer (Thermo Scientific) and the Agilent Bioanalyzer 2100 system. The sequencing was performed on the Illumina NovaSeq platform, generating 250-bp paired-end reads at the Boaojingdian Company in Beijing, China. The Tax4Fun2 package was used for potential functional annotation and metabolic prediction. A total of 5,534,451 reads were obtained from 66 samples, ranging from 64,032 to 95,768 reads per sample.
2.5. Statistical analyses
The alpha diversity indices, namely the Shannon index, observed species, and Chao1, were computed using the Mothur software 8 (Montagna et al., 2018). The beta diversity among different groundwater depths was calculated using Mothur software, analyzed with Bray–Curtis dissimilarity matrix and ANOSIM function. The p-values were modified using Bonferroni correction in anosim function. The sorting of Bray-Curtis distance analysis (PCoA) was performed using the capscale function implemented in the vegan R package (Delgado-Baquerizo et al., 2019; Li et al., 2020b). Soil samples were gathered from different groundwater depths for co-occurrence network analysis. In each of the four groups, taxonomic units that constituted more than half of the samples were used to calculate Spearman correlation coefficients. Benjamini and Hochberg FDR correction was applied to control for false discovery. Statistically strong associations were detected with Spearman’s ρ > 0.8 and FDR-adjusted p < 0.05 and included in the subsequent network construction (Benjamini and Hochberg, 1995). The topological parameters of the network, including the number of nodes, edges, betweenness centrality, and closeness centrality, were calculated using the igraph package (Lu et al., 2018). The impacts of soil properties on bacterial and fungal community composition and diversity at different groundwater depths were assessed using random forest and structural equation models (SEMs) in “randomForest” and “piecewise” packages, respectively (Li et al., 2020a; Kong et al., 2022). First, the relative contributions of groundwater depth, soil properties, and multifunctionality to bacterial and fungal species richness and community composition were measured using a random forest model, to identify the main response variables. Subsequently, a piecewise structural equation model was constructed to examine the direct and indirect impacts of groundwater depth, soil water content, and soil multifunctionality on the diversity and composition of bacteria and fungi. The Shipley’s d-separation test was conducted to check any missing pathways in the model, and a value of p greater than 0.05 indicated no missing pathways. Standardized coefficients for each pathway in each component model were reported, as well as the Fisher’s C statistic, AIC value, and BIC value for the entire model (Shipley, 2009; Kong et al., 2022).
All statistical analyses were performed using IBM SPSS Statistics (version 24) and R 4.0.5 (R Development Core Team, 2021).
3. Results
Significant decrease in groundwater level affected the physicochemical properties of soil and reduce soil multifunctionality (Table 1). Under deep groundwater conditions, the TN and SOC content, C:P, N:P, pH, and moisture in soil were all higher compared to shallow groundwater depths (p < 0.05). The content of NN increased initially and then declined as groundwater depth increased (p < 0.05). However, AN content was increased with the increasing groundwater depth (p < 0.05). Groundwater depth was negatively related to soil multifunctionality and soil multifunctionality was lower in deeper groundwater depth (Table 1).
Soil microbial diversity was significantly altered by increasing groundwater depth in bacteria and fungi (Figure 1). The observed bacterial and fungal species richness, Shannon index, and Chao1 were significantly lower in deeper groundwater depth than in shallower depth. Nevertheless, the association among soil microbial diversity and MF was significant and positive (p < 0.05). Thus, losses of soil microbial diversity in altered groundwater depth were related with downward transitions in soil environmental functions.
Groundwater depth had a significant effect on the structure of soil microbial communities (Figure 2). Principle coordinate analysis (PCoA) suggested that soil microbiota strongly clustered following groundwater depth at bacteria level, which explained 12.71–35.42% of the total variation (Figure 2). Our results demonstrated that groundwater depth heavily affected the soil microbiota.
Figure 2. Soil microbial composition as influenced by groundwater depth on bacteria and fungi levels. Principal Coordinate Analysis (PCoA) plot of Bray–Curtis distances among groundwater depths (GWD = 50, 100, 200 cm); the relative abundance (%) of the phyla and family levels in soil dominant bacterial and fungal communities in three groundwater depth treatments.
At the phyla level, the bacterial communities in the soil were primarily composed of Proteobacteria (29.86%), Actinobacteria (20.14%), Firmicutes (5.98%), Actinobacteriota (12.04%), Acidobacteriota (9.81%), and Chloroflexi (6.18%), while about 15.42–16.57% of the bacteria were unidentified; soil fungal communities were analyzed and found to be mainly composed of members of Ascomycota (65.18%), Basidiomycota (13.85%), Mortierellomycota (15.14%), Chytridiomycota (3.77%), Rozellomycota (0.71%), Glomeromycota (1.01%), and Olpidiomycota (0.19%), and these phyla accounted for more than 85% of the relative abundance of all identified phyla (Figure 2). With increasing groundwater depth, the relative abundance of Proteobacteria decreases in bacteria, while the relative abundance of Ascomycota decreases in fungi.
At the family level, Sphingomonadaceae (26.22%), Nocardioidaceae (23.71%), Bacillaceae (12.56%), Micrococcaceae (12.28%), and Pyrinomonadaceae (13.08%) emerged as a dominant family in the bacterial communities, constituting over 85% of the total bacterial sequences (Figure 2). The relative abundance of Sphingomonadaceae, Nocardioidaceae, and Bacillaceae was declined first then increased with the increasing groundwater depth. The dominant fungal family was Nectriaceae (44.71%), which was followed by Mortierellaceae (26.29%), Aspergillaceae (12.19%), Marasmiaceae (5.78%), Clavicipitaceae (3.96%), Agaricaceae (3.61%), and Cordycipitaceae (3.44%). The relative abundance of Nectriaceae and Aspergillaceae decreased strongly with growing groundwater depth, whereas the relative abundance of Mortierellaceae and Marasmiaceae boosted significantly with increasing groundwater depth.
The soil microbial community network exhibits apparent co-occurrence patterns (Figure 3). Groundwater depth had a negative influence on the complexity of soil bacterial co-occurring network, but increased the complexity of fungi (Figures 3C,D). For bacteria, the betweenness and closeness centrality were significantly greater in deeper depth than others, and the closeness centrality of fungi was significantly greater in 100 cm groundwater depth, indicating that the soil microbial network was less robust in the deeper groundwater depth than in the shallower (Figures 3E,F). These results suggested that increased groundwater depth impacted microbial correlations, and decreased the complication of soil bacterial community networks. Additionally, the effects of groundwater depth on soil microbial co-occurrence network were closely related to community composition and sensitive strains.
Figure 3. Co-occurrence networks for bacteria and fungi in different groundwater depth. Nodes in the network correspond to amplicon sequence variants (ASVs), and the color and size of each node indicate its phylum affiliation and the degree in the network. Red and green edges represent positive and negative interactions, respectively. Network properties for bacteria (A) and fungi (B), comprised of node number (C), edge number (D), betweenness (E), and closeness (F). Asterisks denote significant differences based on the Wilcoxon test. *p < 0.05, **p < 0.01 and, ***p < 0.001.
Forty-two potential metabolic functional groups of the bacterial communities were predicted using Tax4Fun2 package based on KEGG pathway (Level 2) analysis (Figure 4). Increased groundwater depth decreased the values of most bacterial metabolic functional groups. For instance, the relative abundances of groups related to Amino acid metabolism, Biosynthesis of other secondary metabolism, Energy metabolism, Nucleotide metabolism, and Xenobiotics biodegradation and metabolism decreased after increased groundwater depth (p < 0.05). Meanwhile, most bacterial functional groups showed a strong negative relationship with soil properties.
Figure 4. Variation in bacterial functions with altered groundwater depth. The relative abundance of the 44 bacterial functions predicted by Tax4Fun2-KEGG pathway (level 2). The Spearman’s correlation between predicted bacterial functions and soil properties.
Derived from the results obtained from FUNGuild, the affected fungal functions are displayed in Figure 5. The relative abundance of fungal function which associated to Pathotroph decreased with increasing groundwater depth (p < 0.05). However, increased of groundwater depth improved the relative abundance of group associated to Saprotroph function (p < 0.05). Moreover, the relative abundance of fungal function related to Pathotroph was positively related to soil chemical properties (p < 0.05), while the relative abundance of Saprotroph group was negatively related to soil chemical properties (p < 0.05).
Figure 5. Variation in functions of fungal communities with altered groundwater depth. The relative abundance of the 10 functions predicted by FUNGuild. The Spearman’s correlation between predicted fungal functions and soil properties.
Random forest analyses were employed to identify the key properties and multifunctionality of soil measures that were critical in predicting the soil microbial composition and diversity (Figure 6). In bacteria community, random forest models explained 25.47% of the variance in diversity and 12.46% in composition. The models based on random forest analyses were able to account for 14.75 and 10.75% of the variance in diversity and fungal composition, respectively. AN, N:P, and SOC were found to be the main predictors of bacterial diversity, while N:P, C:N, SM, MF, and TN were identified as the dominant predictors of bacterial composition. The main predictors of fungal diversity were found to be AN, MF, pH, SOC, and TP, while pH, NN, SM, and TN were the main predictors of fungal composition.
Figure 6. Primary predictors of the soil bacterial and fungal composition and diversity rely on random forest models. MSE, mean squared error; TN, total nitrogen; NN, nitrate nitrogen; AN, ammonium nitrogen; SOC, soil organic carbon; TP, total phosphorus; SM, soil moisture; MF, multifunctionality of soils.
Structural equation modeling demonstrated that alterations in soil properties induced by changes in groundwater depth had a primary influence on the microbial diversity and composition (Figure 7). In contrast, the microbial composition was not significantly impacted by changes in soil moisture and multifunctionality (Figures 7B,D). Moreover, soil moisture and soil multifunctionality had significant and positive effects on microbial community diversity (Figures 7A,C). Increased of groundwater depth significantly reduced soil moisture and soil multifunctionality, whereas it had not direct effects on microbial community diversity and composition. Taken together, the increase in groundwater depth, which induces soil desiccation and degradation, can contribute to a decrease in microbial community diversity and composition in this semi-arid region.
Figure 7. Impacts of groundwater depth, soil moisture, and soil multifunctionality on the diversity (A,C) and composition (B,D) of soil bacterial (A,B) and fungal (C,D) communities rely on a piecewise structural equation model (SEM). Solid lines indicate significant paths (*p < 0.05; **p < 0.01); dashed lines indicate nonsignificant paths (p ≥ 0.05); black lines indicate positive effects; red lines indicate negative effects.
4. Discussion
4.1. microbial diversity and soil multifunctionality
Our research on soil properties and microbial diversity with contrasting three groundwater depths revealed a decline in soil biodiversity and multifunctionality induced by the decreasing groundwater level. This finding has been extensively documented that the increase in groundwater depth can result in the reduction of soil nutrients, degradation of soil structure, reduction in soil moisture content, and decline in soil functionality, all of which have a detrimental impact on soil microbial diversity and composition (Prieto et al., 2012; Sun et al., 2020; Chen et al., 2021). Whereas, the impact of rising groundwater depth on soil microbiota has been relatively overlooked (Zhang et al., 2022), and it has not been clearly defined whether the alterations in soil microbial diversity caused by increasing groundwater depth are associated with a decline in multiple soil functions. The depletion in soil total nitrogen, soil organic carbon, and soil moisture by increasing groundwater depth could directly resulted in the loss of microbial diversity because reduced soil resource availability constrains the microbial metabolism and composition, furthermore decreased their supports on soil multifunctionality (Chen et al., 2020; Ye et al., 2023). Meanwhile, a previous study suggested that decreased soil organic matter content may improve soil thermal conductivity and a decline in soil heat capacity, then increase the daily variance of soil temperature (Abu-Hamdeh and Reeder, 2000; Balashov et al., 2023). Therefore, the indirect increase in soil thermal variability due to the deficiency of microbial diversity and multifunctionality is well-established, as many soil microbes are responsive to alterations in soil temperature (Karimi et al., 2018; Li J. et al., 2023).
4.2. Composition of soil bacteria and fungi
At the phyla level, the result showed that the relative abundance of Proteobacteria in bacteria decreased with increasing groundwater depth. The phylum Proteobacteria was found to be the most dominant bacterial group in the majority of the soil samples collected in our study. Zhang et, al (2016) found that the Proteobacteria is highly abundant and comprises fast-growing bacteria that can use a variety of resources, which are typically considered copiotrophic (Zhang et al., 2016). Meanwhile, the results indicated that the relative abundance of Actinobacteria increased with increasing groundwater depth. We also found increasing of groundwater depth significantly reduced soil moisture of 0–50 cm layers. The relative abundance of Actinobacteria was reported that decrease by the addition of water (Van Horn et al., 2014; Monteiro et al., 2023). In the fungal composition analysis, we observed a decline in the relative abundance of Ascomycota with increasing groundwater depth. The Ascomycota had previously been reported to have a strong positive correlation with the metabolism of organic substrates in rhizodeposition (Liu et al., 2020; Zhong et al., 2020). Consequently, the loss of available substrate and nutrients reduced the relative abundances of these microbial community phyla which was associated to increasing of groundwater depth. The relative abundance of Sphingomonadaceae increased with increasing groundwater depth at the family level. Our results also suggested that soil pH was negatively correlated with groundwater depth. A previous study demonstrated that Sphingomonadaceae was able to secrete dehydrogenase to balance soil pH (Gong et al., 2021). Increasing groundwater depth improve the relative abundance of Mortierellaceae in fungal community. The observed response pattern is likely attributed to the longer plant-water transport distance associated with deeper groundwater depth, which may prompt the plant to acquire available water resources by increasing the proportion of roots. The Mortierellaceae have previously been reported to have a positive relationship with root growth (Li et al., 2017). Overall, this study suggested that the alterations in soil microbial composition following an increase in groundwater depth were closely linked to the response of soil environmental conditions.
4.3. Microbial co-occurrence network
Our study also showed that the decline of groundwater level reduced the complexity of microbial co-occurring network. The observed decrease in soil microbial diversity may be attributed to the increasing groundwater depth leading to soil desiccation. Drought conditions are known to decrease microbial complexity as a result of their lower resistance to abiotic stress, and thus the loss of soil moisture caused by increasing groundwater depth could be a key driver in the observed decrease in soil microbial diversity (de Vries et al., 2018; Banerjee et al., 2019; Li Z. et al., 2023). Many studies have demonstrated that the complexity of soil microbial networks is positively associated with soil resources, such as soil moisture and fertility, and this interpretation aligns with the observed alterations in microbial composition and the response to soil conditions identified in our study (Guo et al., 2020; Hernandez et al., 2021; Kong et al., 2022; Zou et al., 2023). Besides, the robustness of microbial co-occurrence network was previously reported to be negatively characterized by betweenness and closeness (Figures 3C–F; Wu M.-H. et al., 2021). Our study showed that the responsive taxa of deeper groundwater depth had higher betweenness and closeness. This result suggested that these taxa may be responsible for driving the alterations in the microbial co-occurrence network (de Vries et al., 2018; Yang et al., 2023). Taken together, the results of our study provide additional evidence supporting the notion that deeper groundwater depth could have an indirect effect on the soil microbial co-occurrence network via soil desiccation.
4.4. Microbial functionality
For bacterial functionality, a decrease in the relative abundance of functional groups associated to amino acid metabolism and carbohydrate metabolism was showed in Figure 4. Many studies suggested that soil drought by increasing groundwater depth impaired these functional groups (Chen et al., 2012; Bromke, 2013; Kong et al., 2022). At a groundwater depth of 50 cm, the increase in soil moisture and soil nutrients could potentially mitigate the negative impact of drought on soil bacterial function. However, for deeper groundwater depth, drought caused by soil water stress may restrict microbial growth and hinder the availability of nutrients to bacteria by directly constraining bacterial dispersal within the pores of arid and sandy soil (Giannetta et al., 2018). We also found that the relative abundance of fungal Pathotroph decreased with increasing groundwater depth. The relative abundance of Pathotroph were suggested to be related with plant and root biomass, and deeper groundwater depth restrained growth of herb root (Levi et al., 2022). Although the relative abundance of Saprotroph was higher in the treatment with deeper groundwater depth, it is likely attributed to the increased soil heat capacity caused by drought. Consequently, the enhancement of fungal functionality, which is temperature-sensitive, may be affected (Qu et al., 2021; Li et al., 2022).
5. Conclusion
In summary, we explored the connections between groundwater depth, soil microbial diversity, microbial composition, soil moisture, and soil multifunctionality using the same soil and species at the same site in a semi-arid region. We found that soil moisture and multifunctionality changed by groundwater depth had different effects on soil microbial community (Figure 7). Careful consideration of such changes was necessary because in most arid and semi-arid climates, groundwater depth exceeds 50 cm (Chen et al., 2021; Trinidad Torres-Garcia et al., 2022). The results presented in our study suggest that deeper groundwater depth led to a depletion in soil moisture and soil properties, as well as a shift in microbial community structure and functionality. Essentially, soil microorganisms play a crucial role in supporting soil multifunctionality by promoting decomposition, nutrient cycling, and resource availability in different microenvironments (Delgado-Baquerizo et al., 2016, 2020). The present study highlights the significance of soil microbial communities in facilitating multifunctionality in arid environments, with diverse relationships between alterations in soil microbial diversity and composition and modifications in soil quality triggered by changes in groundwater depth. Additionally, as complex microbial communities are typically more resilient to environmental challenges than simpler ones, the observed reduction in microbial diversity and composition resulting from increased groundwater depth may have long-term undesirable impacts on soil functions (Fierer et al., 2003; Allison and Martiny, 2008; Crowther et al., 2015).
The importance of soil moisture and nutrient availability for the growth of drought-tolerant plants, especially in semi-arid and arid regions, where soil moisture plays a more significant role than rainfall, is widely recognized (McCulley et al., 2004; Li et al., 2021). Besides, the growing depth of groundwater is closely related to species mortality in response to drought (Goulden and Bales, 2019). Consequently, changes in water management may have a stronger effect on microbiota than changes in environmental factors. Furthermore, our results suggested that increasing groundwater depth which reduced soil moisture impaired soil microbial functionality. In semi-arid and arid regions, soil microbiota played a crucial role in regulating the availability of soil nutrients and water to plant roots. This study revealed that the reduction in dominant species in semi-arid and arid regions may be attributed not only to the depletion in soil water content but also to the impairment of multiple microbial functions resulting from drought (Peng et al., 2011; Fu et al., 2017; Stewart et al., 2017). Therefore, reduction of soil multifunction, soil water content, and soil microbial community caused by increasing groundwater depth together contribute to desertification. The investigation of the restoration of the soil microbiota and its functions in sandy soil under future groundwater conservation strategies is necessary.
Data availability statement
The raw data has been successfully released with accession number PRJNA968046 and is available at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA968046.
Author contributions
SZ: conceptualization, data curation, formal analysis, investigation, methodology, validation, visualization, writing—original draft, and writing—review and editing. XZ: funding acquisition, project administration, resources, supervision, validation, visualization, writing—original draft, and writing—review and editing. YL and RZ: conceptualization, formal analysis, investigation, methodology, resources, supervision, validation, visualization, writing—original draft, and writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (No. 42177456), Transformation Program of Scientific and Technological Achievements of Inner Mongolia Autonomous Region (No. 2021CG0012), and the National Project on Science and Technology Basic Resources Survey of China (No. 2017FY100200).
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.
References
Abu-Hamdeh, N. H., and Reeder, R. C. (2000). Soil thermal conductivity: Effects of density, moisture, salt concentration, and organic matter. Soil Sci. Soc. Am. J. 64, 1285–1290. doi: 10.2136/sssaj2000.6441285x
Allison, S. D., and Martiny, J. B. H. (2008). Resistance, resilience, and redundancy in microbial communities. Proc. Natl. Acad. Sci. U. S. A. 105, 11512–11519. doi: 10.1073/pnas.0801925105
Balashov, E. V., Dobrokhotov, A. V., and Kozyreva, L. V. (2023). Effects of thermal and hydrophysical properties of sandy Haplic Podzol on actual evapotranspiration of spring wheat. J. Hydrol. Hydromech. 71, 125–131. doi: 10.2478/johh-2023-0013
Banerjee, S., Walder, F., Buechi, L., Meyer, M., Held, A. Y., Gattinger, A., et al. (2019). Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. ISME J. 13, 1722–1736. doi: 10.1038/s41396-019-0383-2
Bardgett, R. D., and Van Der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi: 10.1038/nature13855
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
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate-a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57, 289–300.
Bromke, M. A. (2013). Amino acid biosynthesis pathways in diatoms. Meta 3, 294–311. doi: 10.3390/metabo3020294
Chai, Y., Cao, Y., Yue, M., Tian, T., Yin, Q., Dang, H., et al. (2019). Soil abiotic properties and plant functional traits mediate associations between soil microbial and plant communities during a secondary forest succession on the Loess Plateau. Front. Microbiol. 10:895. doi: 10.3389/fmicb.2019.00895
Chen, Q.-L., Ding, J., Zhu, D., Hu, H.-W., Delgado-Baquerizo, M., Ma, Y.-B., et al. (2020). Rare microbial taxa as the major drivers of ecosystem multifunctionality in long-term fertilized soils. Soil Biol. Biochem. 141:107686. doi: 10.1016/j.soilbio.2019.107686
Chen, L., Yang, Y., Deng, S., Xu, Y., Wang, G., and Liu, Y. (2012). The response of carbohydrate metabolism to the fluctuation of relative humidity (RH) in the desert soil cyanobacterium Phormidium tenue. Eur. J. Soil Biol. 48, 11–16. doi: 10.1016/j.ejsobi.2011.10.002
Chen, G., Yue, D., Zhou, Y., Wang, D., Wang, H., Hui, C., et al. (2021). Driving factors of community-level plant functional traits and species distributions in the desert-wetland ecosystem of the Shule River Basin, China. Land Degrad. Dev. 32, 323–337. doi: 10.1002/ldr.3624
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
Crowther, T. W., Thomas, S. M., Maynard, D. S., Baldrian, P., Covey, K., Frey, S. D., et al. (2015). Biotic interactions mediate soil microbial feedbacks to climate change. Proc. Natl. Acad. Sci. U. S. A. 112, 7033–7038. doi: 10.1073/pnas.1502956112
Crowther, T. W., Van Den Hoogen, J., Wan, J., Mayes, M. A., Keiser, A. D., Mo, L., et al. (2019). The global soil community and its influence on biogeochemistry. Science 365:772.
Dai, A. G., and Zhao, T. B. (2017). Uncertainties in historical changes and future projections of drought. Part I: estimates of historical drought changes. Clim. Chang. 144, 519–533. doi: 10.1007/s10584-016-1705-2
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
Delgado-Baquerizo, M., Bardgett, R. D., Vitousek, P. M., Maestre, F. T., Williams, M. A., Eldridge, D. J., et al. (2019). Changes in belowground biodiversity during ecosystem development. Proc. Natl. Acad. Sci. U. S. A. 116, 6891–6896. doi: 10.1073/pnas.1818400116
Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 7:10541. doi: 10.1038/ncomms10541
Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., Abades, S., Alfaro, F. D., et al. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat. Ecol. Evol. 4, 210–220. doi: 10.1038/s41559-019-1084-y
Fanin, N., Gundale, M. J., Farrell, M., Ciobanu, M., Baldock, J. A., Nilsson, M.-C., et al. (2018). Consistent effects of biodiversity loss on multifunctionality across contrasting ecosystems. Nat. Ecol. Evol. 2:269.
Fierer, N., Schimel, J. P., and Holden, P. A. (2003). Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 35, 167–176. doi: 10.1016/S0038-0717(02)00251-1
Freilich, M. A., Wieters, E., Broitman, B. R., Marquet, P. A., and Navarrete, S. A. (2018). Species co-occurrence networks: Can they reveal trophic and non-trophic interactions in ecological communities? Ecology 99, 690–699. doi: 10.1002/ecy.2142
Fu, B., Wang, S., Liu, Y., Liu, J., Liang, W., and Miao, C. (2017). “Hydrogeomorphic ecosystem responses to natural and anthropogenic changes in the Loess Plateau of China” in Annual Review of Earth and Planetary Sciences, Vol 45. eds. R. Jeanloz and K. H. Freeman
Garland, G., Banerjee, S., Edlinger, A., Miranda Oliveira, E., Herzog, C., Wittwer, R., et al. (2021). A closer look at the functions behind ecosystem multifunctionality: A review. J. Ecol. 109, 600–613. doi: 10.1111/1365-2745.13511
Garrido, M., Silva, P., and Acevedo, E. (2016). Water Relations and Foliar Isotopic Composition of Prosopis tamarugo Phil., an Endemic Tree of the Atacama Desert Growing at Three Levels of Water Table Depth. Front. Plant Sci. 7:375. doi: 10.3389/fpls.2016.00375
Giannetta, B., Plaza, C., Vischetti, C., Cotrufo, M. F., and Zaccone, C. (2018). Distribution and thermal stability of physically and chemically protected organic matter fractions in soils across different ecosystems. Biol. Fertil. Soils 54, 671–681. doi: 10.1007/s00374-018-1290-9
Gong, W.-J., Niu, Z.-F., Wang, X.-R., and Zhao, H.-P. (2021). How the soil microbial communities and activities respond to long-term heavy metal contamination in electroplating contaminated site. Microorganisms 9:362. doi: 10.3390/microorganisms9020362
Goulden, M. L., and Bales, R. C. (2019). California forest die-off linked to multi-year deep soil drying in 2012-2015 drought. Nat. Geosci. 12:632.
Guo, X., Feng, J., Shi, Z., Zhou, X., Yuan, M., Tao, X., et al. (2018). Climate warming leads to divergent succession of grassland microbial communities. Nat. Clim. Chang. 8:813.
Guo, J., Ling, N., Chen, Z., Xue, C., Li, L., Liu, L., et al. (2020). Soil fungal assemblage complexity is dependent on soil fertility and dominated by deterministic processes. New Phytol. 226, 232–243. doi: 10.1111/nph.16345
Hernandez, D. J., David, A. S., Menges, E. S., Searcy, C. A., and Afkhami, M. E. (2021). Environmental stress destabilizes microbial networks. ISME J. 15, 1722–1734. doi: 10.1038/s41396-020-00882-x
Hu, Y., Zhang, Z., Huang, L., Qi, Q., Liu, L., Zhao, Y., et al. (2019). Shifts in soil microbial community functional gene structure across a 61-year desert revegetation chronosequence. Geoderma 347, 126–134. doi: 10.1016/j.geoderma.2019.03.046
Huang, W.-D., He, Y.-Z., Wang, H.-H., and Zhu, Y.-Z. (2022). Leaf Physiological Responses of Three Psammophytes to Combined Effects of Warming and Precipitation Reduction in Horqin Sandy Land, Northeast China. Front. Plant Sci. 12:785653. doi: 10.3389/fpls.2021.785653
Jia, X., Shao, M. A., Wei, X., Zhu, Y., Wang, Y., 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
Karimi, B., Terrat, S., Dequiedt, S., Saby, N. P. A., Horriguel, W., Lelievre, M., et al. (2018). Biogeography of soil bacteria and archaea across France. Science. Advances 4:eaat1808. doi: 10.1126/sciadv.aat1808
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
Kong, W., Wei, X., Wu, Y., Shao, M., Zhang, Q., Sadowsky, M. J., et al. (2022). Afforestation can lower microbial diversity and functionality in deep soil layers in a semiarid region. Glob. Chang. Biol. 28, 6086–6101. doi: 10.1111/gcb.16334
Kooch, Y., Ghorbanzadeh, N., Kuzyakov, Y., Praeg, N., and Ghaderi, E. (2022). Investigation of the effects of the conversion of forests and rangeland to cropland on fertility and soil functions in mountainous semi-arid landscape. Catena 210:105951. doi: 10.1016/j.catena.2021.105951
Levi, M., Applebaum, I., Sherman, C., Doniger, T., and Steinberger, Y. (2022). Soil fungal community of wheat Triticum aestivum rhizosphere at different phenological stages under a rain-fed management. Rhizosphere 24:100605. doi: 10.1016/j.rhisph.2022.100605
Li, S., Du, P., Wu, X., He, H., Zhou, L., Dong, F., et al. (2022). Trifluralin Impacts Soil Microbial Community and Functions. Front. Environ. Sci. 10:813871. doi: 10.3389/fenvs.2022.813871
Li, S., Huang, X., Lang, X., Shen, J., Xu, F., and Su, J. (2020a). Cumulative effects of multiple biodiversity attributes and abiotic factors on ecosystem multifunctionality in the Jinsha River valley of southwestern China. For. Ecol. Manag. 472:118281. doi: 10.1016/j.foreco.2020.118281
Li, J., Huang, X., Li, S., Tang, R., and Su, J. (2023). Microbial network complexity and diversity together drive the soil ecosystem multifunctionality of forests during different woodland use intensity in dry and wet season. For. Ecol. Manag. 542
Li, S., Huang, X., Shen, J., Xu, F., and Su, J. (2020b). Effects of plant diversity and soil properties on soil fungal community structure with secondary succession in the Pinus yunnanensis forest. Geoderma 379
Li, F., Li, M., Liu, J., Hu, Y., Zhang, Q., and Zhao, H. (2017). Effect of biochar on fungal abundance of rhizosphere soil and cucumber root growth in greenhouse. Trans. Chin. Soc. Agric. Machinery 48, 265–270–265–341.
Li, B.-B., Li, P.-P., Zhang, W.-T., Ji, J.-Y., Liu, G.-B., and Xu, M.-X. (2021). Deep soil moisture limits the sustainable vegetation restoration in arid and semi-arid Loess Plateau. Geoderma 399. doi: 10.1016/j.geoderma.2021.115122
Li, Z., Yang, Y., Zheng, H., Hu, B., Dai, X., Meng, N., et al. (2023). Environmental changes drive soil microbial community assembly across arid alpine grasslands on the Qinghai-Tibetan Plateau, China. Catena 228
Liang, M., Johnson, D., Burslem, D. F. R. P., Yu, S., Fang, M., Taylor, J. D., et al. (2020). Soil fungal networks maintain local dominance of ectomycorrhizal trees. Nature. Communications 11. doi: 10.1038/s41467-020-16507-y
Liu, J., Jia, X., Yan, W., Zho, Y., and Shangguan, Z. (2020). Changes in soil microbial community structure during long-term secondary succession. Land Degrad. Dev. 31, 1151–1166. doi: 10.1002/ldr.3505
Lu, B., Sun, H., Harris, P., Xu, M., and Charlton, M. (2018). Shp2graph: tools to convert a spatial network into an Igraph Graph in R. ISPRS Int. J. Geo Inf. 7. doi: 10.3390/ijgi7080293
Luo, Y., Du, Z., Yan, Z., Zhao, X., Li, Y., Jiang, H., et al. (2020). Artemisia halodendron litters have strong negative allelopathic effects on earlier successional plants in a semi-arid sandy dune region in China. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00961
Luo, Y., Zhao, X., Li, Y., and Wang, T. (2017). Effects of foliage litter of a pioneer shrub (Artemisia halodendron) on germination from the soil seedbank in a semi-arid sandy grassland in China. J. Plant Res. 130, 1013–1021. doi: 10.1007/s10265-017-0954-0
Ma, W., Wei, F., Zhang, J., Karthe, D., and Opp, C. (2022). Green water appropriation of the cropland ecosystem in China. Sci. Total Environ. 806. doi: 10.1016/j.scitotenv.2021.150597
Maestre, F. T., Quero, J. L., Gotelli, N. J., Escudero, A., Ochoa, V., Delgado-Baquerizo, M., et al. (2012). Plant species richness and ecosystem multifunctionality in global drylands. Science 335, 214–218. doi: 10.1126/science.1215442
Maihemuti, B., Simayi, Z., Alifujiang, Y., Aishan, T., Abliz, A., and Aierken, G. (2021). Development and evaluation of the soil water balance model in an inland arid delta oasis_ Implications for sustainable groundwater resource management. Global Ecol. Conserv. 25:e01408. doi: 10.1016/j.gecco.2020.e01408
Mao, W., Zhang, T., Li, Y., Zhao, X., and Huang, Y. (2012). Allometric response of perennial Pennisetum centrasiaticum Tzvel to nutrient and water limitation in the Horqin Sand Land of China. J. Arid. Land 4, 161–170. doi: 10.3724/SP.J.1227.2012.00161
Martinez, S. (2023). Soil microbial community structure in rice rotation systems with crops and pastures in temperate Uruguay. Appl. Soil Ecol. 188. doi: 10.1016/j.apsoil.2023.104929
Mcculley, R. L., Jobbagy, E. G., Pockman, W. T., and Jackson, R. B. (2004). Nutrient uptake as a contributing explanation for deep rooting in arid and semi-arid ecosystems. Oecologia 141, 620–628. doi: 10.1007/s00442-004-1687-z
Montagna, M., Berruti, A., Bianciotto, V., Cremonesi, P., Giannico, R., Gusmeroli, F., et al. (2018). Differential biodiversity responses between kingdoms (plants, fungi, bacteria and metazoa) along an Alpine succession gradient. Mol. Ecol. 27, 3671–3685. doi: 10.1111/mec.14817
Monteiro, M. R., Marshall, A. J., Lee, C. K., Mcdonald, I. R., and Cary, S. C. (2023). Bringing Antarctica to the lab: a polar desert environmental chamber to study the response of Antarctic microbial communities to climate change. Polar Biol. 46, 445–459. doi: 10.1007/s00300-023-03135-7
Peng, C., Ma, Z., Lei, X., Zhu, Q., Chen, H., Wang, W., et al. (2011). A drought-induced pervasive increase in tree mortality across Canada's boreal forests. Nat. Clim. Chang. 1, 467–471. doi: 10.1038/nclimate1293
Perez-Harguindeguy, N., Diaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234. doi: 10.1071/BT12225
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
Qu, Z.-L., Santalahti, M., Koster, K., Berninger, F., Pumpanen, J., Heinonsalo, J., et al. (2021). Soil fungal community structure in Boreal Pine forests: from Southern to Subarctic Areas of Finland. Front. Microbiol. 12:653896. doi: 10.3389/fmicb.2021.653896
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, eap.2677. doi: 10.1002/eap.2677
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, 1–9. doi: 10.1038/ncomms14196
Shipley, B. (2009). Confirmatory path analysis in a generalized multilevel context. Ecology 90, 363–368. doi: 10.1890/08-1034.1
Stewart, C. E., Roosendaal, D., Denef, K., Pruessner, E., Comas, L. H., Sarath, G., et al. (2017). Seasonal switchgrass ecotype contributions to soil organic carbon, deep soil microbial community composition and rhizodeposit uptake during an extreme drought. Soil Biol. Biochem. 112, 191–203. doi: 10.1016/j.soilbio.2017.04.021
Su, B., Huang, J., Fischer, T., Wang, Y., Kundzewicz, Z. W., Zhai, J., et al. (2018). Drought losses in China might double between the 1.5 degrees C and 2.0 degrees C warming. Proc. Natl. Acad. Sci. U. S. A. 115, 10600–10605. doi: 10.1073/pnas.1802129115
Su, T., Jia, B., Hu, Y., Yang, Q., and Mao, W. (2021). Effects of groundwater depth on soil environmental factors and root biomass of typical plant communities in sandy grassland. Pratacult. Sci. 38, 1694–1705.
Sun, H., Chen, Y., Chen, Y., Zhang, Y., and He, Z. (2020). Groundwater Evapotranspiration in Desert Riparian Forest in the Lower Reaches of the Tarim River. Arid Zone Res. 37, 116–125.
Trinidad Torres-Garcia, M., Oyonarte, C., Cabello, J., Guirado, E., Rodriguez-Lozano, B., and Jacoba Salinas-Bonillo, M. (2022). The potential of groundwater-dependent ecosystems to enhance soil biological activity and soil fertility in drylands. Sci. Total Environ. 826:154111. doi: 10.1016/j.scitotenv.2022.154111
Van Der Heijen, M. G. A. (2008). The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems (vol 11, pg 296, 2008). Ecol. Lett. 11:651.
Van Horn, D. J., Okie, J. G., Buelow, H. N., Gooseff, M. N., Barrett, J. E., and Takacs-Vesbach, C. D. (2014). Soil microbial responses to increased moisture and organic resources along a salinity gradient in a Polar Desert. Appl. Environ. Microbiol. 80, 3034–3043. doi: 10.1128/AEM.03414-13
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, 1–10. 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. Workshop of the SCOPE-Committee-on-Soil-and-Sediment-Biodiversity-and-Ecosystem-Function ing, 2004 Oct 2002 Colorado State Univ, Estes Pk, CO. Washington, DC: Island Press. pp. 1–12.
Wang, T. (2016). Study on the coordinated development of ecosystem and socio-economic system in desertification control: a case study of desertification control in semiarid area in North China. Acta Ecol. Sin. 36, 7045–7048.
Wang, Z., Zhang, Q., Staley, C., Gao, H., Ishii, S., Wei, X., et al. (2019). Impact of long-term grazing exclusion on soil microbial community composition and nutrient availability. Biol. Fertil. Soils 55, 121–134. doi: 10.1007/s00374-018-01336-5
Wu, M.-H., Chen, S.-Y., Chen, J.-W., Xue, K., Chen, S.-L., Wang, X.-M., et al. (2021b). Reduced microbial stability in the active layer is associated with carbon loss under alpine permafrost degradation. Proc. Natl. Acad. Sci. U. S. A. :118:e2025321118. doi: 10.1073/pnas.2025321118
Wu, H., Hao, B., Cai, Y., Liu, G., and Xing, W. (2021a). Effects of submerged vegetation on sediment nitrogen-cycling bacterial communities in Honghu Lake (China). Sci. Total Environ. 755:142541. doi: 10.1016/j.scitotenv.2020.142541
Yang, L., Sun, R., Li, J., Zhai, L., Cui, H., Fan, B., et al. (2023). Combined organic-inorganic fertilization builds higher stability of soil and root microbial networks than exclusive mineral or organic fertilization. Soil Ecol. Lett. 5, 2. doi: 10.1007/s42832-022-0142-6
Ye, C., Gong, Y., Chen, M., Delgado-Baquerizo, M., Che, R., Liu, S., et al. (2023). Revegetation promotes soil microbial network stability in a novel riparian ecosystem. J. Appl. Ecol. doi: 10.1111/1365-2664.14449
Yu, Y., Lee, C., Kim, J., and Hwang, S. (2005). Group-specific primer and probe sets to detect methanogenic communities using quantitative real-time polymerase chain reaction. Biotechnol. Bioeng. 89, 670–679. doi: 10.1002/bit.20347
Zhang, D., Cui, R., Fu, B., Yang, Y., Wang, P., Mao, Y., et al. (2020). Shallow groundwater table fluctuations affect bacterial communities and nitrogen functional genes along the soil profile in a vegetable field. Appl. Soil Ecol. 146:103368. doi: 10.1016/j.apsoil.2019.103368
Zhang, C., Liu, G., Xue, S., and Wang, G. (2016). Soil bacterial community dynamics reflect changes in plant community and soil properties during the secondary succession of abandoned farmland in the Loess Plateau. Soil Biol. Biochem. 97, 40–49. doi: 10.1016/j.soilbio.2016.02.013
Zhang, M., Wang, K., Shi, C., Li, X., Qiu, Z., and Shi, F. (2023). Responses of fungal assembly and co-occurrence network of rhizosphere soil to Amaranthus palmeri invasion in Northern China. J. Fungi 9:5. doi: 10.3390/jof9050509
Zhang, Y., Wu, Z., Singh, V. P., Lin, Q., Ning, S., Zhou, Y., et al. (2023). Agricultural drought characteristics in a typical plain region considering irrigation, crop growth, and water demand impacts. Agric. Water Manag. 282:108266. doi: 10.1016/j.agwat.2023.108266
Zhang, B., Zeng, F., Gao, X., Shareef, M., Zhang, Z., 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
Zhong, Z., Zhang, X., Wang, X., Fu, S., Wu, S., Lu, X., et al. (2020). Soil bacteria and fungi respond differently to plant diversity and plant family composition during the secondary succession of abandoned farmland on the Loess Plateau, China. Plant Soil 448, 183–200. doi: 10.1007/s11104-019-04415-0
Zou, X., Yao, K., Zeng, F., Zhang, C., Zeng, Z., and Zhang, H. (2023). Diversity and Assembly of Bacteria Community in Lime Soil under Different Karst Land-Use Types. Forests 14:672. doi: 10.3390/f14040672
Keywords: semi-arid region, groundwater depth, soil multifunctionality, soil microbial community, relationship
Citation: Zhao S, Zhao X, Li Y, Zhang R, Zhao Y, Fang H and Li W (2023) Impact of altered groundwater depth on soil microbial diversity, network complexity and multifunctionality. Front. Microbiol. 14:1214186. doi: 10.3389/fmicb.2023.1214186
Edited by:
Maria Alicja Niklińska, Jagiellonian University, PolandReviewed by:
Marika Pellegrini, University of L'Aquila, ItalyCésar Marín, Santo Tomás University, Chile
Copyright © 2023 Zhao, Zhao, Li, Zhang, Zhao, Fang and Li. 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: Xueyong Zhao, emhhb3h5QGx6Yi5hYy5jbg==; Yulin Li, bHl1bGluQGdtYWlsLmNvbQ==