Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 19 November 2025

Sec. Functional Plant Ecology

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1705626

This article is part of the Research TopicUncovering plant adaptation mechanisms for effective ecological restorationView all 3 articles

Soil quality enhancement drives tree growth and broadleaf dominance in fir-broadleaf mixed plantations

Zhigao Fu,Zhigao Fu1,2Yihua Xiao,*Yihua Xiao1,2*Shirong Liu,Shirong Liu3,4Han XuHan Xu1Yan Wang,Yan Wang1,2Huosheng ZhuHuosheng Zhu5
  • 1Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou, China
  • 2Pearl River Delta Farmland Shelterbelt Ecosystem Research Station, National Forestry and Grassland Administration, Guangzhou, China
  • 3Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, China
  • 4Baotianman Forest Ecosystem Research Station, National Forestry and Grassland Administration, Nanyang, Henan, China
  • 5Lechang Forest Farm of Guangdong Province, Forestry Administration of Guangdong Province, Lechang, Guangdong,  China

Introduction: Evaluating soil quality is essential for guiding reforestation and land management strategies, particularly in degraded Chinese fir plantations where long-term productivity and successional dynamics remain poorly understood.

Methods: This study assessed ten mixed-species planting patterns to quantify the Soil Quality Index (SQI) using a Minimum Data Set (MDS) approach, which reduces data redundancy by statistically identifying key indicators from a larger dataset, thereby effectively capturing essential soil functions, and subsequently explored the relationships between SQI and stand growth, structural diversity, biomass, net primary productivity (NPP), as well as percentage of broadleaf species (PBS).

Results: Significant differences were observed across planting patterns in diameter at breast height (DBH), tree height (TH), stand biomass (FB), structural diversity (variation in DBH [CVD] and Gini coefficient [GiniD]), and PBS. Soil properties—including physical (soil moisture), chemical (soil organic carbon [SOC], total nitrogen [TN], total phosphorus [TP], ammonium nitrogen [NH₄⁺], nitrate nitrogen [NO₃⁻], available phosphorus [AP]), microbial (microbial biomass carbon [MBC], nitrogen [MBN], and phosphorus [MBP]), and enzymatic (e.g., peroxidase [POD], alkaline phosphatase [ALP], urease [URE])—also varied significantly. SQI values ranged from 0.42 to 0.65, with patterns Fir–Mytilaria laosensis mixed (ML), Fir–Castanopsis hystrix mixed (CH), Fir–Michelia chapensis mixed (MC), and Fir–Schima superba mixed (SS) associated with both high SQI and greater biomass. Sensitivity analysis identified Fir–Cinnamomum porrectum mixed (CP), ML, and SS as particularly responsive to hybridization. Among soil factors, URE, AP, and MBC were key drivers of productivity, while URE, AP, MBC, and POD significantly predicted the proportion of broadleaf trees. Enhanced soil quality was positively associated with increases in DBH, TH, and PBS, accelerating the successional transition from fir-dominated to broadleaf-dominated stands. However, SQI was not significantly correlated with structural diversity metrics.

Discussion: These results underscore the importance of rational species selection in restoring degraded fir plantations and demonstrate that improving soil quality is a critical mechanism promoting near-natural forest succession.

1 Introduction

Soil is fundamental for plant growth and development (Li et al., 2020). Soil quality directly influences the biological productivity and functioning of forests (Liang et al., 2022). Furthermore, soil quality has been defined as "the ability of soil to sustain plant and animal productivity, enhance water and air quality, and support human health and habitation within natural ecosystems" (Li et al., 2019a). Consequently, evaluating soil quality is crucial for the advancement of sustainable forestry practices (Yu et al., 2023).

Since the release of the United States Department of Agriculture (USDA) classification system for land potential in 1961, numerous soil quality assessment methods have been developed (Atalay, 2016). The Soil Quality Index (SQI) is a widely utilized method for assessing soil function (Li et al., 2019b; Qian et al., 2023). Soil function generally can be evaluated through indicators of soil physicochemical and biochemical properties (Bünemann et al., 2018). However, the extensive range of physical, chemical, and biological properties complicates the measurement process, making it impractical to consider all these properties (Li et al., 2019b). Currently, the minimum data set (MDS), developed by Larson and Pierce, is arguably the most commonly employed method for evaluating SQI (Li et al., 2019a; Raiesi and Beheshti, 2022; Qian et al., 2023, Qian et al., 2023). This approach comprehensively accounts for the combined influence of measurements, weights, and interactions between indicators on the assessment results (Rezaei et al., 2006). Mixed forests markedly enhance soil quality by improving soil structure, augmenting nutrient availability, and increasing fungal diversity along with the stability of microbial networks (Ge et al., 2025). This advantage over pure coniferous forests is consistently observed across various ecosystems, with coniferous-broadleaf mixtures generally proving superior in maintaining and enhancing soil quality (Shao et al., 2020). This trend is also evident in karst-degraded regions, where mixed forests significantly increase soil organic carbon (SOC) and total nitrogen (TN) levels, thereby improving overall soil fertility and quality (Guan and Fan, 2020). Collectively, the evidence indicates that mixed forests achieve significantly higher Soil Quality Index (SQI) values compared to coniferous forests, highlighting the importance of mixed afforestation in sustaining soil quality.

Planted forests are recognized as an effective strategy for improving the quality of degraded soils and mitigating soil erosion throughout the world (Guo et al., 2021). Chinese Fir, as a major plantation species in southern China, is characterized by its straight stem shape, fast-growing and productive nature, strong adaptability, and ease of management and harvesting (Guo et al., 2022). Chinese Fir is extensively cultivated across 16 provinces in China (Li et al., 2024), with fir plantation forests comprising 24% of the total area of plantation forests in the country and 6% of plantation forests worldwide (Hemati et al., 2020). However, studies have reported that continuous cultivation and short-rotation harvesting associated with fir plantations are diminishing soil quality (Farooq, 2019). To date, extensive literature indicates that monoculture plantation forests encounter issues such as low biodiversity, as well as declining soil fertility and water-holding capacity (Chen et al., 2021; Li et al., 2021), which severely impede the advancement of modern sustainable forestry (Guo et al., 2023). Consequently, the establishment of mixed forests has become a central focus of various forest management strategies (Bauhus et al., 2017). The introduction of broadleaf trees into coniferous forests represents a crucial approach to enhancing carbon stocks within forest stands (Feng et al., 2022).

Incorporating broadleaf species into Fir plantations can improve soil physicochemical properties, accelerate soil nutrient turnover, and enhance stand productivity (Hou et al., 2021; Gao et al., 2022; Guo et al., 2023). Moreover, while the benefits of species mixing in enhancing forest structural stability and productivity are well documented (Hua et al., 2022), the impact on soil quality and fertility (as measured by the soil quality index) resulting from the conversion of cedar monoculture plantations to mixed forests remains largely unexplored in southern China (Guo et al., 2023). Therefore, there is an urgent need to evaluate soil quality, as maintaining or improving soil fertility is essential for sustainable forestry (Ye et al., 2022).

