Elucidating the correlation between soil microbial communities and forest productivity is the focus of research in the field of forest ecology. Nonetheless, the relationship between stand age, soil quality, soil microorganisms, and their combined influence on productivity is still unclear. In this study, five development stages (14, 25, 31, 39, and >80 years) of larch (Larix principis-rupprechtii) forests were investigated in Inner Mongolia and Shanxi provinces of China. We evaluated soil quality using the Integrated Soil Quality Index (SQI) and analyzed changes in bacterial and fungal communities using high-throughput sequencing. Regression models were also established to examine the impacts of stand age, microbial diversity, and SQI on productivity. The findings revealed an ascending trend in soil organic matter (SOM), total nitrogen (TN), total phosphorus (TP), available potassium (AK), and SQI in 14, 25, 31, and 39-year-old stands. The abundance of oligotrophic bacteria Acidobacteria exhibited a gradual decline with increasing forest age, whereas copiotroph bacteria Proteobacteria displayed a progressive increase. Stands older than 80 years exhibited a higher abundance of both the saprophytic fungus Ascomycota and mycorrhizal fungus Basidiomycota. Forest age had a significant impact on microbial diversity, particularly in terms of bacterial diversity, impacting both α and β diversity. The soil bacterial community structure was influenced by AK, SOM, TN, TP, and pH. Conversely, the fungal community structure was regulated by crucial factors including SOM, TN, TP, TK, AK, and pH. Fungal diversity demonstrated a significant and positive correlation with the basal area increment (BAI) of larch. Furthermore, microbial diversity accounted for 23.6% of the variation in BAI. In summary, the findings implied a robust association between microbial composition, diversity, and soil chemical properties throughout the chronosequence of larch forests. These factors collectively played a crucial role in influencing the productivity of larch forest.

Understanding the factors influencing tree productivity is central to forest ecology. However, the relative contributions of neighborhood interactions, tree species diversity, and tree size to larch (Larix principis-rupprechtii) productivity require further study. Three plots in the Guandi Mountains, Shanxi Province, were set up for each of the following forest types: natural pure larch forest (PL), mixed larch and birch (Betula platyphylla) forest (LB), and mixed larch and spruce (Picea asperata) forest (LS). Based on the tree size-stratified sampling method, a total of 318 tree core samples were collected. A linear mixed model was used to analyze the effects of tree size, dominance, mixing, and neighborhood competition on larch productivity. Birch and spruce promoted larch growth at the stand and individual tree levels, and birch exhibited a more significant facilitating effect. Intraspecific competition was the main factor affecting larch growth. When the intensity of competition among trees was low, the basal area increment (BAI) of larch in the mixed forests was higher than that in the pure forest. However, with increasing competition, the BAI of larch was lower in the mixed forests than in the pure forest. Factors including tree size, dominance, and mingling were positively correlated with the BAI of larch. With increasing tree size, the BAI of larch was higher in the mixed forests than in the pure forest and higher in LB than in LS. When the dominance was less than 0.5, the BAI of larch was higher in the pure forest than in the mixed forests and higher in LS than in LB. With increasing dominance, the BAI of larch was higher in the mixed forests than in the pure forest. The BAI of larch increased with an increasing mixing degree in the mixed forests, and the increasing trend of BAI was larger in LB than in LS. Larch productivity was influenced mainly by neighborhood interactions and stand structure. Improving neighborhood tree diversity and increasing the large tree proportion and dominance of larch will be helpful for improving larch productivity in mixed forests.

