Shuangbo Bi1
Xuerong Shi1
Chunju Peng2
Tianyi Hu1
Jing Chen1Jingchen Xie1Haicheng Li1
Tingting Cao1Man Shi1
Zhikang Wang1
Quan Li1*
Xinzhang Song1- 1State Key Laboratory for Development and Utilization of Forest Food Resources, Zhejiang A&F University, Hangzhou, China
- 2Wenzhou Key Laboratory of Early Sprouting Tea Breeding, Wenzhou Academy of Agricultural Sciences, Wenzhou, China
Introduction: Accurate understanding of soil phosphorus (P) fractions is crucial for enhancing plant productivity and deciphering forest succession patterns; however, the dynamics of rhizosphere soil P fractions and their influencing factors during forest succession or land-type conversion, particularly in highly weathered tropical and subtropical regions, have not been comprehensively elucidated.
Methods: Using a space-for-time replacement strategy, in this study, we examined how Moso bamboo (Phyllostachys edulis) expansion into Chinese fir (Cunninghamia lanceolata) forests affects P fractions in rhizosphere soil across various seasons within a subtropical region. The research focused on seasonal variations in soil P dynamics resulting from this invasive expansion. We further evaluated key drivers, encompassing soil physicochemical characteristics and microbial traits.
Results and discussion: Compared to pure Chinese fir forests, mixed bamboo–fir stands had significant reductions in total P (excluding spring), CaCl2-P, Citrate-P, Enzyme-P (excluding spring), and HCl-P (excluding winter) throughout the seasonal cycle (p < 0.05). Pure bamboo forests showed further reduction in total P, Citrate-P, Enzyme-P, and HCl-P, along with reduced CaCl2-P (except summer and winter) (p < 0.05), with most P fractions (except CaCl2-P in summer, Citrate-P and HCl-P in summer and autumn, and Enzyme-P in summer) being lower in these stands than in mixed forests, which showed a decreasing trend with increasing expansion intensity. CaCl2-P, citrate-P, and HCl-P levels were consistently higher in summer and autumn than in winter and spring across Moso bamboo, Chinese fir, and mixed forest stands. Variations in P fractions were under the major control of nitrogen components and soil pH. This study highlights the importance of clarifying P fraction dynamics to understand forest succession mechanisms and informing P management strategies for enhancing forest productivity.
Introduction
Phosphorus (P) is a critical limiting factor for ecosystem productivity and function, playing fundamental roles in key physiological processes, including anabolism, catabolism, energy transfer, and cellular homeostasis-that collectively underpin plant development and soil fertility (Bai et al., 2024; Jian et al., 2022; Paz-Ares et al., 2022; Yang et al., 2024; Ye et al., 2021). Plants acquire P mainly from the soil, whose P composition is highly complex (Bai et al., 2023; Zhang et al., 2025b). Therefore, accurate characterization of soil P fractions is crucial for improving plant P-use efficiency and understanding P cycling dynamics (Guan et al., 2024; Park et al., 2022). However, traditional approaches, such as the Hedley fractionation method, have emphasized chemical extraction but often overlooked plant acquisition pathways (DeLuca et al., 2015; Hedley and Stewart, 1982). Recently, the bioavailable P fractionation methodology has been crafted, allowing for the sorting of soil P according to how plants absorb it; this includes CaCl2-P, which plants readily take up, Citrate-P, which is inorganic P activated by organic acids secreted by roots, Enzyme-P, which is organic P broken down by phosphatase enzymes, and HCl-P, the pool of P that can be released through proton shedding (DeLuca et al., 2015). Applying this approach to different forest types offers a more ecologically relevant understanding of P speciation and insights into plant P-acquisition strategies.
Tropical and subtropical forests make up roughly 45% of the global woodlands and store approximately 72% of the carbon (C) in woodland biomass, underscoring their essential role in mitigating anthropogenic CO2 emissions and supporting C neutrality (FAO, 2020; Wang et al., 2024). In these regions, the warm and humid climate promotes the immobilization of soil P through fixation with iron and aluminum oxides, resulting in limited P availability for plants (Kochian, 2012). Moso bamboo (Phyllostachys edulis) and Chinese fir (Cunninghamia lanceolata) are important forest resources in the subtropical region of China (Li et al., 2024; Wu et al., 2019; Wang et al., 2025; Zhang et al., 2024, Zhang et al., 2025a). The area of Moso bamboo forests is approximately 5.28 million hectares, with an annual C sequestration capacity of approximately 6.1–7.3 t C·ha-1 (Li et al., 2025; Lv et al., 2025; Song et al., 2017, Song et al., 2020). Chinese fir plantations cover over 9.9 million hectares, with a total standing volume of 755 million cubic meters (Lv et al., 2024). Driven by both natural dispersal and weakened human management (Li et al., 2022; Xu et al., 2020), Moso bamboo, with its clonal growth characteristics (Kotangale et al., 2025; Zhang et al., 2024, 2025a), continues to expand into Chinese fir forests and other forest types (Bai et al., 2016b; Jiang et al., 2024). Therefore, elucidating the changes in soil P fractions during Moso bamboo expansion is crucial not only for understanding the evolution of P cycling in the context of forest succession but also for assessing and enhancing regional forest C sink functions.
The expansion of Moso bamboo significantly alters soil properties, including pH (Ouyang et al., 2022), enzyme activities (Liu et al., 2021), nitrogen (N) mineralization rates (Song et al., 2016a; Wu et al., 2019), C and N cycling (Bai et al., 2016a; Chen et al., 2021; Li et al., 2017), and microbial community structure and function (Liu et al., 2021). Such biological and environmental variable shifts are expected to profoundly influence soil P fractions. Furthermore, the expansion of Moso bamboo has been found to alter the contents of CaCl2-P, HCl-P, and Enzyme-P in bulk soil (Song et al., 2016b; Yang et al., 2024). However, these studies have largely overlooked the rhizosphere soil. Although rhizosphere P pools are generally smaller in size than those in bulk soil, they exhibit higher turnover rates (Qetrani et al., 2024). Consequently, even minor changes in these rhizosphere pools induced by plant expansion may significantly affect plant-available P supply. This highlights the need for in-depth research into rhizosphere P cycling processes in the context of plant expansion. Additionally, Moso bamboo, as a fast-growing clonal plant, possesses a dense, shallow rhizome-root system and an aggressive growth strategy with high nutrient demands (Peng et al., 2021; Wang et al., 2025). In contrast, Chinese fir is a deep-rooted tree species with relatively slow growth and a more conservative nutrient use strategy (Yan et al., 2019). This sharp contrast in growth rate, resource acquisition strategy, and root spatial distribution makes their rhizospheres a highly representative model system for investigating how plant expansion differentially drives soil P cycling processes and underlying microbial mechanisms. Studying this specific interface can provide key insights into the belowground interactions and resource competition between expanding and native tree species. Moreover, given that plants may employ seasonally dependent P-acquisition strategies, temporal variations in rhizospheric P dynamics have not been comprehensively elucidated. Therefore, investigating seasonal variations in rhizosphere soil P fractions amidst Moso bamboo’s expansion of Chinese fir forests is essential for elucidating the P-driven mechanisms underlying subtropical forest succession.
To address this gap, we chose a space-for-time substitution methodology (representing a chronosequence from Chinese fir forest to the bamboo expansion front and mature bamboo stands) to investigate the seasonal variations of soil P fractions and their driving factors. We hypothesized that (1) P fractions would diminish along the expansion gradient, driven by bamboo’s highly efficient P acquisition mechanisms (e.g., organic acid exudation and phosphatase activation), and (2) seasonal variations would dominate the spatiotemporal heterogeneity of P fractions by modulating soil properties and plant demand. Overall, the research offers a crucial understanding of how bamboo expansion affects biogeochemical processes, offering valuable knowledge for crafting strategies to minimize its ecological footprint and promote responsible forest stewardship.
Materials and methods
Study site and experimental design
The research was performed at Banqiao Town (30°10’N, 119°45’E), Lin’an District, China (Figure 1), situated in a subtropical monsoon weather zone, characterized by a mean annual temperature of 17.5 °C and annual precipitation ranging from 1350 to 1500 mm. Soils are slightly acidic red soils (Ultisols) derived from siltstone. Around 2004, pure Chinese fir stands replaced local evergreen broad-leaved forests. The subsequent Moso bamboo–Chinese fir mixed forest developed through bamboo expansion into fir stands. The selected pure Chinese fir forests (PCF) were even-aged stands approximately 15 years old, representing a mature stage prior to significant bamboo invasion. The Moso bamboo-Chinese fir mixed forests (MMC) represented an intermediate invasion stage, with bamboo having established for approximately 15 years. The pure Moso bamboo forests (PMB) represented a late invasion stage, where bamboo has been dominant for over 15 years (old bamboo naturally die or be cut down). This age sequence underpins our space-for-time substitution approach.
Figure 1. Schematic illustrating the experimental design: site locations and basic stand characteristics. DBH, diameter at breast height; PCF, the pure Chinese fir forest; MMC, the mixed Moso bamboo - Chinese fir forest; PMB, pure Moso bamboo forest.
In June 2021, sample transects were established utilizing a space-for-time substitution approach, aligned with the directional extension of Moso bamboo plantations, adopting pure Chinese fir stands and pure Moso bamboo stands as reference boundaries. Three forest types were selected along each transect: (1) PCF, (2) MMC (fir-to-bamboo density ratio: 3:2), and (3) PMB. Three typical sample plots measuring 20 m × 20 m were developed for each forest species under comparable site conditions, resulting in nine sites. A 10-m-wide buffer zone was preserved between neighboring plots. Figure 1 illustrates the configuration of the experimental plots and the features of the stand.
Sampling collection
Samples of the root systems and rhizosphere soil were gathered from Moso bamboo and Chinese fir forest stands during summer (2021.07), autumn (2021.10), winter (2022.01), and spring (2022.04). Within each plot, three Chinese fir in PCF (Mean diameter at breast height (DBH)= 10.5cm) and MMC (Mean DBH = 8.2 cm) and Moso bamboo in MMC (Mean DBH = 7.6 cm) and PMB (Mean DBH = 9.6 cm) were selected based on comparable DBH (Figure 1). After removing surface litter, roots and rhizosphere soil were excavated from four cardinal directions per tree. Rhizosphere soil was collected using the root-shaking method. All specimens were conveyed to the lab in portable coolers, and the root materials were used to determine arbuscular mycorrhizal fungi (AMF) colonization rates. After screening at 2 mm to eliminate pebbles and plant detritus, the soil specimens were partitioned into two groups. A subset was stored frozen at -20 °C for subsequent physicochemical and microbial analyses. The alternative subset was air-dried to evaluate AMF spore density and diverse soil physicochemical characteristics.
Measurement of soil P fractionation
Rhizosphere soil P fractions were determined using a sequential biological-based method (DeLuca et al., 2015), which defines operationally extracted fractions and their plant acquisition pathways (Table 1). Briefly, 0.5 g of fresh soil was individually subjected to treatment with 10 mL of 0.01 mol·L-1 CaCl2 for CaCl2-P, 0.01 mol·L-1 citrate for Citrate-P, 0.02 EU·mL-1 enzyme cocktail for Enzyme-P, and 1 mol·L-1 HCl for HCl-P. After shaking at 180 rpm and 25 °C for 3 h, the mixtures were centrifuged at 10,000 rpm for 1 min. P in the collected supernatants was measured using the malachite green colorimetric technique at 630 nm.
Table 1. Correspondence between the operationally defined phosphorus (P) fractions extracted by the sequential biological-based method and their associated plant acquisition pathways.
Soil physicochemical properties
SWC was assessed using the oven-drying technique. pH was assessed utilizing a pH meter with a soil-to-water ratio of 1:2.5 (w/v) (Peng and Wang, 2016). The soil total N (TN) was measured using a modified Kjeldahl method (Li et al., 2024). Total P (TP) was analyzed by molybdenum blue colorimetry following H2SO4-H2O2 digestion. Soil dissolved organic carbon (DOC) and organic nitrogen (DON) were quantified using a total organic carbon analyzer (Shimadzu TOC-V CPH, Japan). Soil ammonium nitrogen (NH4+-N) and nitrate nitrogen (NO3–-N) were quantified via KCl extraction (soil-to-solution ratio = 1:10), they were determined by the azo-cyanine blue colorimetric method and ultraviolet spectrophotometry, respectively (Lu, 2000).
Soil microbial traits
Microbial biomass C and N (MBC/MBN) were quantified through chloroform fumigation-extraction methodology (Brookes et al., 1985). MBC and MBN were quantified by a TOC analyzer. Soil acid phosphatase activity (ACP) was assayed by p-nitrophenyl phosphate colorimetry (Schneider et al., 2000). Arbuscular mycorrhizal fungi (AMF) colonization rates were determined by acid fuchsin staining and the grid-line intersect method (Kormanik et al., 1980). AMF spore density (SD) was quantified using wet sieving and decanting followed by sucrose density gradient centrifugation (An et al., 1990).
Data analysis
All data underwent homogeneity of variance testing; mathematical transformations were applied when necessary. One-way analyses of variance (ANOVA) with Tukey post-hoc tests assessed differences in the soil P fractions, soil physicochemical properties, microbial biomass, ACP, AMF colonization rate, and spore density among the forest stands. Two-way ANOVA evaluated the effects of forest stands and season on soil P fractions, physicochemical properties, and microbial traits. Pearson correlation analysis assessed the relationships among soil P fractions, physicochemical features, and microbiological traits. Additionally, random forest and linear mixed-effects models were utilized to assess soil and microbial drivers relative importance and contributions to different P fractions using the “rfPermute” package (version 2.5.5) and “lmerTest” package, respectively. PLS-PM assessed the direct and indirect impacts of biotic and abiotic factors on P fractions across forest types by using the ‘plspm’ package (version 4.5.1).
Results
Rhizosphere soil properties
Compared with PCF, MMC and PMB decreased TN (29.0%–59.5%) and NO3–-N (21.7%–95.7%) in all four seasons (p< 0.05, Figure 2). MMC and PMB significantly reduced SWC by 11.5%–31.1% in summer and winter; DON in summer, autumn, and winter; NH4+-N content in spring; and significantly increased soil pH in spring, summer, and autumn (p < 0.05, Figure 2). In spring, soil pH was significantly higher in all three forest stands than in summer and autumn (p < 0.05, Figure 2). Conversely, SWC, DOC, and DON were highest in autumn (p < 0.05, Figure 2). Additionally, TN and ammonium N content were significantly higher in summer than in other seasons, specifically in pure Chinese fir and Moso bamboo forests (p < 0.05, Figure 2). Repeated-measures two-way ANOVA revealed that bamboo expansion, season, and their interaction profoundly influenced soil pH, SWC, TN, TN: TP, NH4+-N, NO3–-N, DON, and DOC (p < 0.05, Table 2).
Figure 2. Influence of Moso bamboo expansion on rhizospheric characteristics of Chinese fir forests across seasons. (a) SWC; (b) pH; (c) TN; (d) NO3--N; (e) NH4+-N; (f) DON; (g) TN:TP; (h) DOC. Uppercase letters denote significant seasonal variations within a forest stand, whereas lowercase letters represent differences between forest stands within a given season. SWC, soil water content; TN, soil total nitrogen; NO3--N, nitrate nitrogen; NH4+-N, ammoniacal nitrogen; DON, dissolved organic nitrogen; TN:TP, the ratio of total nitrogen to total phosphorus; DOC, dissolved organic carbon.
Table 2. A repeated-measures two-way ANOVA was performed to assess the effects of bamboo expansion, season, and their interaction on rhizosphere soil physicochemical properties, microbial characteristics, and P fractions.
Rhizosphere microbial traits
The rhizosphere soil of both the mixed forests and pure Moso bamboo forests exhibited significantly low AMF spore density (45.1%–70.1%) during the summer and autumn seasons compared to pure Chinese fir forests (p< 0.05, Figure 3). However, the colonization rate of AMF in the Moso bamboo plots decreased compared to the Chinese fir forests (p< 0.05, Figure 3). The MBC in pure Moso bamboo forests was significantly lower than that in pure Chinese fir forests during spring, while no significant change was observed in mixed forests. The MBN in pure Moso bamboo forests was significantly higher than that in pure Chinese fir forests during summer and autumn but showed a significant decrease in spring (p< 0.05, Figure 3). There were no significant differences in ACP activity among the other tree species (p > 0.05, Figure 3). AMF spore density peaked in summer across all stands. In contrast, MBC and MBN were highest in autumn (Figure 3). The repeated-measures two-way ANOVA revealed that bamboo expansion, season, and their interaction significantly affected MBN and AMF spore density and colonization rate (p < 0.05, Table 2).
Figure 3. Influence of Moso bamboo expansion on rhizospheric microbial characteristics of Chinese fir forests across seasons. (a) AMF colonization rate; (b) SD; (c) ACP; (d) MBC; (e) MBN; (f) MBC:MBN. Capital letters denote seasonal marked variations within a forest stand; lowercase letters signify notable variances between forest stands within a season. SD, AMF spore density; ACP, acid phosphatase activity; MBC, microbial biomass carbon; MBN, microbial biomass nitrogen.
Rhizosphere P fractions
Compared to the PCF, both MMC and PMB significantly reduced soil TP content in summer and winter, with the lowest values observed in PMB (p < 0.05; Figure 4). Across all three forest stands, HCl-P constituted the largest proportion of bioavailable P, followed by Citrate-P, while CaCl2-P and Enzyme-P accounted for relatively minor proportions—a pattern that remained consistent across seasons (Figure 4).
Figure 4. Influence of Moso bamboo development into Chinese fir forests on (a) TP; individual phosphorus fractions (b–e), ratios of other P fractions to CaCl2-P (f–h), and the proportional contributions of these fractions to total phosphorus (i–l). Significant variations between forest stands within a season are indicated by lowercase letters, while significant differences between seasons within a forest stand are indicated by uppercase letters. TP, soil total phosphorus; MCF, Mixed Chinese fir forest; MMB, Mixed Moso bamboo forest.
During spring and autumn, MMC significantly decreased CaCl2-P by 13.6%–49.3%, while Citrate-P (24.9%–53.1%), HCl-P (36% - 62.1%) contents were significantly decreased in summer, autumn, and spring relative to the PCF. PMB significantly reduced Enzyme-P (51.5%–77.5%) and HCl-P (30.8%–70.9%) across all four seasons, significantly decreased Citrate-P (36.1% - 53.5%) during spring, summer, and autumn, while only significantly lowered CaCl2-P in spring (75.0%) and autumn (89.8%). (p < 0.05; Figure 4). Compared with PCF, PMB significantly reduced the ratios of Citrate−P to CaCl2−P, Enzyme−P to CaCl2−P, and HCl−P to CaCl2−P in summer, whereas these ratios were significantly increased in autumn and spring (p < 0.05; Figure 4).
CaCl2-P, Citrate-P, and HCl-P were markedly elevated during summer–autumn than winter–spring across the PCF and PMF stands, while Enzyme-P showed an inverse seasonal pattern with peak values in winter–spring across the PCF and PMB stands (p < 0.05; Figure 4). In terms of seasonal dynamics, CaCl2-P, Citrate-P, and HCl-P were markedly elevated during summer–autumn compared to winter–spring across both the PCF and PMB stands. In contrast, Enzyme-P exhibited an inverse pattern, with peak values occurring in winter–spring across the PCF and PMB stands (p < 0.05; Figure 4). The repeated-measures two-way ANOVA revealed that bamboo expansion, season, and their interaction significantly influenced all P fractions (p < 0.05, Table 2).
Drivers of soil P fraction dynamics
All soil P fractions exhibited substantial positive associations with TN, NO3–-N, and DON (p < 0.05; Figure 5). Moreover, CaCl2-P, Citrate-P, and HCl-P exhibited an advantageous connection with AMF spore density while demonstrating a negative correlation with soil pH (p < 0.05; Figure 5).
Figure 5. Correlations between rhizosphere soil phosphorus fractions, soil physicochemical properties, and microbial traits across different seasons and forest stands. *p < 0.005, **p < 0.001; ***p < 0.0001. SWC, soil water content; TN, soil total nitrogen; TP, soil total phosphorus; NO3--N, nitrate nitrogen; NH4+-N, ammoniacal nitrogen; TN:TP, the ratio of total nitrogen to total phosphorus; DON, dissolved organic nitrogen; DOC, dissolved organic carbon; SD, AMF spore density; ACP, acid phosphatase activity; MBC, microbial biomass carbon; MBN, microbial biomass nitrogen.
Hierarchical partitioning within the mixed-effects model revealed that soil physicochemical properties accounted for substantially greater variance (62.7%–85.9%) in all four P fractions compared to microbial factors (14.1%–37.3%), highlighting the predominant role of abiotic drivers over biological processes in regulating P speciation (Figure 6). The random forest analysis further identified NO3–-N and pH as the most important factors influencing CaCl2-P, while AMF spore density and pH were the dominant predictors for Citrate-P. TP and TN exerted the strongest effects on Enzyme-P. In contrast, TN, TP, AMF spore density, and pH were the primary factors governing HCl-P dynamics (Figure 6). PLS-PM indicated that both abiotic (pH, SWC, and N forms) and biotic factors (microbial biomass andAMF spore density) acted synergistically to regulate soil P fractions.Among these, soil N components (TN, NO3–-N, TN: TP, DON, and NH4+-N) served as the strongest positive direct drivers, while Moso bamboo expansion exhibited consistent negative direct effects (Figure 7).
Figure 6. Analyses applying random forest and mixed-effects models to assess soil physicochemical and microbial influences on rhizosphere phosphorus fractions: (a) CaCl2-P, (b) Citrate-P, (c) Enzyme-P, and (d) HCl-P. *p < 0.005. TP, soil total phosphorus; TN, soil total nitrogen; NH4+-N, ammoniacal nitrogen; DON, dissolved organic nitrogen; DOC, dissolved organic carbon; MBC, microbial biomass carbon; NO3--N, nitrate nitrogen; SD, AMF spore density; SWC, soil water content; MBN, microbial biomass nitrogen; ACP, acid phosphatase activity.
Figure 7. PLS-PM analysis revealing the direct and indirect impacts of biotic and abiotic factors on rhizosphere soil phosphorus fractions during Moso bamboo forests’ expansion into Chinese fir forests. In this diagram, red and blue lines denote positive and negative correlations, respectively, while solid and dashed lines indicate significant and non-significant pathways. Adjacent to each arrow, numerical values display standardized path coefficients, and the arrows’ thickness corresponds to the association’s magnitude *p < 0.05, **p < 0.01, ***p < 0.001. SWC, soil water content; TN, soil total nitrogen; NO3--N, nitrate nitrogen; DON, dissolved organic nitrogen.
Discussion
Reduction of rhizosphere P pools accompanying bamboo expansion
This study provides the first direct evidence that Moso bamboo expansion into Chinese fir forests markedly reduces total P and all measured bioavailable P fractions in rhizosphere soil, fully supporting the first hypothesis. Consistent with this, Wu et al. (2018) documented a reduction in total P in bulk soil following bamboo growth within coniferous forests. In contrast, expansion into broad-leaved forests generally does not affect total P levels and yields varying outcomes in bioavailable P (Yang et al., 2024). These discrepancies underscore the pivotal role of forest types and soil substrates role in governing P dynamics. Consequently, systematic comparisons across forest and soil types are urgently needed to clarify the generality of bamboo-driven P reduction.
The expansion of Moso bamboo forests caused a substantial decline in soil CaCl2-P content (Figure 4), presumably owing to heightened nutritional requirements from both bamboo and fir in mixed stands (Wu et al., 2019). Additionally, the elevated biomass of Moso bamboo compared to other subtropical tree types augments its absorption of labile P (Song et al., 2017; Xu et al., 2020), further reducing rhizosphere available P. However, in summer, the CaCl2-P content in the rhizosphere of pure Moso bamboo forests was significantly higher than that in pure Chinese fir forests and Moso bamboo–Chinese fir mixed forests. This may be attributed to higher rhizospheric microbial activity in pure bamboo stands during the summer, which enhances the mobilization of this P fraction. Therefore, we acknowledge that this conclusion is based on a limited set of measured microbial attributes (MBC, MBN, ACP, AMF). A more comprehensive profiling of the microbial community (e.g., via high-throughput sequencing) and its functional potential might reveal a greater contribution of microbial processes to P cycling than currently estimated. Therefore, the role of microbial factors in our system may be underrepresented.
The reduction in Citrate-P and HCl-P may be attributed to Moso bamboo’s extensive rhizome–root system; substantial belowground biomass improves P-acquisition efficiency through fine roots (Song et al., 2017, Song et al., 2020; Ye et al., 2021), thereby diminishing investment in exudate-based P-mobilizing strategies and resulting in the decline of these P fractions. The observed rise in soil pH following bamboo expansion supports this interpretation (Figures 2, 5). This was because soil pH impacts the composition and variety of P-solubilizing microbes and affects iron and aluminum oxides, which in turn influence soil P fractions (Hou et al., 2018b; Wu et al., 2020; Wan et al., 2021; Zhou et al., 2020). However, some studies have found that Moso bamboo expansion can lead to a decrease in soil pH (Yang et al., 2024). The divergence in these findings may be attributed to variations in geographical environments or differences in the conditions of the invaded forest stands. The pH increase in our study may be related to specific rhizosphere processes of Moso bamboo, such as distinct exudate inputs, alterations in microbial community structure, or differences in organic acid metabolism pathways, the precise mechanisms of which require further investigation. Nevertheless, the observed pH increase was consistent with the trends in Citrate-P and HCl-P reported in our study, jointly supporting the core hypothesis that bamboo expansion influences phosphorus availability by modifying the rhizosphere microenvironment.
Decline in Enzyme-P may be attributed to the higher C:P ratio of Moso bamboo litter compared to Chinese fir litter (Song et al., 2015), which stimulates microbial P immobilization. Additionally, rich in recalcitrant lignin, bamboo litter hinders the breakdown of organic substances, reducing enzymatically accessible P renewal (Li et al., 2019; Liu et al., 2024). Diminished ACP activity and microbial populations across pure bamboo and mixed soils, versus pure Chinese fir, partially reinforce this process (Figure 3). In summary, the reduction in both total and bioavailable P during Moso bamboo expansion indicates a shift toward a “high turnover, low storage” P-cycling strategy, aligning with earlier observations (Li et al., 2021).
In all three forest stands, the high proportions of HCl-P and Citrate-P in TP, along with their ratios to CaCl2-P, exceeded 1 (Figure 4), suggesting that both Moso bamboo and Chinese fir primarily rely on root-exudate-mediated P activation for P acquisition—a conclusion aligned with Yang et al. (2024). Furthermore, the significantly lower AMF colonization rate in Moso bamboo than in Chinese fir indicates that the latter relies more heavily on the mycorrhizal pathway for P acquisition. However, as these results only indirectly reflect plant P uptake pathways, more direct methodologies are needed to clarify P acquisition strategies and better understand the P-cycling dynamics during bamboo forest expansion. Furthermore, if some mixed stands were not formed by bamboo invasion, the conclusion that “bamboo expansion leads to a decrease in rhizosphere phosphorus content” requires further verification, and the underlying mechanisms warrant further investigation.
Moso bamboo expansion drives persistent rhizospheric P decline despite seasonal dynamics
CaCl2-P, Citrate-P, and HCl-P levels were notably elevated in summer and fall compared to winter and spring across Moso bamboo, Chinese fir, and mixed forest stands (Figure 4), corroborating the second hypothesis. This seasonal pattern is driven by two interconnected mechanisms: (1) Optimal summer–autumn climatic conditions and elevated microbial activity, accelerating organic matter decomposition and subsequent release of labile P (e.g., CaCl2-P). In this study, the higher MBC content in autumn, which showed a significantly positive correlation with CaCl2–P (Figures 3, 4), supports this mechanism. (2) During the growing season (summer–autumn), root organic acid exudation rates are higher than they are in winter–spring, and high humidity reduced soil pH, thereby increasing the abundance of Fe(III)-reducing bacteria (e.g., Geobacter) (Hou et al., 2018a); this enhances the dissolution of Citrate-P and HCl-P. Lower pH observed in summer and autumn compared to winter and spring further supports this mechanism (Figure 2).
In all three forest types, the rhizosphere Enzyme-P concentration was markedly decreased in summer compared to other seasons (Figure 4). This pattern can be attributed to two key mechanisms: (1) concentrated litterfall in autumn and winter increases organic P input, whereas summer—being the peak growing season—exhibits reduced organic P input, leading to lower Enzyme-P; and (2) lower temperatures suppress acid phosphatase activity, thereby limiting the mineralization of Enzyme-P. This study also found higher ACP activity in summer compared to other seasons (Figure 3), further supporting the second mechanism. Despite distinct seasonal patterns in rhizosphere P fractions across the forest types, the overall trend of expansion-induced P reduction from bamboo to fir stands remained unaltered.
However, we acknowledge that the seasonal effects discussed here also incorporate between-year variability and potential influences from management activities (e.g., bamboo shoot harvesting typically occurring in winter–spring) (Liu et al., 2023), which may confound the interpretation of purely seasonal drivers.
Conclusion
In this study, a space-for-time replacement method was used to demonstrate, for the first time, the dynamic evolution of rhizosphere P pools throughout the spread of Moso bamboo inside Chinese fir woodlands (Figure 8). Total P and bioavailable fractions demonstrated a significant decrease along the expansion gradient. This pattern remained stable across seasons. Reduction was directly linked to bamboo’s high root biomass and fine-root foraging strategy for efficient P acquisition. Rhizosphere Citrate-P and HCl-P peaked in summer–autumn while Enzyme-P peaked in winter–spring, demonstrating synchronous P-fraction seasonality across forest types. These findings advance niche competition theory by elucidating P-form turnover during bamboo expansion and guide seasonal P management for subtropical plantations, supporting sustainable forestry under “dual carbon” goals. To resolve critical knowledge gaps, subsequent research must dissect interspecific variation in root-foraging precision and leaf-P resorption efficiency between Moso bamboo and Chinese fir, thereby refining P cycling models for subtropical forests.
Figure 8. A conceptual graphic depicting the alterations in soil phosphorus fractions and their principal driving forces after the expansion of Moso bamboo into Chinese fir woods. TP, soil total phosphorus; SWC, soil water content; TN, soil total nitrogen; NO3--N, nitrate nitrogen; ACP, acid phosphatase activity; SD, AMF spore density.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
SBB: Data curation, Formal Analysis, Software, Writing – original draft. XRS: Data curation, Formal Analysis, Investigation, Software, Visualization, Writing – original draft. CJP: Conceptualization, Methodology, Software, Writing – review & editing. TYH: Investigation, Formal Analysis, Writing – review & editing. JC: Investigation, Writing – review & editing. JCX: Investigation, Writing – review & editing. HCL: Investigation, Writing – review & editing. TTC: Writing – review & editing. MS: Writing -review & editing. ZKW: Writing – review & editing. QL: Conceptualization, Project administration, Funding acquisition, Project administration, Writing – review & editing. XZS: Conceptualization, Project administration, Resources, Supervision, Methodology, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication.This work was financially supported by the Zhejiang Provincial National Natural Science Foundation of China (LQ23C160006, LQ24C160005) and National Natural Science Foundation of China (32401344, 32125027, 32301674).
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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
An, Z.-Q., Hendrix, J. W., Hershman, D. E., and Henson, G. T. (1990). Evaluation of the “most probable number” (mpn) and wet-sieving methods for determining soil-borne populations of endogonaceous mycorrhizal fungi. Mycologia 82, 576–581. doi: 10.1080/00275514.1990.12025932
Bai, J., Chen, R., Men, X., and Cheng, X. (2023). Divergent linkages of soil phosphorus fractions to edaphic properties following afforestation in the riparian zone of the upper yangtze river, China. Chemosphere 313, 137452. doi: 10.1016/j.chemosphere.2022.137452
Bai, S., Conant, R. T., Zhou, G., Wang, Y., Wang, N., Li, Y., et al. (2016a). Effects of moso bamboo encroachment into native, broad-leaved forests on soil carbon and nitrogen pools. Sci. Rep. 6, 31480. doi: 10.1038/srep31480
Bai, S., Wang, Y., Conant, R. T., Zhou, G., Xu, Y., Wang, N., et al. (2016b). Can native clonal moso bamboo encroach on adjacent natural forest without human intervention? Sci. Rep. 6, 31504. doi: 10.1038/srep31504
Bai, K., Wang, W., Zhang, J., Yao, P., Cai, C., Xie, Z., et al. (2024). Effects of phosphorus-solubilizing bacteria and biochar application on phosphorus availability and tomato growth under phosphorus stress. BMC Biol. 22, 211. doi: 10.1186/s12915-024-02011-y
Brookes, P. C., Landman, A., Pruden, G., and Jenkinson, D. S. (1985). Chloroform fumigation and the release of soil nitrogen: a rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol. Biochem. 17, 837–842. doi: 10.1016/0038-0717(85)90144-0
Chen, Z., Li, Y., Chang, S. X., Xu, Q., Li, Y., Ma, Z., et al. (2021). Linking enhanced soil nitrogen mineralization to increased fungal decomposition capacity with moso bamboo invasion of broadleaf forests. Sci. Tot. Environ. 771, 144779. doi: 10.1016/j.scitotenv.2020.144779
DeLuca, T. H., Glanville, H. C., Harris, M., Emmett, B. A., Pingree, M. R. A., De Sosa, L. L., et al. (2015). A novel biologically-based approach to evaluating soil phosphorus availability across complex landscapes. Soil Biol. Biochem. 88, 110–119. doi: 10.1016/j.soilbio.2015.05.016
FAO (2020). Global forest resources assessment 2020: Main report. Rome: Food and Agriculture Organization of the United Nations. doi: 10.4060/ca8753en
Guan, X., Chen, J., Liu, G., and Wang, X. (2024). Soil phosphorus forms in saline soil after the application of biomass materials. Agronomy 14, 255. doi: 10.3390/agronomy14020255
Hedley, M. J. and Stewart, J. W. B. (1982). Method to measure microbial phosphate in soils. Soil Biol. Biochem. 14, 377–385. doi: 10.1016/0038-0717(82)90009-8
Hou, E., Tan, X., Heenan, M., and Wen, D. (2018a). A global dataset of plant available and unavailable phosphorus in natural soils derived by hedley method. Sci. Data 5, 180166. doi: 10.1038/sdata.2018.166
Hou, E., Wen, D., Kuang, Y., Cong, J., Chen, C., He, X., et al. (2018b). Soil pH predominantly controls the forms of organic phosphorus in topsoils under natural broadleaved forests along a 2500 km latitudinal gradient. Geoderma 315, 65–74. doi: 10.1016/j.geoderma.2017.11.041
Jian, Z., Ni, Y., Lei, L., Xu, J., Xiao, W., and Zeng, L. (2022). Phosphorus is the key soil indicator controlling productivity in planted masson pine forests across subtropical China. Sci. Tot. Environ. 822, 153525. doi: 10.1016/j.scitotenv.2022.153525
Jiang, R., Lin, J., Zhang, X., and Kang, M. (2024). Investigating changes of forest aboveground biomass induced by moso bamboo expansion with terrestrial laser scanner. Ecol. Inf. 83, 102812. doi: 10.1016/j.ecoinf.2024.102812
Kormanik, P. P., Bryan, W. C., and Schultz, R. C. (1980). Procedures and equipment for staining large numbers of plant root samples for endomycorrhizal assay. Can. J. Microbiol. 26, 536–538. doi: 10.1139/m80-090
Kotangale, P., Agashe, A., Sawarkar, R., Dewangan, C., Tijare, G., and Singh, L. (2025). Unlocking the hidden power of bamboo rhizomes: a comprehensive review of their role in nutrient storage, water retention, and plant growth. Adv. Bamb. Sci. 10, 100122. doi: 10.1016/j.bamboo.2025.100122
Li, Y., Li, Y., Chang, S. X., Xu, Q., Guo, Z., Gao, Q., et al. (2017). Bamboo invasion of broadleaf forests altered soil fungal community closely linked to changes in soil organic C chemical composition and mineral N production. Plant Soil 418, 507–521. doi: 10.1007/s11104-017-3313-y
Li, F.-R., Liu, L.-L., Liu, J.-L., and Yang, K. (2019). Abiotic and biotic controls on dynamics of labile phosphorus fractions in calcareous soils under agricultural cultivation. Sci. Tot. Environ. 681, 163–174. doi: 10.1016/j.scitotenv.2019.05.091
Li, Q., Lv, J., Peng, C., Xiang, W., Xiao, W., and Song, X. (2021). Nitrogen -addition accelerates phosphorus cycling and changes phosphorus use strategy in a subtropical moso bamboo forest. Environ. Res. Lett. 16, 24023. doi: 10.1088/1748-9326/abd5e1
Li, Y., Wang, N., and Latiff, A. R. A. (2025). Development of the bamboo forest economy: reviewing China’s “bamboo as a substitute for plastic initiative” and its development. Adv. Bamb. Sci. 11, 100130. doi: 10.1016/j.bamboo.2025.100130
Li, Q., Zhang, C., Shi, M., Lv, J., Peng, C., Zhang, J., et al. (2024). Long-term nitrogen addition has a positive legacy effect on soil respiration in subtropical moso bamboo forests. Geoderma 452, 117092. doi: 10.1016/j.geoderma.2024.117092
Li, C., Zhong, Q., Yu, K., and Li, B. (2022). Carbon, nitrogen, and phosphorus stoichiometry between leaf and soil exhibit the different expansion stages of moso bamboo (phyllostachys edulis (carriere) J. Houzeau) into chinese fir (cunninghamia lanceolata (lamb.) hook.) forest. Forests 13, 1830. doi: 10.3390/f13111830
Liu, J., Compson, Z. G., Gui, X., Yang, Q., Song, Q., Huang, D., et al. (2024). Weak responses of soil microorganisms to leaf litter inputs after native phyllostachys edulis invasion into adjacent native forests. Plant Soil 494, 685–699. doi: 10.1007/s11104-023-06311-0
Liu, C., Zheng, C., Wang, L., Zhang, J., Wang, Q., Shao, S., et al. (2023). Moso bamboo invasion changes the assembly process and interactive relationship of soil microbial communities in a subtropical broadleaf forest. For. Ecol. Manag. 536, 120901. doi: 10.1016/j.foreco.2023.120901
Liu, C., Zhou, Y., Qin, H., Liang, C., Shao, S., Fuhrmann, J. J., et al. (2021). Moso bamboo invasion has contrasting effects on soil bacterial and fungal abundances, co-occurrence networks and their associations with enzyme activities in three broadleaved forests across subtropical China. For. Ecol. Manag. 498, 119549. doi: 10.1016/j.foreco.2021.119549
Lu, R. K. (2000). Soil and Agro-Chemistry Analytical Methods (Beijing, China: China Agricultural Sciences Press).
Lv, Z., Duan, A., and Zhang, J. (2024). Influence of forest age, tree size, and climate factors on biomass and carbon storage allocation in chinese fir forests. Ecol. Indic. 163, 112096. doi: 10.1016/j.ecolind.2024.112096
Lv, J., Li, Q., Cao, T., Shi, M., Peng, C., Deng, L., et al. (2025). A compartmentation approach to deconstruct ecosystem carbon fluxes of a moso bamboo forest in subtropical China. For. Ecosyst. 12, 100286. doi: 10.1016/j.fecs.2024.100286
Ouyang, M., Tian, D., Pan, J., Chen, G., Su, H., Yan, Z., et al. (2022). Moso bamboo (phyllostachys edulis) invasion increases forest soil pH in subtropical China. Catena 215, 106339. doi: 10.1016/j.catena.2022.106339
Park, Y., Solhtalab, M., Thongsomboon, W., and Aristilde, L. (2022). Strategies of organic phosphorus recycling by soil bacteria: acquisition, metabolism, and regulation. Environ. Microbiol. Rep. 14, 3–24. doi: 10.1111/1758-2229.13040
Paz-Ares, J., Puga, M. I., Rojas-Triana, M., Martinez-Hevia, I., Diaz, S., Poza-Carrión, C., et al. (2022). Plant adaptation to low phosphorus availability: core signaling, crosstalks, and applied implications. Mol. Plant 15, 104–124. doi: 10.1016/j.molp.2021.12.005
Peng, C., Song, Y., Li, C., Mei, T., Wu, Z., Shi, Y., et al. (2021). Growing in mixed stands increased leaf photosynthesis and physiological stress resistance in moso bamboo and mature chinese fir plantations. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.649204
Peng, X. and Wang, W. (2016). Stoichiometry of soil extracellular enzyme activity along a climatic transect in temperate grasslands of northern China. Soil Biol. Biochem. 98, 74–84. doi: 10.1016/j.soilbio.2016.04.008
Qetrani, S., Bouray, M., and Oukarroum, A. (2024). Phosphorus mobilization and acquisition in the alkaline-calcareous rhizosphere: a synthesis. Rhizosphere 30, 100907. doi: 10.1016/j.rhisph.2024.100907
Schneider, K., Turrión, M., and Gallardo, J. (2000). Modified method for measuring acid phosphatase activities in forest soils with high organic matter content. Commun. Soil Sci. Plant Anal. 31, 3077–3088. doi: 10.1080/00103620009370651
Song, X., Chen, X., Zhou, G., Jiang, H., and Peng, C. (2017). Observed high and persistent carbon uptake by moso bamboo forests and its response to environmental drivers. Agric. For. Meteorol. 247, 467–475. doi: 10.1016/j.agrformet.2017.09.001
Song, X., Gu, H., Wang, M., Zhou, G., and Li, Q. (2016b). Management practices regulate the response of moso bamboo foliar stoichiometry to nitrogen deposition. Sci. Rep. 6, 24107. doi: 10.1038/srep24107
Song, Q., Ouyang, M., Yang, Q., Lu, H., Yang, G., Chen, F., et al. (2016a). Degradation of litter quality and decline of soil nitrogen mineralization after moso bamboo (phyllostachys pubscens) expansion to neighboring broadleaved forest in subtropical China. Plant Soil 404, 113–124. doi: 10.1007/s11104-016-2835-z
Song, X., Peng, C., Ciais, P., Li, Q., Xiang, W., Xiao, W., et al. (2020). Nitrogen addition increased CO2 uptake more than non-CO2 greenhouse gases emissions in a moso bamboo forest. Sci. Adv. 6, eaaw5790. doi: 10.1126/sciadv.aaw5790
Song, X., Zhou, G., Gu, H., and Qi, L. (2015). Management practices amplify the effects of N deposition on leaf litter decomposition of the moso bamboo forest. Plant Soil 395, 391–400. doi: 10.1007/s11104-015-2578-2
Wan, W., Hao, X., Xing, Y., Liu, S., Zhang, X., Li, X., et al. (2021). Spatial differences in soil microbial diversity caused by pH -driven organic phosphorus mineralization. Land. Degrad. Dev. 32, 766–776. doi: 10.1002/ldr.3734
Wang, Z., Li, Q., Shi, M., Leite, M. F. A., Chen, X., Kuramae, E. E., et al. (2025). Compartmentalized homeostasis drives high bamboo forest productivity under nutrient imbalance. Adv. Sci., e17442. doi: 10.1002/advs.202517442
Wang, Z., Tian, H., Pan, S., Shi, H., Yang, J., Liang, N., et al. (2024). Phosphorus limitation on CO2 fertilization effect in tropical forests informed by a coupled biogeochemical model. For. Ecosyst. 11, 100210. doi: 10.1016/j.fecs.2024.100210
Wu, C., Mo, Q., Wang, H., Zhang, Z., Huang, G., Ye, Q., et al. (2018). Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) invasion affects soil phosphorus dynamics in adjacent coniferous forests in subtropical China. Ann. For. Sci. 75, 24. doi: 10.1007/s13595-018-0703-0
Wu, H., Xiang, W., Chen, L., Ouyang, S., Xiao, W., Li, S., et al. (2020). Soil phosphorus bioavailability and recycling increased with stand age in chinese fir plantations. Ecosystems 23, 973–988. doi: 10.1007/s10021-019-00450-1
Wu, H., Xiang, W., Ouyang, S., Forrester, D. I., Zhou, B., Chen, L., et al. (2019). Linkage between tree species richness and soil microbial diversity improves phosphorus bioavailability. Funct. Ecol. 33, 1549–1560. doi: 10.1111/1365-2435.13355
Xu, Q.-F., Liang, C.-F., Chen, J.-H., Li, Y.-C., Qin, H., and Fuhrmann, J. J. (2020). Rapid bamboo invasion (expansion) and its effects on biodiversity and soil processes. Global Ecol. Conserv. 21, e00787. doi: 10.1016/j.gecco.2019.e00787
Yan, X.-L., Wang, C., Ma, X., and Wu, P. (2019). Root morphology and seedling growth of three tree species in southern China in response to homogeneous and heterogeneous phosphorus supplies. Trees 33, 1283–1297. doi: 10.1007/s00468-019-01858-x
Yang, D., Shi, F., Fang, X., Zhang, R., Shi, J., and Zhang, Y. (2024). Effect of the moso bamboo pyllostachys edulis (carrière) J.houz. on soil phosphorus bioavailability in a broadleaf forest (jiangxi province, China). Forests 15, 328. doi: 10.3390/f15020328
Ye, J., Yue, C., Hu, Y., and Ma, H. (2021). Spatial patterns of global-scale forest root-shoot ratio and their controlling factors. Sci. Tot. Environ. 800, 149251. doi: 10.1016/j.scitotenv.2021.149251
Zhang, X., Huang, Z., Zhong, Z., Li, Q., and Bian, F. (2025b). Forest management impacts on soil phosphorus cycling: insights from metagenomics in moso bamboo plantations. J. Environ. Manag. 373, 123735. doi: 10.1016/j.jenvman.2024.123735
Zhang, W., Shi, M., Yang, K., Zhang, J., Gao, Z., El-Kassaby, Y. A., et al. (2024). Regulatory networks of senescence-associated gene-transcription factors promote degradation in moso bamboo shoots. Plant Cell Environ. 47, 3654–3667. doi: 10.1111/pce.14950
Zhang, J., Shi, M., Zhu, C., Yang, K., Li, Q., Song, X., et al. (2025a). Stable isotope labelling and gene expression analysis reveal dynamic nitrogen-supply mechanisms for rapid growth of moso bamboo. Hortic. Res. 12, uhaf062. doi: 10.1093/hr/uhaf062
Keywords: forest succession, phosphorus fractions, rhizosphere soil, soil pH, subtropical forest
Citation: Bi S, Shi X, Peng C, Hu T, Chen J, Xie J, Li H, Cao T, Shi M, Wang Z, Li Q and Song X (2026) Expansion of Moso bamboo into Chinese fir stands persistently depletes rhizosphere bioavailable P pools: a seasonal, space-for-time approach. Front. Plant Sci. 16:1718574. doi: 10.3389/fpls.2025.1718574
Received: 04 October 2025; Accepted: 18 December 2025; Revised: 14 December 2025;
Published: 20 January 2026.
Edited by:
Ming Dong, Hangzhou Normal University, ChinaCopyright © 2026 Bi, Shi, Peng, Hu, Chen, Xie, Li, Cao, Shi, Wang, Li and Song. 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: Quan Li, cXVhbmxpQHphZnUuZWR1LmNu