Forest stand modification through the introduction of different tree species is a common strategy in forestry production and management (Liu et al., 2016). Such modifications may lead to changes in vegetation diversity and shifts in ecological niche utilization, which can, in turn, affect forest structure and ecosystem stability (Guo et al., 2023). Stand structure plays a crucial role in predicting forest growth and productivity (Wang et al., 2021). The diversity of stand structure, characterized by variability in individual tree size, such as the Gini coefficient and the coefficient of variation in tree diameter (Ali, 2019) is increasingly recognized as a significant driver of ecosystem functioning (LaRue et al., 2019; Ouyang et al., 2023). Tree mortality is a fundamental process in forest succession (Zhang et al., 2020). Stand structure influences variations in competitive intensity and resource utilization, thereby affecting both tree mortality (Kweon and Comeau, 2019) and forest succession (Zhang et al., 2020). Current studies have raised concerns regarding the relationship between stand structure and productivity in mixed fir stands (Wang et al., 2021), the impacts of changes in stand density on stand structure (Chen et al., 2021; Liu et al., 2022; Sun et al., 2023), and the effects of stand structure and climatic factors on tree mortality (Zhang et al., 2020). However, the influence of changes in soil quality on stand structure and forest succession has been infrequently investigated. Exploring the relationship between soil quality, stand structure, and forest succession will not only provide a foundation for the scientific management of Fir plantation forests and the adjustment of local forestry policies, but will also serve as a reference for the sustainable development of the plantation forest industry globally.

Different tree species configurations were selected for stand rejuvenation at a fourth-generation Fir forest logging site in the southern subtropics. This paper employs the model for soil quality evaluation (MDS) to assess the differences in soil quality among various stand types. Additionally, it explores the relationships between soil quality and stand growth, biomass, structural diversity, and community succession to broadleaf species, expressed as a percentage of broadleaf species. The study hypothesizes that improved soil quality is positively associated with stand growth, biomass, and structural diversity, and soil quality improvement mediates tree size dynamics (DBH, TH) to accelerate competitive exclusion of fir by broadleaf species in mixed stands.

2 Materials and methods

2.1 Study area descriptions

The study was conducted in Lechang Forest Farm, located in Guangdong Province, China (25°06'–25°19'N, 113°11'–113°23'E; Figure 1). The site spans 5,398.3 hectares and experiences a meso-subtropical monsoon climate, with a mean annual temperature of 19.6 °C, average annual precipitation of 1,522.3 mm, and mean annual relative humidity of 76%.

Figure 1
Map showing Lechang City in China, with detailed terrain features, highlighting sampling points, major railways, national highways, expressways, and main rivers. Insets display the city's location within Guangdong province and China.

Figure 1. Map of the study area.

The Fir plantation forest in the study area has been cultivated for four generations, resulting in significant soil fertility decline and general forest degradation. The area was last logged in 2003, after which sprouting vegetation was retained. Following the logging event, the affected areas were replanted with a variety of broadleaf species, including Mytilaria laosensis, Michelia macclurei, Michelia maudiae, Cinnamomum porrectum, Castanopsis hystrix, Cinnamomum burmanni, Liquidambar formosana, Michelia chapensis, Schima superba, and Michelia odora, among others. In total, ten broadleaf and coniferous species were introduced to establish mixed stands, with broadleaf species deliberately interplanted among the remaining firs.

Land preparation involved band reclamation and the construction of weed-branch berms, while weeding and maintenance were conducted during the first three years after planting. This process led to the development of a mixed forest stand composed of both modified tree species and naturally regenerating fir trees, with an initial planting density of 3,975 trees per hectare.

In July 2023, representative stands of each of the ten tree species were selected, and three permanent plots (20 m × 30 m) were established for each species, with a minimum spacing of 20 m between plots. Within each plot, all trees were surveyed individually. Recorded data included species identity, diameter at breast height (DBH), tree height, height to first branch, spatial coordinates, and general growth condition. Further plot details are provided in Table 1.

Table 1
www.frontiersin.org

Table 1. Basic characteristics of different tree species plantations for sample plots.

2.2 Soil sampling and soil properties analysis

Each 20 × 30 m sample plot was subdivided into six 10 × 10 m subplots. After removing surface litter, soil was sampled. Using a 5 cm diameter auger, nine surface soil cores (0–10 cm) were collected in an S-shaped pattern within each subplot. This sampling depth was selected as it represents the primary zone influenced by litter decomposition, root activity, and microbial processes, and it is the standard for soil quality assessments in subtropical plantations, where fine roots and nutrient cycling are predominantly concentrated. These cores were combined to form a composite sample. Additionally, three undisturbed soil samples were collected per subplot using a ring knife to determine bulk density and gravimetric moisture content. All soil samples were passed through a 2 mm sieve after removing visible roots and gravel. The sieved samples were stored at low temperature in a mobile refrigerator and transported to the laboratory, where each sample was divided into three portions: (1) air-dried for physicochemical analysis and amino sugar quantification; (2) refrigerated at –4 °C for biological analysis; and (3) retained as backup. Detailed analytical procedures for soil physicochemical and biological measurements are provided in Appendix 1.

2.3 Structural diversity and biomass of different mixed planting patterns

Structural diversity was quantified using variation in tree diameter at breast height (DBH), a standard metric for stand structure (Brejda et al., 2000; Schnabel et al., 2019). DBH variation was assessed by calculating the coefficient of variation (CV), defined as the standard deviation of DBH divided by the mean DBH (Wan et al., 2024). The Gini coefficient was also calculated to evaluate inequality in DBH distribution (Cordonnier and Kunstler, 2015).

GiniD=100%1n(DBHiDBH)2DBH

where GiniD​ is the Gini coefficient, n is the number of trees in the plot, DBHi is the DBH of the ith tree, and DBH is the mean DBH.

Due to large variation in tree size, species composition was quantified using biomass proportion rather than tree counts per hectare. Aboveground biomass was estimated using species-specific allometric equations incorporating DBH and tree height (Návar, 2009):

W=a(D2H)b

Where W is biomass (kg), D is DBH (cm), H is height (m), and a and b are species-specific coefficients (see Appendix II). Total plot biomass was calculated by summing individual tree biomass and scaled to per-hectare values. Net primary productivity was estimated as the average annual increment of aboveground biomass, defined as total biomass divided by stand age (Luo and Chen, 2011).

2.4 Soil quality assessment

Soil quality was assessed using the Soil Management Assessment Framework (SMAF), with the Soil Quality Index (SQI) calculated based on a Minimum Data Set (MDS) derived from a Total Data Set (TDS) of 20 soil indicators (Yu et al., 2023). Indicators were chosen to represent physical, chemical, and biological soil functions relevant to nutrient cycling, carbon dynamics, and structural stability (Nosrati, 2013; Guo et al., 2017; He et al., 2022). The TDS included soil moisture (SM), bulk density (BD), soil organic carbon (SOC), total nitrogen (TN), total and available phosphorus (TP, AP), ammonium (NH4+) and nitrate (NO3-) nitrogen, microbial biomass carbon (MBC), nitrogen (MBN), and phosphorus (MBP), and nine enzymes involved in C, N, and P cycling: peroxidase (POD), β-xylosidase (BX), cellobiohydrolase (CBH), β-glucosidase (BG), urease (URE), N-acetyl-β-D-glucosaminidase (NAG), leucine aminopeptidase (LAP), acid phosphatase (ACP), and alkaline phosphatase (ALP).

Kaiser–Meyer–Olkin (KMO = 0.704) and Bartlett’s test (p < 0.01) confirmed the suitability of the dataset for Principal Component Analysis (PCA). Indicators with loadings ≥ 0.5 on principal components (PCs) with eigenvalues ≥ 1 were grouped by PC and ranked using a normalized contribution metric:

Nik=i=1k(Uik2·ek)

Where Nik is the normalized loading of indicator I on PC k, Uik is the indicator loading, and ek​ is the PC eigenvalue. Indicators within 10% of the highest Nik were identified as highly weighted (Yu et al., 2023). If multiple high-weight indicators were not significantly correlated (Pearson correlation), all were retained in the MDS; if correlated, only the indicator with the highest norm was retained (Li et al., 2013; Jin et al., 2021). MDS indicators were scored using a nonlinear scoring function (SNL), preferred for its ability to reflect soil functional responses (Raiesi and Beheshti, 2022):