The size–growth relationship (SGR) quantifies growth partitioning among different sized trees in a stand and helps to elucidate stand growth dynamics during stand development. SGR strongly correlates with stand density, stand age, and drought severity. This study focused on larch (Larix principis-rupprechtii (Mayr)) forests with different stand ages (17–19 years, 20–29 years, 30–39 years, and 40–46 years) and stand densities (300–1400 trees·ha⁻¹ and 1400–3300 trees·ha⁻¹) as the research subjects. Employing a linear mixed model, we aimed to quantify the effects of stand density, stand age, and drought index on SGR. The results revealed that the Gini coefficient had a significant variation between high-density and low-density larch stands (p < 0.01). Stand age, stand density, and drought index exerted varying degrees of influence on the growth dominance coefficient (GD) and SGR of larch forests. In stands aged less than 39 years, a pronounced growth advantage of large trees over small trees was evident (typically SGR > 1 and GD > 0), indicating a positive growth dominance stage where large trees dominated. Compared to high-density stands, low-density stands exhibited a notably greater positive growth dominance. After 40 years, the growth stage transitioned to a relatively symmetric stage (SGR approximately equal to 1), with a discernible shift towards a reverse growth dominance stage (GD < 0). Compared to stand density and stand age, the drought index had a more influential effect on SGR. As drought severity increased, SGR increased, amplifying the growth advantage of large trees over small ones. The findings underscored the significance of adjusting stand density and optimizing tree size structure to enhance larch resilience against the warming and drying effects.

Larix principis-rupprechtii Mayr (larch) is one of the main afforestation and timber production species used in North China. Climate change has led to a change in its suitable distribution and growth. However, the impact of climate change on its growth suitability is not clear. In this study, using forest resource inventory data and spatially continuous environmental factor data (temperature, precipitation, topography, and soil) in Hebei and Shanxi Provinces, China, the random forest model (RF) was used to simulate the larch site index (SI) and growth suitability under three shared socioeconomic pathways (SSPs: SSP1-2.6, SSP2-4.5, and SSP5-8.5) for the current and future (2021–2040, 2041–2060 and 2080–2100). The results revealed that (1) RF had excellent performance in predicting the regional SI (R² = 0.73, MAE = 0.93 m, RMSE = 1.35 m); (2) the main factors affecting the productivity of larch were the mean temperature of the warmest quarter (BIO10), elevation (ELEV), mean diurnal range (BIO2), and annual precipitation (BIO12); and (3) larch currently had a higher SI in the Bashang areas and in the high-altitude mountains. The areas characterized as unsuitable, poorly suitable, moderately suitable, and highly suitable accounted for 15.45%, 42.12%, 31.94%, and 10.49% of the total area, respectively. (4) Future climate warming had an obvious inhibitory effect on the SI, and the effect strengthened with increasing radiation intensity and year. (5) The moderately suitable and highly suitable areas of larch growth showed a downward trend under future climate scenarios. By the end of this century, the suitable growth areas would decrease by 14.14% under SSP1-2.6, 15.17% under SSP2-4.5, and 19.35% under SSP5-8.5. The results revealed the impact of climate change on larch growth suitability, which can provide a scientific basis for larch forest management.

Forests can simultaneously provide a variety of ecosystem functions and services (ecosystem multifunctionality, EMF). Different stand types, influenced by biotic and abiotic factors, play a key role in determining EMF. To clarify the impact of stand type, as well as biotic and abiotic factors, on EMF, this study quantified EMF across three stand types: Betula platyphylla pure forest (BP), B. platyphylla–Larix principis-rupprechtii mixed forest (BL), and L. principis-rupprechtii pure forest (LP). The multiple-threshold approach was employed to quantify EMF, while structural equation modeling was used to analyze the primary factors influencing EMF. The results indicated the following: (1) BL had higher stand productivity than both BP and LP; (2) BL exhibited significantly higher functional diversity and soil fertility maintenance compared to LP, with no significant difference between BP and BL; (3) BP demonstrated a significantly stronger nutrient cycling function than LP, with no significant difference between LP and BL; (4) the ranking of EMF at all threshold levels was BL > BP > LP; (5) soil was an effective driver of EMF across all threshold levels; and (6) both the niche complementarity effect and the mass ratio effect jointly drove EMF at the low threshold (<50%), with the influence of both effects diminishing as the threshold increased. This study enhances our understanding of the key drivers of EMF in forest ecosystems and provides valuable insights for informing multifunctional forest management practices.