Soil organic carbon (SOC) in the active layer (0-2 m) of the Tibetan Plateau (TP) permafrost region is sensitive to climate change, with significant implications for the global carbon cycle. Environmental factors-including parent material, climate, vegetation, topography, soil, and human activities-inevitably drive SOC variations. However, vegetation and climate are likely the two most influential factors impacting SOC variations. To test this hypothesis, we conducted experiments using 31 environmental variables combined with the recursive feature elimination (RFE) algorithm. These experiments showed that RFE retained all vegetation variables [Land cover types (LCT), normalized difference vegetation index (NDVI), leaf area index (LAI), and gross primary productivity (GPP)] as well as two climate variables [Moisture index (MI) and drought index (DI)], supporting our hypothesis. We then analyzed the relationship between SOC and the retained vegetation and climate variables using random forest (RF), Shapley additive explanations (SHAP), and GeoDetector models to quantify the independent and interactive drivers of SOC distribution and to identify the optimal conditions for SOC accumulation. The RF model explained 68% and 42% of the spatial variability in SOC at depths of 0-1 m and 1-2 m, respectively, with SOC stocks higher in the southeast and lower in the northwest. Additionally, SOC stock at 0-1 m was significantly higher (p 0.05). Spearman correlation coefficients results indicated that NDVI, LAI, GPP, and MI had highly significant positive correlations with SOC (p < 0.01), whereas DI had a highly significant negative correlation with SOC (p < 0.01). SHAP analysis revealed environmental thresholds for SOC variations, with notable shifts at NDVI (0.40), LAI (7), GPP (250 g C m(-)(2) year(-)(1)), MI (0.40), and DI (0.50). The spatial distribution of these thresholds aligns with the 400 mm equivalent precipitation line. Additionally, GeoDetector results emphasized that interactions between climate and vegetation factors enhance the explanatory power of individual variables on SOC variations. The swamp meadow type, with an NDVI range of 0.73-0.84, LAI range of 11.06-15.94, and MI range of 0.46-0.56, was identified as the most favorable environment for SOC accumulation. These findings are essential for balancing vegetation and climate conditions to sustain SOC levels and mitigate climate change-driven carbon release.
Background Soil microorganisms in the thawing permafrost play key roles in the maintenance of ecosystem function and regulation of biogeochemical cycles. However, our knowledge of patterns and drivers of permafrost microbial communities is limited in northeastern China. Therefore, we investigated the community structure of soil bacteria in the active, transition and permafrost layers based on 90 soil samples collected from 10 sites across the continuous permafrost region using high-throughput Illumina sequencing. Results Proteobacteria (31.59%), Acidobacteria (18.63%), Bacteroidetes (9.74%), Chloroflexi (7.01%) and Actinobacteria (6.92%) were the predominant phyla of the bacterial community in all soil layers; however, the relative abundances of the dominant bacterial taxa varied with soil depth. The bacterial community alpha-diversity based on the Shannon index and the phylogenetic diversity index both decreased significantly with depth across the transition from active layer to permafrost layer. Nonmetric multidimensional scaling analysis and permutation multivariate analysis of variance revealed that microbial community structures were significantly different among layers. Redundancy analysis and Spearman's correlation analysis showed that soil properties differed between layers such as soil nutrient content, temperature and moisture mainly drove the differentiation of bacterial communities. Conclusions Our results revealed significant differences in bacterial composition and diversity among soil layers. Our findings suggest that the heterogeneous environmental conditions between the three soil horizons had strong influences on microbial niche differentiation and further explained the variability of soil bacterial community structures. This effort to profile the vertical distribution of bacterial communities may enable better evaluations of changes in microbial dynamics in response to permafrost thaw, which would be beneficial to ecological conservation of permafrost ecosystems.
The Qaidam Basin of the Qinghai-Tibet Plateau is a cold, hyper-arid desert that presents extreme challenges to microbial communities. As little is known about variations between surface and subsurface microbial communities, high-throughput DNA sequencing was used in this study to profile bacterial communities of the soil samples collected at different depths in three regions in the Qaidam Basin. The alpha-diversity indices (Chao, Shannon, and Simpson) indicated that bacterial abundance and diversity were higher in the east and the high-elevation regions compared to the west region. In general,Firmicuteswas dominant in the west region, whileProteobacteriaandAcidobacteriawere dominant in the east and the high-elevation regions. The structure of the bacterial communities differed greatly across regions, being strongly correlated with total organic carbon (TOC) and total nitrogen (TN) content. The differences in bacterial communities between the surface and the subsurface soil samples were smaller than the differences across the regions. Network analyses of environmental factors and bacterial genera indicated significant positive correlations in all regions. Overall, our study provides evidence that TOC and TN are the best predictors of both surface and subsurface bacterial communities across the Qaidam Basin. This study concludes that the bacterial community structure is influenced by both the spatial distance and the local environment, but environmental factors are the primary drivers of bacterial spatial patterns in the Qaidam Basin.
Earth's cryosphere and biosphere are extremely sensitive to climate changes, and transitions in states could alter the carbon emission rate to the atmosphere. However, little is known about the climate sensitivities of frozen soil and vegetation production. Moreover, how does climate heterogeneity control the spatial patterns of such sensitivities, and influence regional vulnerability of both frozen soil and vegetation production? Such questions are critical to be answered. We compiled long-time-series dataset including frozen soil depth (FD), normalized difference vegetation index (NDVI), and temperature and precipitation across Tibetan Plateau to quantify their sensitivities. Results reveal large spatial heterogeneity in FD and NDVI sensitivities. Precipitation alleviated FD sensitivities to warming in the cold northeast zone but accelerated FD sensitivities to precipitation in the warm south and southeast. Meanwhile, the positive warming effect on the NDVI was largely offset by slow increase of precipitation. Areas with high FD decreasing rate were found in northeast, inland, and south and southeast zones. Predominate area across the nine eco-regions are characterized as medium FD decreasing rate, and are synchronized with positive NDVI response in inland and west Himalayas, but negative in northeast and south and southeast. Precipitation restriction on NDVI would be pronounced in moist south and southeast. Our study provides new information that makes a much-needed contribution to advancing our understandings of the effects of global climate change on cryosphere and biosphere, which has important implications for global climate and our ability to predict, and therefore prepare for, future global climatic changes. Our attempt confirms that the method we used could be used to identify climate sensitivity of permafrost based on substantial observation data on active layer dynamics in future.
Frozen ground degradation plays an important role in vegetation growth and activity in high-altitude cold regions. This study estimated the spatiotemporal variations in the active layer thickness (ALT) of the permafrost region and the soil freeze depth (SFD) in the seasonally frozen ground region across the Three Rivers Source Region (TRSR) from 1980 to 2014 using the Stefan equation, and differentiated the effects of these variations on alpine vegetation in these two regions. The results showed that the average ALT from 1980 to 2014 increased by 23.01 cm/10a, while the average SFD decreased by 3.41 cm/10a, and both changed intensively in the transitional zone between the seasonally frozen ground and permafrost. From 1982-2014, the increase in the normalized difference vegetation index (NDVI) and the advancement of the start of the vegetation growing season (SOS) in the seasonally frozen ground region (0.0078/10a, 1.83d/10a) were greater than those in the permafrost region (0.0057/10a, 0.39d/10a). The results of the correlation analysis indicated that increases in the ALT and decreases in the SFD in the TRSR could lead to increases in the NDVI and advancement of the SOS. Surface soil moisture played a critical role in vegetation growth in association with the increasing ALT and decreasing SFD. The NDVI for all vegetation types in the TRSR except for alpine vegetation showed an increasing trend that was significantly related to the SFD and ALT. During the study period, the general frozen ground conditions were favorable to vegetation growth, while the average contributions of ALT and SFD to the interannual variation in the NDVI were greater than that of precipitation but less than that of temperature.
Permafrost regions with high soil organic carbon (SOC) storage are extremely vulnerable to global warming. However, our understanding of the temperature sensitivity of SOC decomposition in permafrost regions remains limited, leading to considerable uncertainties in predicting carbon-climate feedback magnitude and direction in these regions. Here, we investigate general patterns and underlying mechanisms of SOC decomposition rate and its temperature sensitivity (Q(10)) at different soil depths across Tibetan permafrost regions. Soils were collected at two depths (0-10 and 20-30 cm) from 91 sites across Tibetan permafrost regions. SOC decomposition rate and Q(10) value were estimated using a continuous-flow incubation system. We found that the SOC decomposition rate in the upper layer (0-10 cm) was significantly greater than that in the lower layer (20-30 cm). The SOC content governed spatial variations in decomposition rates in both soil layers. However, the Q(10) value in the upper layer was significantly lower than that in the lower layer. Soil pH and SOC decomposability had the greatest predictive power for spatial variations in Q(10) value within the upper and lower layers, respectively. Owing to the greater temperature sensitivity in the lower layer, our results imply that subsurface soil carbon is at high risk of loss, and that soil carbon sequestration potential might decrease in these regions in a warming world.