SNL=a1+(xxm)b

Where a = 1 is the maximum score, x is the observed value, xm is the mean value, and b = −2.5 or 2.5 for "more is better" or "less is better" indicators, respectively (Zhang et al., 2011; Yu et al., 2018). Most indicators were scored as "more is better," except for bulk density (BD), which was scored as "less is better."

The overall Soil Quality Index (SQI) was calculated using a weighted additive model:

SQI=i=1nWi·Si

Where Wi is the weight of indicator i based on its contribution to PC variance, and Si is its SNL score. Indicator weights were calculated as:

Wi=Ci1nCi

With Ci representing the common factor variance of indicator i (Rahmanipour et al., 2014). The SQI was normalized to a 0–1 scale and categorized into five classes: very high (>0.55), high (0.46–0.55), medium (0.37–0.46), low (0.28–0.37), and very low (<0.28), based on national standards for soil fertility classification (ISMAPRC, 1996).

To evaluate the sensitivity of SQI to different hybridization models, a sensitivity index (SI) was calculated:

SI=SQImaxSQImin

Higher SI values indicate greater responsiveness of soil quality to hybridization pattern (Mamehpour et al., 2021).

2.5 Statistical analysis

Data were processed and analyzed using Microsoft Excel 2010 and SPSS 25.0 (Statistical Graphics Corp., Princeton, USA). The Shapiro–Wilk test was used to assess data normality, and Levene’s test was applied to evaluate the homogeneity of variances. When assumptions of normality or homoscedasticity were violated, data were log-transformed prior to analysis.

One-way analysis of variance (ANOVA) was conducted to examine differences in diameter at breast height (DBH), tree height (TH), percentage of broadleaf species (PBS), coefficient of variation of DBH (CVD), Gini coefficient of DBH (GiniD), forest biomass (FB), net primary productivity (NPP), and soil properties across different mixed planting patterns. Duncan’s multiple range test was used for post hoc comparisons. All results are presented as mean ± standard error (SE), based on three replicates (n = 3). Pearson correlation and Mantel tests were employed to evaluate the relationships between soil properties, soil quality, stand structure, and productivity. Partial least squares regression (PLSR) was used to assess the relationship between the soil quality index (SQI), stand structural attributes, and the proportion of broadleaf species. Random forest analysis was conducted to determine the relative importance of soil properties in explaining variation in stand productivity and the proportion of broadleaf species. A partial least squares path modeling (PLS-PM) analysis was performed using the "plspm" package in R to investigate the key drivers influencing net primary productivity.

3 Results

3.1 Structural diversity and forest productivity

The values of DBH (6.52–21.34 cm), TH (7.96–22.55 m), GiniD (0.83–4.54), NPP (1.62–14.38 t·hm-²·a-¹), PBS (27.13%–100%), FB (82.78–274.18 t·hm-²), and CVD (21.22%–46.17%) varied significantly across different mixed planting patterns (p< 0.05; Figures 2a, b). Notably, the highest values of FB, DBH, TH, and PBS were all observed in the ML planting pattern, while the highest CVD and GiniD values were recorded in CH. The lowest values of DBH and PBS were found in MME, whereas the lowest FB, TH, CVD, and GiniD values were observed in CB, MO, ML, and MM, respectively. The NPP of the ML planting pattern was significantly higher than that of MO (Figures 2a, b).

Figure 2
Circular bar charts and a correlation matrix visualize data on forest metrics, including diameter and height. Charts (a) and (b) use different colors to represent diameter, height, productivity, species percentage, biomass, and variation coefficients. Correlation matrix (c) highlights statistical significance with asterisks.

Figure 2. Effects of ten mixed planting patterns (ML, Fir–Mytilaria laosensis mixed; MM, Fir–Michelia macclurei mixed; MME, Fir–Michelia maudiae mixed; CP, Fir–Cinnamomum porrectum mixed; CH, Fir–Castanopsis hystrix mixed; CB, Fir–Cinnamomum burmanni mixed; LF, Fir–Liquidambar formosana mixed; MC, Fir–Michelia chapensis mixed; SS, Fir–Schima superba mixed; MO, Fir–Michelia odora mixed.) on forest structural diversity and productivity. (a) Net primary productivity (NPP), diameter at breast height (DBH), tree height (TH), and the Gini coefficient of DBH (GiniD); (b) percentage of broadleaf species (PBS), coefficient variation of DBH (CVD), and forest biomass (FB); (c) Pearson correlation matrix depicting the relationships among all stand structural and productivity-related variables. * indicates that a correlation is significant (P < 0.05); ** indicates that a correlation is extremely significant (P < 0.01).

Correlation analysis revealed that FB was significantly positively correlated with DBH, TH, and PBS. In addition, DBH showed significant positive correlations with both TH and PBS, and PBS was also positively correlated with TH. A strong positive correlation was found between GiniD and CVD. In contrast, CVD was significantly negatively correlated with DBH, TH, and FB, and GiniD showed a significant negative correlation with FB. NPP was significantly correlated with all variables except GiniD. Specifically, NPP showed significant positive correlations with DBH, TH, PBS, and FB, and a significant negative correlation with CVD (Figure 2c).

3.2 Soil properties of different mixed planting patterns

The results of the one-way ANOVA indicated that most physical (SM), chemical (TN, TP, SOC, NH4+, NO3-, AP), microbial (MBC, MBN, MBP), and enzyme activity (POD, LAP, BX, ALP, URE, CBH, NAG, ACP, and BG) characteristics, with the exception of BD, exhibited significant differences among the 10 mixed planting patterns (Table 2).

Table 2
www.frontiersin.org

Table 2. The soil properties of ten 10 mixed planting patterns (mean ± stand and deviation, n=3).

3.3 Soil quality assessment

Principal component analysis (PCA) identified six principal components with eigenvalues greater than 1, which were selected for constructing the Minimum Data Set (MDS). The eigenvalues of the components were as follows: PC1 = 6.366, PC2 = 2.196, PC3 = 1.948, PC4 = 1.759, PC5 = 1.368, and PC6 = 1.101. Collectively, these six components explained 73.68% of the total variance in the dataset (Table 3).

Table 3
www.frontiersin.org

Table 3. Load matrix communality value and norm values of assessment indicators.

Based on PCA results, all 20 soil indicators were grouped by principal component as follows:

● Group 1: SM, BD, TN, SOC, MBC, MBN, MBP, AP, LAP, ALP, URE

● Group 2: NO3-, CBH, ACP, BG

● Group 3: NH4+

● Group 4: POD, BX

● Group 5: NAG

● Group 6: TP

From each group, indicators with high loadings (within 10% of the group’s highest norm value) were initially retained. These included:

● Group 1: ALP, SOC, TN, MBN, MBC

● Group 2: NO3-, CBH, ACP, BG

● Groups 3–6: NH4+, POD, BX, NAG, TP (each the sole indicator in its group)

Correlation analysis (Figure 3a) was then used to further refine indicator selection. In Group 1, ALP showed significant correlations with SOC (r = 0.40*), TN (r = 0.37*), MBC (r = 0.59**), and MBN (r = 0.62**), and was therefore selected as the representative indicator due to its highest norm value. In Group 2, NO3- and CBH were significantly negatively correlated (r = –0.38*), so NO3-, ACP, and BG were retained. In Group 4, POD and BX were not significantly correlated and were both included. In Groups 3, 5, and 6, NH4+, NAG, and TP were each retained by default as the only high-loading indicators in their respective groups.

