Does stand density affect understory vegetation and soil properties of differently aged Robinia pseudoacacia plantations?

Robinia pseudoacacia is the most important tree species used for afforestation on the Chinese Loess Plateau. However, most R. pseudoacacia plantations face environmental problems due to inappropriate initial densities. The sustainable development of R. pseudoacacia plantations is challenging due to significant knowledge gaps regarding the appropriate stand densities for maximizing ecological benefits. Herein, R. pseudoacacia plantations of three age classes—young, middle-aged, and mature—were sampled in the loess hilly-gully region of northern Shaanxi Province, China. Based on stand density, R. pseudoacacia plantations in each age class were categorized as low (<1000 trees ha−1), medium (1000–1500 trees ha−1), and high (>1500 trees ha−1) density. We analyzed the effects of stand density on soil physicochemical properties, understory vegetation community structure, and ecosystem multifunctionality. In young stands, a medium stand density significantly increased the Shannon-Wiener index, species richness, and above-ground biomass of herbaceous communities. In middle-aged stands, stand density had no significant effects on understory diversity, whereas medium-density stands showed increased understory above-ground biomass. In mature stands, low-density stands significantly increased understory above-ground biomass and species richness of herbaceous communities. In all age classes, stand density did not significantly affect shrubby community composition, but herbaceous communities shifted compositionally towards shade-tolerant and drought-tolerant species as stand density increased. In the soils of young and middle-aged stands, organic carbon, ammonium nitrogen, nitrate nitrogen, available phosphorus, and water content were highest in medium-density stands, whereas in mature stands, these variables were highest in low-density plantations. In young, middle-aged, and mature stands, the ecosystem multifunctionality indices were highest at densities of 1000–1500, 1000–1500, and < 1000 trees ha−1, respectively. We concluded that the effects of stand density on ecosystem functions in R. pseudoacacia plantations varied with stand age. Our findings highlight a need to implement age-based stand density management strategies for R. pseudoacacia plantations to improve their ecological benefits.