Figure 3
A complex diagram with four panels illustrating soil analysis data. (a) A correlation matrix with colored squares indicating Pearson's r values for various soil indicators, with a legend showing significance levels. (b) A radar chart displaying weight assignment of soil indicators in MDS, with values for each indicator. (c) A radar chart depicting the soil sensitivity index with specific measurements labeled. (d) A radial plot of soil quality index, showing different categories and statistical significance with labeled sections. The image uses color coding and lines to convey relationships and data significance.

Figure 3. The Pearson correlation between soil indicators and the Mantel test concerning the soil quality index, tree growth metrics (DBH, TH, CVD, and GiniD), and net primary productivity (NPP) in relation to soil properties (a), weight assignment soil indicators in MDS (b), soil sensitivity index (c) and Soil quality index (d). The figure (d) value reported as mean ± standard deviation (n = 3). Significant differences at p < 0.05 are indicated by different lowercase letters. Classification of soil quality, “very high” (SQI>0.55), “high” (0.46-0.55), “medium” (0.37-0.46), “low” (0.28-0.37) and “very low” (SQI<0.28) (Guo et al., 2017).

Ultimately, the final MDS consisted of nine indicators: ALP, NO3-, ACP, BG, POD, BX, NH4+, NAG, and TP. The assigned weights for these indicators were 0.120, 0.106, 0.086, 0.114, 0.108, 0.110, 0.123, 0.121, and 0.112, respectively (Figure 3b).

The soil sensitivity indices (SI) for the ten mixed planting patterns ranged from 1.002 to 1.496, with CP exhibiting the highest sensitivity and CH the lowest. Other planting patterns with relatively high SI values included ML (1.340), CB (1.327), and SS (1.217) (Figure 3c). Significant differences in Soil Quality Index (SQI) were observed among the ten planting patterns (p< 0.05), with values ranging from 0.42 to 0.65. CH had a significantly higher SQI than both MME and MO (Figure 3d). According to the SQI classification (Guo et al., 2017), all ten mixed planting patterns were rated at or above the medium level: ML, MM, CH, MC, and SS were classified as very high; CP, CB, and LF as high; and MME and MO as medium.

Among the soil properties, urease activity (URE), available phosphorus (AP), and microbial biomass carbon (MBC) had the strongest influence on stand productivity. Additionally, URE, AP, MBC, and Peroxidase (POD) were significant predictors of the proportion of broadleaf species in mixed stands (Figures 4a, b).

Figure 4
Two bar charts and six scatter plots illustrate ecological data relationships. Charts (a) and (b) show factors affecting net primary productivity and broad-leaved species percentage, respectively, with significant values marked. Scatter plots (c) to (h) show relationships between SQI and various ecological parameters, including DBH, TH, CV, Girin, NPP, and PBS, each fitted with a regression line, displaying equations, R-squared, and P-values. Each plot is color-coded for clarity.

Figure 4. Relative importance of soil variables in predicting variations in (a) net primary productivity and (b) the percentage of broad-leafed species in mixed plantations. The least-squares regression analysis includes SQI and the following variables: (c) DBH, (d) TH, (e) CVD; (f) GiniD, (g) NPP, (h) PBS. All data were logarithmically transformed prior to the fitting analysis. Statistically significant relationships (P < 0.05) are indicated, with the 95% confidence intervals represented by the shaded areas. Asterisks denote significant correlations at the * P < 0.05, ** P < 0.01, and *** P < 0.001 levels, respectively.

Stand structural diversity, as measured by the coefficient of variation in DBH (CVD) and the Gini coefficient (GiniD), showed no significant correlation with SQI (Figures 4e, f). In contrast, tree height (TH), DBH, net primary productivity (NPP), and the proportion of broadleaf species increased significantly with higher SQI values (p< 0.01; Figures 4c, d, g, h).

The PLS-PM model confirms that mixed planting patterns do not directly enhance net primary productivity to a great extent but rather function by initiating a key ecological succession process. Specifically, these patterns improve the soil quality and promote tree growth, which in turn drives the stand composition toward a higher percentage of broadleaf species. Ultimately, this increased percentage of broadleaf species acts as the most powerful direct driver, significantly boosting the forest's net primary productivity. This finding underscores the pivotal role of species composition, particularly the proportion of broadleaf species, in shaping the productivity of mixed forests (Figure 5).

Figure 5
A structural equation model diagram on the left illustrates relationships between mixed planting patterns, soil quality index, tree growth, percentage of broadleaf species, and net primary productivity, with arrows showing direct paths and standardized coefficients. The right side is a bar chart comparing direct, indirect, and total effects for each variable. The goodness-of-fit is 0.74, with significant paths marked by asterisks indicating different levels of statistical significance.

Figure 5. Partial least squares path modeling (PLS-PM) describing the direct and indirect effects of mixed planting patterns on net primary productivity. The model demonstrates a good goodness of fit (GOF = 0.74). Path coefficients are displayed alongside the arrows, with significance levels indicated as ** (P < 0.01) and *** (P < 0.001). The values represent the proportion of variance explained for each endogenous variable. Gray single-headed arrows denote direct effects on net primary productivity, while blue arrows represent the direct effects of mixed planting patterns on mediating variables (soil quality index, tree growth [DBH, TH, CVD, GiniD], percentage of broadleaf species), which subsequently exert indirect effects on net primary productivity.

4 Discussion

4.1 Growth and soil characteristics of forest stands

Previous studies have shown that mixed-species plantations can enhance ecosystem productivity and improve ecosystem functions and services relative to the original single-species plantation (Guo et al., 2023). Near-naturalization of coniferous forests—through structural transformation of homogeneous monocultures into mixed stands—has been proposed as a strategy to achieve this (Hou et al., 2021; Guo et al., 2022). In this study, different combinations of fir and broadleaf species resulted in substantial differences in stand growth (DBH and TH), structural diversity (CVD and GiniD), and total biomass (Figure 2). These differences likely stem from interspecific variation in shade tolerance and growth rates, which leads to vertical stratification of the forest canopy (Thurm and Pretzsch, 2016).

Interestingly, stand structural diversity in this study was negatively correlated with both stand growth and biomass (Figure 2). This contrasts with findings from other studies reporting positive relationships between structural diversity and aboveground biomass, including in temperate mixed forests of northeastern China (Yuan et al., 2018), subtropical secondary forests (Ali et al., 2016), and large-scale tropical forests (Ali et al., 2019). For example, Zhang and Chen (2015) reported that tree diameter diversity was positively correlated with aboveground biomass, and Dănescu et al. (2016) highlighted the direct contribution of structural diversity to productivity. The discrepancy between our results and previous findings may be attributed to differences in community composition, species-specific traits, and environmental conditions (Tetemke et al., 2021; Ullah et al., 2021). It has also been proposed that negative relationships between structural diversity and productivity can emerge under competitive exclusion and selection effects, particularly in less-disturbed or high-productivity environments (Paquette and Messier, 2011).

Most soil physical, chemical, microbial, and enzymatic indicators—excluding bulk density (BD)—exhibited significant variation among the mixed planting types (Table 2), but no single forest type consistently exhibited the highest soil values. This suggests that the effects of species mixing on soil properties are not uniform. One potential explanation is that some of the broadleaf evergreens introduced have relatively low litter inputs, influencing soil structure and function (Zhou et al., 2020). Past studies have similarly shown that shifts in soil physicochemical characteristics depend heavily on the vegetation types used to restore degraded lands (Hou et al., 2019), and species-specific effects in mixed stands can further mediate these changes (Stewart et al., 2011). Overall, this underscores the importance of rational species selection in reforestation strategies to enhance soil improvement in degraded fir plantations.

4.2 Soil quality index and forest succession

Monitoring and evaluating soil quality is a critical component of sustainable forest management (Huang et al., 2018). The Soil Quality Index (SQI), developed from a multivariate indicator set, provides a quantitative framework to evaluate the impact of management practices on soil health (Chen et al., 2021; He et al., 2022). In this study, an MDS of nine indicators—including three physicochemical properties (TP, NO3-, NH4+) and six enzyme activity metrics (ALP, ACP, BG, POD, BX, NAG)—was selected to simplify and strengthen the SQI assessment. The inclusion of enzyme indicators involved in C, N, and P cycling reflects their value in capturing biologically mediated changes in soil function. These findings refine and support previous calls to include enzyme activities in soil quality assessments for cedar forests.

The observed SQI values (0.42–0.65) were lower than those previously reported in other fir-broadleaf plantations (0.66–0.85; Huang et al., 2018), likely due to differences in climate, stand age, or species composition (Shao et al., 2020; Cao et al., 2023). Among the mixed planting patterns, CH, ML, MM, MC, and SS exhibited significantly higher SQI values (Figure 3d). Notably, the ML, MC, and SS patterns also had high stand biomass (Figure 2a), suggesting that these combinations not only enhance soil quality but also promote productivity. Although CH had lower biomass, it featured the highest SQI and exhibited large DBH and tree height values (Figures 2b, c), likely due to low stand density and dominance by larger individuals (Table 1). This suggests that CH may be a suitable planting model for cultivating large-diameter timber in degraded fir plantations. In contrast, the MM pattern exhibited higher soil quality but relatively low DBH and tree height (Figure 2), making it less suitable for timber production despite good soil performance. These differences highlight the importance of selecting mixed-species strategies tailored to specific management goals—whether soil restoration or timber yield.

Forest structural attributes, including variation in tree size (DBH and height), are known to influence species diversity, regeneration dynamics, and forest function (Clark, 2010; Yuan et al., 2018; Hui et al., 2019; Põldveer et al., 2021). In this study, the proportion of broadleaf species differed significantly across mixed stands (Figure 2d), likely due to species-specific requirements for light and water. High structural diversity, as seen in CH and MC, can enhance ecological niche complementarity through improved light utilization. Conversely, lower structural complexity, such as in ML, may reduce resource-use efficiency (Ali et al., 2016).

The transition from fir- to broadleaf-dominant stands was closely tied to both tree size and soil quality. Improved soil quality significantly promoted DBH and height (Figure 4), which in turn increased the percentage of broadleaf trees in mixed stands. These findings are consistent with the idea that nutrient-rich soils facilitate growth and canopy development (Poorter et al., 2017), leading to changes in vertical structure and light interception (Onoda et al., 2014).

Forest structure also mediates competition and regeneration. As tree size variation increases, dominant individuals face reduced competition while suppressed individuals experience intensified resource limitation (Ali and Yan, 2017; Zhang et al., 2017). Size-asymmetric competition—particularly for light—is a major mechanism driving tree mortality and skewed size distributions (Bourdier et al., 2016; Soares et al., 2016). In our study, although the initial planting densities were uniform, stand density diverged considerably over time (Table 1), highlighting the dynamic nature of post-establishment competition. Fir, being shade-intolerant, initially grows faster than broadleaf species, but is eventually outcompeted under closed-canopy conditions as broadleaf trees capture more light and water due to their larger canopies (Zhang et al., 2020).

As upper-canopy trees escape light limitation, suppressed fir individuals face increasingly severe shading. This can lead to carbon starvation, where photosynthetic gains are insufficient to cover respiratory costs (Luo and Chen, 2011; Zhang et al., 2020). Larger Gini coefficients reflect increased size inequality and intensified asymmetric competition, which can accelerate mortality of smaller individuals. This is supported by our observation that ML, which had larger average DBH and tree height, also exhibited lower Gini coefficients—suggesting a more uniform, less competitive size structure.

One limitation of this study is the absence of baseline soil data from the plantation establishment period (circa 2003). This precludes us from definitively ruling out the potential influence of initial soil heterogeneity on the observed community trajectories. Moreover, the strong correlations between current soil properties (e.g., URE, AP, MBC) and stand attributes (DBH, TH, PBS) likely reflect bidirectional plant-soil feedbacks rather than a unidirectional soil effect. For instance, introducing broadleaf species can initiate positive feedback cycles by enhancing nutrient availability and microbial activity via litter and root exudates, which in turn promotes their growth and dominance (Zhou et al., 2020; Guo et al., 2023). Nevertheless, significant differences in SQI and stand structure emerged after two decades from a common starting point—a uniformly degraded fir plantation. Coupled with the identification of key, dynamic soil drivers (e.g., URE, MBC) that are responsive to vegetation change, this provides compelling evidence that soil quality enhancement is a critical mechanism driving successional transition. Thus, our findings reveal a strong, functionally important plant-soil feedback loop that governs succession in these mixed plantations, rather than a simple unidirectional pathway.

5 Conclusions

This study evaluated soil properties, stand growth, structural diversity, biomass, and the proportion of broadleaf species across various fir–broadleaf mixing patterns. The results demonstrated that differences in soil properties were primarily influenced by species composition. Among the planting patterns, ML, CH, MC, and SS significantly improved both soil quality and stand biomass. Enhanced soil quality played a key role in promoting tree growth, thereby accelerating the near-natural successional transition from fir-dominated stands to broadleaf forests. For the conversion of degraded Chinese fir plantations, the species MC, ML, and CH are proposed as suitable candidates to enhance ecosystem productivity and soil quality. These findings underscore the importance of rational species selection in rehabilitating degraded cedar plantations and highlight soil quality as a critical driver of successional dynamics in mixed stands. This work offers new insights into the management and restoration of degraded fir plantations and presents a foundation for the development of ecologically informed hybridization strategies.

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

ZF: Investigation, Visualization, Writing – original draft, Writing – review & editing. YX: Investigation, Visualization, Writing – original draft, Writing – review & editing. SL: Conceptualization, Writing – review & editing. HX: Software, Writing – review & editing. YW: Investigation, Writing – review & editing. HZ: Investigation, Validation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was jointly supported by the National Key Research and Development Program of China (2022YFF1303003-06) and the Guangdong Forestry Sci-Tech Innovation Platform Project (2022CG646 & 2024CG259).

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.

Generative AI statement

The author(s) declare that no Generative AI was 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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1705626/full#supplementary-material

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

Ali, A. (2019). Forest stand structure and functioning: Current knowledge and future challenges. Ecol. Indic. 98, 665–677. doi: 10.1016/j.ecolind.2018.11.017

Crossref Full Text | Google Scholar

Ali, A., Lin, S.-L., He, J.-K., Kong, F.-M., Yu, J.-H., and Jiang, H.-S. (2019). Climate and soils determine aboveground biomass indirectly via species diversity and stand structural complexity in tropical forests. For. Ecol. Manage. 432, 823–831. doi: 10.1016/j.foreco.2018.10.024

Crossref Full Text | Google Scholar

Ali, A. and Yan, E.-R. (2017). The forest strata-dependent relationship between biodiversity and aboveground biomass within a subtropical forest. For. Ecol. Manage. 401, 125–134. doi: 10.1016/j.foreco.2017.06.056

Crossref Full Text | Google Scholar

Ali, A., Yan, E.-R., Chen, H. Y. H., Chang, S. X., Zhao, Y.-T., Yang, X.-D., et al. (2016). Stand structural diversity rather than species diversity enhancesaboveground carbon storage in secondary subtropical forests in Eastern China. Biogeosciences 13, 4627–4635. doi: 10.5194/bg-13-4627-2016

Crossref Full Text | Google Scholar

Atalay, I. (2016). A new approach to the land capability classification: case study of Turkey. Proc. Environ. Sci. 32, 264–274. doi: 10.1016/j.proenv.2016.03.031

Crossref Full Text | Google Scholar

Bauhus, J., Forrester, D. I., and Pretzsch, H. (2017). “Mixed-species forests: the development of a forest management paradigm,” in Mixed-species forests. Eds. Pretzsch, H., Forrester, D. I., and Bauhus, J. (Springer Berlin Heidelberg), 1–25, Berlin, Heidelberg). doi: 10.1007/978-3-662-54553-9_1

Crossref Full Text | Google Scholar

Bourdier, T., Cordonnier, T., Kunstler, G., Piedallu, C., Lagarrigues, G., and Courbaud, B. (2016). Tree size inequality reduces forest productivity: an analysis combining inventory data for ten european species and a light competition model. PloS One 11, e0151852. doi: 10.1371/journal.pone.0151852

PubMed Abstract | Crossref Full Text | Google Scholar

Brejda, J. J., Moorman, T. B., Karlen, D. L., and Dao, T. H. (2000). Identification of regional soil quality factors and indicators I. Central and southern high plains. Soil Sci. Soc Am. J. 64, 2115–2124. doi: 10.2136/sssaj2000.6462115x

Crossref Full Text | Google Scholar

Bünemann, E. K., Bongiorno, G., Bai, Z., Creamer, R. E., De Deyn, G., De Goede, R., et al. (2018). Soil quality – A critical review. Soil Biol. Biochem. 120, 105–125. doi: 10.1016/j.soilbio.2018.01.030

Crossref Full Text | Google Scholar

Cao, X., Zhang, Z., Yan, W., Mo, Y., Wu, S., Wang, M., et al. (2023). Assessment of soil quality in the transformation from pure chinese fir plantation to mixed broad-leaved and cunninghamia lanceolata plantation in subtropical China. Forests 14, 1867. doi: 10.3390/f14091867

Crossref Full Text | Google Scholar

Chen, Q., Liu, Z., Zhou, J., Xu, X., and Zhu, Y. (2021). Long-term straw mulching with nitrogen fertilization increases nutrient and microbial determinants of soil quality in a maize–wheat rotation on China’s Loess Plateau. Sci. Total Environ. 775, 145930. doi: 10.1016/j.scitotenv.2021.145930

Crossref Full Text | Google Scholar

Clark, J. S. (2010). Individuals and the variation needed for high species diversity in forest trees. Science 327, 1129–1132. doi: 10.1126/science.1183506

PubMed Abstract | Crossref Full Text | Google Scholar

Cordonnier, T. and Kunstler, G. (2015). The Gini index brings asymmetric competition to light. Perspect. Plant Ecol. Evol. Syst. 17, 107–115. doi: 10.1016/j.ppees.2015.01.001

Crossref Full Text | Google Scholar

Dănescu, A., Albrecht, A. T., and Bauhus, J. (2016). Structural diversity promotes productivity of mixed, uneven-aged forests in southwestern Germany. Oecologia 182, 319–333. doi: 10.1007/s00442-016-3623-4

PubMed Abstract | Crossref Full Text | Google Scholar

Farooq, T. H. (2019). Chinese fir (cunninghamia lanceolata) a green gold of China with continues decline in its productivity over the successive rotations: a review. Appl. Ecol. Environ. Res. 17, 11055–11067. doi: 10.15666/aeer/1705_1105511067

Crossref Full Text | Google Scholar

Feng, Y., Schmid, B., Loreau, M., Forrester, D. I., Fei, S., Zhu, J., et al. (2022). Multispecies forest plantations outyield monocultures across a broad range of conditions. Science 376, 865–868. doi: 10.1126/science.abm6363

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, G., Huang, X., Xu, H., Wang, Y., Shen, W., Zhang, W., et al. (2022). Conversion of pure Chinese fir plantation to multi-layered mixed plantation enhances the soil aggregate stability by regulating microbial communities in subtropical China. For. Ecosyst. 9, 100078. doi: 10.1016/j.fecs.2022.100078

Crossref Full Text | Google Scholar

Ge, Z., Zhang, X., Liu, C., Li, M., Wang, R., Zhang, Y., et al. (2025). Microbial determinants of soil quality in mixed larch and birch forests: network structure and keystone taxa abundances. Front. Plant Sci. 16. doi: 10.3389/fpls.2025.1491038

PubMed Abstract | Crossref Full Text | Google Scholar

Guan, H. and Fan, J. (2020). Effects of vegetation restoration on soil quality in fragile karst ecosystems of southwest China. PeerJ 8, e9456. doi: 10.7717/peerj.9456

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, J., Feng, H., McNie, P., Liu, Q., Xu, X., Pan, C., et al. (2023). Species mixing improves soil properties and enzymatic activities in Chinese fir plantations: A meta-analysis. CATENA 220, 106723. doi: 10.1016/j.catena.2022.106723

Crossref Full Text | Google Scholar

Guo, J., Feng, H., Roberge, G., Feng, L., Pan, C., McNie, P., et al. (2022). The negative effect of Chinese fir (Cunninghamia lanceolata) monoculture plantations on soil physicochemical properties, microbial biomass, fungal communities, and enzymatic activities. For. Ecol. Manage. 519, 120297. doi: 10.1016/j.foreco.2022.120297

Crossref Full Text | Google Scholar

Guo, L., Sun, Z., Ouyang, Z., Han, D., and Li, F. (2017). A comparison of soil quality evaluation methods for Fluvisol along the lower Yellow River. CATENA 152, 135–143. doi: 10.1016/j.catena.2017.01.015

Crossref Full Text | Google Scholar

Guo, Y., Abdalla, M., Espenberg, M., Hastings, A., Hallett, P., and Smith, P. (2021). A systematic analysis and review of the impacts of afforestation on soil quality indicators as modified by climate zone, forest type and age. Sci. Total Environ. 757, 143824. doi: 10.1016/j.scitotenv.2020.143824

PubMed Abstract | Crossref Full Text | Google Scholar

He, H., Peng, M., Lu, W., Hou, Z., and Li, J. (2022). Commercial organic fertilizer substitution increases wheat yield by improving soil quality. Sci. Total Environ. 851, 158132. doi: 10.1016/j.scitotenv.2022.158132

PubMed Abstract | Crossref Full Text | Google Scholar

Hemati, Z., Selvalakshmi, S., Xia, S., and Yang, X. (2020). Identification of indicators: Monitoring the impacts of rubber plantations on soil quality in Xishuangbanna, Southwest China. Ecol. Indic. 116, 106491. doi: 10.1016/j.ecolind.2020.106491

Crossref Full Text | Google Scholar

Hou, X., Han, H., Tigabu, M., Cai, L., Meng, F., Liu, A., et al. (2019). Changes in soil physico-chemical properties following vegetation restoration mediate bacterial community composition and diversity in Changting, China. Ecol. Eng. 138, 171–179. doi: 10.1016/j.ecoleng.2019.07.031

Crossref Full Text | Google Scholar

Hou, L., Zhang, Y., Li, Z., Shao, G., Song, L., and Sun, Q. (2021). Comparison of soil properties, understory vegetation species diversities and soil microbial diversities between chinese fir plantation and close-to-natural forest. Forests 12, 632. doi: 10.3390/f12050632

Crossref Full Text | Google Scholar

Hua, F., Bruijnzeel, L. A., Meli, P., Martin, P. A., Zhang, J., Nakagawa, S., et al. (2022). The biodiversity and ecosystem service contributions and trade-offs of forest restoration approaches. Science 376, 839–844. doi: 10.1126/science.abl4649

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, Z., Selvalakshmi, S., Vasu, D., Liu, Q., Cheng, H., Guo, F., et al. (2018). Identification of indicators for evaluating and monitoring the effects of Chinese fir monoculture plantations on soil quality. Ecol. Indic. 93, 547–554. doi: 10.1016/j.ecolind.2018.05.034

Crossref Full Text | Google Scholar

Hui, G., Zhang, G., Zhao, Z., and Yang, A. (2019). Methods of forest structure research: a review. Curr. For. Rep. 5, 142–154. doi: 10.1007/s40725-019-00090-7

Crossref Full Text | Google Scholar

ISMAPRC (1996). Classification of type regions and fertility of cultivated land in China (NY/T 309-1993).

Google Scholar

Jin, H., Shi, D., Lou, Y. B., Zhang, J., Ye, Q., and Jiang, N. (2021). Evaluation of the quality of cultivated-layer soil based on different degrees of erosion in sloping farmland with purple soil in China. CATENA 198, 105048. doi: 10.1016/j.catena.2020.105048

Crossref Full Text | Google Scholar

Kweon, D. and Comeau, P. G. (2019). Relationships between tree survival, stand structure and age in trembling aspen dominated stands. For. Ecol. Manage. 438, 114–122. doi: 10.1016/j.foreco.2019.02.003

Crossref Full Text | Google Scholar

LaRue, E. A., Hardiman, B. S., Elliott, J. M., and Fei, S. (2019). Structural diversity as a predictor of ecosystem function. Environ. Res. Lett. 14, 114011. doi: 10.1088/1748-9326/ab49bb

Crossref Full Text | Google Scholar

Li, Y., Dong, S., Wen, L., Wang, X., and Wu, Y. (2013). Assessing the soil quality of alpine grasslands in the Qinghai-Tibetan Plateau using a modified soil quality index. Environ. Monit. Assess. 185, 8011–8022. doi: 10.1007/s10661-013-3151-1

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Li, M., Li, X., Liu, Z., Ming, A., Lan, H., et al. (2021). The abundance and structure of deadwood: A comparison of mixed and thinned chinese fir plantations. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.614695

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X., Liu, Y., Wu, G., Lie, Z., Sheng, H., Aguila, L. C. R., et al. (2024). Mixed plantations do not necessarily provide higher ecosystem multifunctionality than monoculture plantations. Sci. Total Environ. 914, 170156. doi: 10.1016/j.scitotenv.2024.170156

PubMed Abstract | Crossref Full Text | Google Scholar

Li, P., Shi, K., Wang, Y., Kong, D., Liu, T., Jiao, J., et al. (2019a). Soil quality assessment of wheat-maize cropping system with different productivities in China: Establishing a minimum data set. Soil Tillage Res. 190, 31–40. doi: 10.1016/j.still.2019.02.019

Crossref Full Text | Google Scholar

Li, X., Wang, D., Ren, Y., Wang, Z., and Zhou, Y. (2019b). Soil quality assessment of croplands in the black soil zone of Jilin Province, China: Establishing a minimum data set model. Ecol. Indic. 107, 105251. doi: 10.1016/j.ecolind.2019.03.028

Crossref Full Text | Google Scholar

Li, Y., Zeng, C., and Long, M. (2020). Variation of soil nutrients and bacterial community diversity of different land utilization types in Yangtze River Basin, Chongqing Municipality. PeerJ 8, e9386. doi: 10.7717/peerj.9386

PubMed Abstract | Crossref Full Text | Google Scholar

Liang, C., Liu, L., Zhang, Z., Ze, S., Ji, M., Li, Z., et al. (2022). Do mixed pinus yunnanensis plantations improve soil’s physicochemical properties and enzyme activities? Diversity 14, 214. doi: 10.3390/d14030214

Crossref Full Text | Google Scholar

Liu, W., Liu, M., Li, W., Zeng, F., and Qu, Y. (2016). Influence of ginseng cultivation under larch plantations on plant diversity and soil properties in Liaoning Province, Northeast China. J. Mt. Sci. 13, 1598–1608. doi: 10.1007/s11629-015-3753-x

Crossref Full Text | Google Scholar

Liu, D., Zhou, C., He, X., Zhang, X., Feng, L., and Zhang, H. (2022). The effect of stand density, biodiversity, and spatial structure on stand basal area increment in natural spruce-fir-broadleaf mixed forests. Forests 13, 162. doi: 10.3390/f13020162

Crossref Full Text | Google Scholar

Luo, Y. and Chen, H. Y. H. (2011). Competition, species interaction and ageing control tree mortality in boreal forests. J. Ecol. 99, 1470–1480. doi: 10.1111/j.1365-2745.2011.01882.x

Crossref Full Text | Google Scholar

Mamehpour, N., Rezapour, S., and Ghaemian, N. (2021). Quantitative assessment of soil quality indices for urban croplands in a calcareous semi-arid ecosystem. Geoderma 382, 114781. doi: 10.1016/j.geoderma.2020.114781

Crossref Full Text | Google Scholar

Návar, J. (2009). Allometric equations for tree species and carbon stocks for forests of northwestern Mexico. For. Ecol. Manage. 257, 427–434. doi: 10.1016/j.foreco.2008.09.028

Crossref Full Text | Google Scholar

Nosrati, K. (2013). Assessing soil quality indicator under different land use and soil erosion using multivariate statistical techniques. Environ. Monit. Assess. 185, 2895–2907. doi: 10.1007/s10661-012-2758-y

PubMed Abstract | Crossref Full Text | Google Scholar

Onoda, Y., Saluñga, J. B., Akutsu, K., Aiba, S., Yahara, T., and Anten, N. P. R. (2014). Trade-off between light interception efficiency and light use efficiency: implications for species coexistence in one-sided light competition. J. Ecol. 102, 167–175. doi: 10.1111/1365-2745.12184

Crossref Full Text | Google Scholar

Ouyang, S., Gou, M., Lei, P., Liu, Y., Chen, L., Deng, X., et al. (2023). Plant functional trait diversity and structural diversity co-underpin ecosystem multifunctionality in subtropical forests. For. Ecosyst. 10, 100093. doi: 10.1016/j.fecs.2023.100093

Crossref Full Text | Google Scholar

Paquette, A. and Messier, C. (2011). The effect of biodiversity on tree productivity: from temperate to boreal forests. Glob. Ecol. Biogeogr. 20, 170–180. doi: 10.1111/j.1466-8238.2010.00592.x

Crossref Full Text | Google Scholar

Põldveer, E., Potapov, A., Korjus, H., Kiviste, A., Stanturf, J. A., Arumäe, T., et al. (2021). The structural complexity index SCI is useful for quantifying structural diversity of Estonian hemiboreal forests. For. Ecol. Manage. 490, 119093. doi: 10.1016/j.foreco.2021.119093

Crossref Full Text | Google Scholar

Poorter, L., van der Sande, M. T., Arets, E. J. M. M., Ascarrunz, N., Enquist, B. J., Finegan, B., et al. (2017). Biodiversity and climate determine the functioning of Neotropical forests. Glob. Ecol. Biogeogr. 26, 1423–1434. doi: 10.1111/geb.12668

Crossref Full Text | Google Scholar

Qian, F., Yu, Y., Dong, X., and Gu, H. (2023). Soil quality evaluation based on a minimum data set (MDS)—A case study of tieling county, northeast China. Land 12, 1263. doi: 10.3390/land12061263

Crossref Full Text | Google Scholar

Rahmanipour, F., Marzaioli, R., Bahrami, H. A., Fereidouni, Z., and Bandarabadi, S. R. (2014). Assessment of soil quality indices in agricultural lands of Qazvin Province, Iran. Ecol. Indic. 40, 19–26. doi: 10.1016/j.ecolind.2013.12.003

Crossref Full Text | Google Scholar

Raiesi, F. and Beheshti, A. (2022). Evaluating forest soil quality after deforestation and loss of ecosystem services using network analysis and factor analysis techniques. CATENA 208, 105778. doi: 10.1016/j.catena.2021.105778

Crossref Full Text | Google Scholar

Rezaei, S. A., Gilkes, R. J., and Andrews, S. S. (2006). A minimum data set for assessing soil quality in rangelands. Geoderma 136, 229–234. doi: 10.1016/j.geoderma.2006.03.021

Crossref Full Text | Google Scholar

Schnabel, F., Schwarz, J. A., Dănescu, A., Fichtner, A., Nock, C. A., Bauhus, J., et al. (2019). Drivers of productivity and its temporal stability in a tropical tree diversity experiment. Glob. Change Biol. 25, 4257–4272. doi: 10.1111/gcb.14792

PubMed Abstract | Crossref Full Text | Google Scholar

Shao, G., Ai, J., Sun, Q., Hou, L., and Dong, Y. (2020). Soil quality assessment under different forest types in the Mount Tai, central Eastern China. Ecol. Indic. 115, 106439. doi: 10.1016/j.ecolind.2020.106439

Crossref Full Text | Google Scholar

Soares, A. A. V., Leite, H. G., Souza, A. L., Silva, S. R., Lourenço, H. M., and Forrester, D. I. (2016). Increasing stand structural heterogeneity reduces productivity in Brazilian Eucalyptus monoclonal stands. For. Ecol. Manage. 373, 26–32. doi: 10.1016/j.foreco.2016.04.035

Crossref Full Text | Google Scholar

Stewart, A., Reedyk, S., Franz, B., Fomradas, K., Hilliard, C., and Hall, S. (2011). A planning tool for design and location of vegetative buffers on watercourses in the Canadian prairies. J. Soil Water Conserv. 66, 97A–103A. doi: 10.2489/jswc.66.4.97A

Crossref Full Text | Google Scholar

Sun, K., Sun, H., Lu, G., Fang, L., Wan, Z., and Tan, Z. (2023). Study on the effects of stand density management of Chinese fir plantation in Northern China. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1130299

PubMed Abstract | Crossref Full Text | Google Scholar

Tetemke, B. A., Birhane, E., Rannestad, M. M., and Eid, T. (2021). Species diversity and stand structural diversity of woody plants predominantly determine aboveground carbon stock of a dry Afromontane forest in Northern Ethiopia. For. Ecol. Manage. 500, 119634. doi: 10.1016/j.foreco.2021.119634

Crossref Full Text | Google Scholar

Thurm, E. A. and Pretzsch, H. (2016). Improved productivity and modified tree morphology of mixed versus pure stands of European beech (Fagus sylvatica) and Douglas-fir (Pseudotsuga menziesii) with increasing precipitation and age. Ann. For. Sci. 73, 1047–1061. doi: 10.1007/s13595-016-0588-8

Crossref Full Text | Google Scholar

Ullah, F., Gilani, H., Sanaei, A., Hussain, K., and Ali, A. (2021). Stand structure determines aboveground biomass across temperate forest types and species mixture along a local-scale elevational gradient. For. Ecol. Manage. 486, 118984. doi: 10.1016/j.foreco.2021.118984

Crossref Full Text | Google Scholar

Wan, P., Zhou, Z., Yuan, Z., Wei, H., Huang, F., Li, Z., et al. (2024). Fungal community composition changes and reduced bacterial diversity drive improvements in the soil quality index during arable land restoration. Environ. Res. 244, 117931. doi: 10.1016/j.envres.2023.117931

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Z., Zhang, X., Chhin, S., Zhang, J., and Duan, A. (2021). Disentangling the effects of stand and climatic variables on forest productivity of Chinese fir plantations in subtropical China using a random forest algorithm. Agric. For. Meteorol. 304–305, 108412. doi: 10.1016/j.agrformet.2021.108412

Crossref Full Text | Google Scholar

Ye, Y., Sun, X., Zhao, J., Wang, M., and Guan, Q. (2022). Establishing a soil quality index to assess the effect of thinning on soil quality in a Chinese fir plantation. Eur. J. For. Res. 141, 999–1009. doi: 10.1007/s10342-022-01471-7

Crossref Full Text | Google Scholar

Yu, P., Liu, J., Tang, H., Sun, X., Liu, S., Tang, X., et al. (2023). Establishing a soil quality index to evaluate soil quality after afforestation in a karst region of Southwest China. CATENA 230, 107237. doi: 10.1016/j.catena.2023.107237

Crossref Full Text | Google Scholar

Yu, P., Liu, S., Zhang, L., Li, Q., and Zhou, D. (2018). Selecting the minimum data set and quantitative soil quality indexing of alkaline soils under different land uses in northeastern China. Sci. Total Environ. 616–617, 564–571. doi: 10.1016/j.scitotenv.2017.10.301

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, Z., Wang, S., Ali, A., Gazol, A., Ruiz-Benito, P., Wang, X., et al. (2018). Aboveground carbon storage is driven by functional trait composition and stand structural attributes rather than biodiversity in temperate mixed forests recovering from disturbances. Ann. For. Sci. 75, 67. doi: 10.1007/s13595-018-0745-3

Crossref Full Text | Google Scholar

Zhang, X., Cao, Q. V., Duan, A., and Zhang, J. (2017). Modeling tree mortality in relation to climate, initial planting density, and competition in Chinese fir plantations using a Bayesian logistic multilevel method. Can. J. For. Res. 47, 1278–1285. doi: 10.1139/cjfr-2017-0215

Crossref Full Text | Google Scholar

Zhang, Y. and Chen, H. Y. H. (2015). Individual size inequality links forest diversity and above-ground biomass. J. Ecol. 103, 1245–1252. doi: 10.1111/1365-2745.12425

Crossref Full Text | Google Scholar

Zhang, X., Wang, Z., Chhin, S., Wang, H., Duan, A., and Zhang, J. (2020). Relative contributions of competition, stand structure, age, and climate factors to tree mortality of Chinese fir plantations: Long-term spacing trials in southern China. For. Ecol. Manage. 465, 118103. doi: 10.1016/j.foreco.2020.118103

Crossref Full Text | Google Scholar

Zhang, C., Xue, S., Liu, G.-B., and Song, Z.-L. (2011). A comparison of soil qualities of different revegetation types in the Loess Plateau, China. Plant Soil 347, 163–178. doi: 10.1007/s11104-011-0836-5

Crossref Full Text | Google Scholar

Zhou, L., Sun, Y., Saeed, S., Zhang, B., and Luo, M. (2020). The difference of soil properties between pure and mixed Chinese fir (Cunninghamia lanceolata) plantations depends on tree species. Glob. Ecol. Conserv. 22, e01009. doi: 10.1016/j.gecco.2020.e01009

Crossref Full Text | Google Scholar

Keywords: soil quality index, stand structure, soil properties, minimum data set, net primary productivity (NPP)

Citation: Fu Z, Xiao Y, Liu S, Xu H, Wang Y and Zhu H (2025) Soil quality enhancement drives tree growth and broadleaf dominance in fir-broadleaf mixed plantations. Front. Plant Sci. 16:1705626. doi: 10.3389/fpls.2025.1705626

Received: 15 September 2025; Accepted: 31 October 2025;
Published: 19 November 2025.

Edited by:

Yuguo Liu, Chinese Academy of Forestry, China

Reviewed by:

Ronghua Li, South China Agricultural University, China
Ziqing Zhao, Chinese Academy of Forestry, China

Copyright © 2025 Fu, Xiao, Liu, Xu, Wang and Zhu. 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: Yihua Xiao, anh4aWFveWlodWFAMTI2LmNvbQ==

Disclaimer: 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.