Intensified Interspecific Competition for Water after Afforestation with Robinia pseudoacacia into a Native Shrubland in the Taihang Mountains, Northern China

Understanding how soil water source is used spatiotemporally by tree species and if native species can successfully coexist with introduced species is crucial for selecting species for afforestation. In the rocky mountainous areas of the Taihang Mountains, alien Robinia pseudoacacia L. has been widely afforested into the native shrublands dominated by Ziziphus jujuba Mill var. spinosa and Vitex negundo L. var. heterophylla to improve forest coverage and soil nutrients. However, little is known about the water relation among species, especially seasonal water use sources in different microsites. We selected the soil and plant xylem samples of two opposite microtopographic sites (ridge and valley) monthly in the growth season to analyze isotope composition. The proportions of water sources were quantified by the MixSIAR model and compared pairwise between species, microsites and seasons. We found that deep subsoil water at a depth of 40–50 cm contributed up to 50% of the total water uptake for R. pseudoacacia and Z. jujuba in the growing season, indicating that they stably used deeper soil water and had intense water competition. However, V. negundo had a more flexible water use strategy, which derived more than 50% of the total water uptake from the soil layer of 0–10 cm in the rainy season, but majorly captured soil water at a depth of 30–50 cm in the dry season. Therefore, high niche overlaps were shown in V. negundo with the other two species in the dry season, but niche segregation was seen in the rainy season. The microtopographic sites did not shift the seasonal dynamic of the water source use patterns of the three studied species, but the water use niche overlap was higher in the valley than in the ridge. Taken together, the introduced species R. pseudoacacia intensified water competition with the native semi-arbor species Z. jujuba, but it could commonly coexist with the native shrub species V. negundo. Therefore, our study on seasonal water use sources in different microsites provides insight into species interaction and site selection for R. pseudoacacia afforestation in the native shrub community in rocky mountainous areas. It is better to plant R. pseudoacacia in the shrubland in the valley so as to avoid intense water competition and control soil erosion.


Introduction
Afforestation is one of the main measures for vegetation restoration [1], and several key national ecological restoration programs have been undertaken in China with the long-running aim of fighting against land degradation, and preventing sand and dust is helpful for understanding and predicting the structure and function of ecosystems in water-stress conditions [49,50].
A large number of studies on water use sources have been conducted in arid and semiarid regions, and coastal sand dune systems [40,43,[51][52][53][54][55]; however, little is known about the plant water sources of the restored and reconstructed forests in the hilly areas of the Taihang Mountains [21]. The Taihang Mountains Afforestation Project has been primarily carried out for the purpose of conserving water and soil and maintaining regional ecological security [2,56]. As it lies in an inland temperate zone affected by monsoons, vegetation restoration often suffers from seasonal drought, which affects the growth and survival of afforested species, and even influences the succession of restored vegetation [7]. Ziziphus jujuba and Vitex negundo are the dominant woody species of native shrub communities, as well as the pioneer species in the water limited areas with strong resistance to water stress in the infertile hilly land [7,56]. R. pseudoacacia is an introduced tree species for afforestation dominated by shrubs to improve soil nutrients in the Taihang Mountains [14,57]. These three woody plants are the common species that constitute hilly plant communities, representing a typical afforestation mode of the ecological restoration program. For this reason, we used the stable isotope method to study seasonal water sources and niche segregation of three co-occurring woody species in two typically microtopographic sites, i.e., in the ridge and valley, and explore whether R. pseudoacacia can stably coexist with the other two species after afforestation. Our main objective was to address whether the seasonal variation of water source supply and environmental heterogeneity could affect plant water use strategies and their interspecific relationship. We hypothesized that: (1) interspecific differentiation of plant water source use would present among the co-occurring species with different lifeforms; (2) plant water source use patterns changed seasonally and spatially with water distribution in the soil profiles, which influenced niche differentiation of the studied species; (3) the heterogeneity of soil water conditions in the microtopographic sites affects the proportion of plant water source use.

Study Sites and Target Plant Species
This study was conducted at Taihang Hilly Ecosystem Experimental Station (37 • Figure 1). The station is characterized by a semiarid continental monsoon climate, with annual mean air temperature of 13.2 • C, coldest monthly temperature of −1.6 • C in January, and warmest of 26.3 • C in July [6]. The mean annual precipitation is 570 mm, with over 65% of the precipitation falling between July and September [58,59]. Mean annual evapotranspiration is 1934 mm, about 3.4 times of the precipitation. The mountainous cinnamon soil, with a lot of pebbles and coarse soil textures, is developed mainly from the weathering of granitic gneiss, partly by limestone and shale. The soil condition is well drained and generally nutrient poor, with a weak capacity for retaining fertilizer and water. Before the 1980s, large areas of mountain forest were logged and reconstructed as secondary forest and plantations through vegetation restoration [6,59].
Considering the topographical difference in soil water conditions, we conducted this research at two microtopographic sites, in the ridge and valley, at a distance of 500 m. At each site, one HOBO micrometeorological station was installed to automatically record meteorological data at a frequency of every half hour. The data showed that there were no differences in temperature and precipitation except soil moisture between the two sites during the experiment period. The growth status of mixed forests of R. pseudoacacia with Z. jujuba and V. negundo was investigated at both microsites (Table 1). R. pseudoacacia grew better in the valley than in the ridge, whereas Z. jujuba favored growth in the ridge. No difference was found for the growth of V. negundo. Z. jujuba and V. negundo was investigated at both microsites (Table 1). R. pseudoacacia grew better in the valley than in the ridge, whereas Z. jujuba favored growth in the ridge. No difference was found for the growth of V. negundo.

Sample Collection and Stable Isotope Analyses
Five plots (10 m × 10 m), separated by at least 100 m, were chosen for sample collection at each of the two microtopographic sites. In each plot, a sampling point with three target species growing together and similar coverage to the corresponding plot was used for sampling. The samples were collected on two consecutive sunny days of each month,  The lowercase letters (a, b) represented significant differences of tree height and diameter at breast height (DBH) for the same species between the ridge and the valley sites (p < 0.05).

Sample Collection and Stable Isotope Analyses
Five plots (10 m × 10 m), separated by at least 100 m, were chosen for sample collection at each of the two microtopographic sites. In each plot, a sampling point with three target species growing together and similar coverage to the corresponding plot was used for sampling. The samples were collected on two consecutive sunny days of each month, once a month, from May to September 2017. For each sampling time, samples were collected in 10 plots (five plots × two microtopographic sites). Therefore, five replicate sampling points were used for each sampling site and month. A total of 8-10 segments of non-photosynthetic and lignified twigs from each target species in each sampling point were clipped and mixed as one sample from the sunward canopy of each species. To maximally reduce the risk of stem water evaporation, the peeled xylem parts of the stem segments were immediately put into screw-cap vials, and then sealed with parafilm [60]. Concurrently, soil samples were collected at the same locations under the canopy of plant sampling. In each point, soil samples were collected with soil auger in 5 layers with an interval of 10 cm and then divided into 2 parts. One part was sealed in screw-cap vials for water extraction and isotopic analyses; the other part was stored in aluminum boxes to measure gravimetric soil water content. All plant stems and soil samples were immediately placed into a portable cooler and then stored in a refrigerator at −20 • C until water extraction and isotopic analysis.
Water from plant xylem and soil samples was extracted by a cryogenic vacuum distillation system (LI-2100, Lica United, Beijing, China) and then filtered using a 0.2-µm filter before stable isotopic analysis. The stable isotope composition of hydrogen (δ 2 H) and oxygen (δ 18 O) in water were determined by a laser absorption water vapor isotope analyzer (Picarro-i2120, Santa Clara, CA, USA) at the Key Laboratory of Agricultural Water Resources, Chinese Academy of Science [61]. The isotope ratio of δ 2 H or δ 18 O expressed in standard δ-notation as per mil (‰) was calculated as where R sample and R standard are the isotope ratio ( 2 H/ 1 H, 18

Quantifying the Proportion of Water Use Sources
A Bayesian stable isotope mixing model (MixSIAR) for R was used to calculate the relative contribution of water use sources to plant growth [62][63][64]. The raw xylem isotope values (δ 2 H and δ 18 O) from three co-occurring species in two microtopographic sites were used as mixture data. The means and standard errors of isotope ratios (δ 2 H and δ 18 O) from five soil layers were used as source data, which had no concentration dependence but varied by species and sites. The discrimination data were set as zero because no fractionation occurred in δ 2 H and δ 18 O during plant water uptake [65]. Fixed effects (species and microtopographic sites), where species were nested within microtopographic sites, were conducted in all analyses. Markov Chain Monte Carlo (MCMC) parameters were specified as "very long" (number of chains = 3, chain length = 1,000,000, burn-in = 500,000, thin = 500), aiming to converge on the posterior distributions for all variables in the model [66]. Before accepting any of the MixSIAR results, we ensured that the model had converged. Gelman-Rubin and Geweke diagnostic tests were conducted to ascertain whether the model was close to convergence [67]. We chose the "Residual * Process" option of error structure and set it to "uninformative prior" in the model [63]. Finally, we ran the MixSIAR model and obtained the process output. The median values were presented as the proportional contribution of water sources to plant species.

Species Niche Overlap of Water Use
Pairwise niche overlap of co-occurring species was calculated using the index of proportional similarity (PS) based on Trogisch et al. [68] and Hoekstra et al. [69]: where p 1i and p 2i represented the proportional contribution (in %) to water uptake of the vertical soil depth intervals within species. The value of PS is minimum (0) when there is no overlap of water use source between two species, and is maximum (1) when there is complete overlap.

Statistical Analysis
Statistical analyses were performed using SPSS 17.0 (SPSS Inc., Chicago, IL, USA). One-way ANOVA was used to calculate the means and standard errors and test significant differences of the data among species and soil layers in different seasonal periods and at different sites. Multi-factor ANOVA was conducted to test the interaction effects of site, season, and soil depth on soil water content (SWC) and soil water isotopic values, and also quantify the effects of site, season, and species on plant xylem water isotopic values. We used the least-significant difference (LSD) test (at level of α = 0.05) nested in multi-factor ANOVA to determine whether there was a significant difference in the periods of season, site, soil depth, or species. All pictures were graphed with OriginPro 2016 (Origin Lab Corp., Northampton, MA, USA).

Variation in Soil Water Content
The seasonal and vertical variations of SWC showed similar trends in two microtopographic sites during the studied season, while the mean values of SWC in different soil depths were significantly higher in the valley than in the ridge ( Figure 2, Table 2). The main effects of season and soil depth presented significant differences in SWC ( Table 2). During the growing season, there was a summer peak of SWC in August and a spring lowest in June. In the dry season, the SWC decreased slowly with soil depth and showed a significant layer difference between 0-10 cm and 40-50 cm. In the rainy season, the SWC decreased more steeply with soil depth, with significantly lower SWC at the soil layer below 30 cm than that at 10 cm ( Figure 2). The SWC of each soil layer continuously decreased after August but showed no significant difference below the depth of 30 cm in May and September.
One-way ANOVA was used to calculate the means and standard errors and test signifi-cant differences of the data among species and soil layers in different seasonal periods and at different sites. Multi-factor ANOVA was conducted to test the interaction effects of site, season, and soil depth on soil water content (SWC) and soil water isotopic values, and also quantify the effects of site, season, and species on plant xylem water isotopic values. We used the least-significant difference (LSD) test (at level of α = 0.05) nested in multi-factor ANOVA to determine whether there was a significant difference in the periods of season, site, soil depth, or species. All pictures were graphed with OriginPro 2016 (Origin Lab Corp., Northampton, MA, USA).

Variation in Soil Water Content
The seasonal and vertical variations of SWC showed similar trends in two microtopographic sites during the studied season, while the mean values of SWC in different soil depths were significantly higher in the valley than in the ridge ( Figure 2, Table 2). The main effects of season and soil depth presented significant differences in SWC (Table 2). During the growing season, there was a summer peak of SWC in August and a spring lowest in June. In the dry season, the SWC decreased slowly with soil depth and showed a significant layer difference between 0-10 cm and 40-50 cm. In the rainy season, the SWC decreased more steeply with soil depth, with significantly lower SWC at the soil layer below 30 cm than that at 10 cm ( Figure 2). The SWC of each soil layer continuously decreased after August but showed no significant difference below the depth of 30 cm in May and September.

Variation in Soil Water Isotopic Values
The isotopic values of δ 18 O and δ 2 H of soil water both varied significantly at different soil layers, in different periods of season, between microtopographical sites, and with interaction among them (Table 2). Spatially, the δ 18 O and δ 2 H values decreased with increasing soil depths because of reduced evaporation and isotope fractionation. However, there was an exception in the dry season (May and June) with more negative isotope values of shallower soil layers above 20 cm and fewer changes of subsoil layers below 20 cm in the valley (Figure 3). Isotopic values of subsoil layers (30-50cm) between the two sites were significantly different in May, June, and July ( Figure 3). Seasonally, the isotopic values of δ 18 O and δ 2 H fluctuated more substantially over time in the shallower soil layer (0-10 cm) than in the other subsoil layers due to the interaction of precipitation and evaporation. The isotopic values of soil water were less negative in the dry season (May to June) than in the rainy season (July to September), showing a more enriched heavy isotope due to strong evaporation and less depletion due to limited precipitation in the dry season ( Figure 3).

Variation in Plant Xylem Water Isotopic Values
The isotopic values of plant xylem water showed differences among species and periods of season, but no differences between microtopographic sites ( Table 3). The average δ 18 O and δ 2 H values of R. pseudoacacia and Z. jujuba were more negative than that of V.

Variation in Plant Xylem Water Isotopic Values
The isotopic values of plant xylem water showed differences among species and periods of season, but no differences between microtopographic sites ( Table 3). The average δ 18 O and δ 2 H values of R. pseudoacacia and Z. jujuba were more negative than that of V. negundo, but there was no difference between R. pseudoacacia and Z. jujuba. Differences in the isotopic values of species were not noticeable in the dry season, but were clearly intensified during the rainy season ( Figure 4). The seasonal plant xylem δ 2 H and δ 18 O values of V. negundo increased from the dry season to the rainy season. However, those of R. pseudoacacia varied significantly only in the dry season, but there was no difference in the rainy season, with a minimum in May. In comparison, those of Z. jujuba changed little between seasons (Figure 4).

Proportion of Plant Water Source Utilization
On the whole, there were more major sources of water uptake in the humid valley than in the relatively drier ridge ( Figure 5).
During the dry season, no marked species differences occurred in the water source of the same soil layer ( Figure 5). However, water uptake proportion of different soil depths changed with seasons and sites. In May, all studied species used a higher proportion, up  (Figure 5d-f). In June, these species utilized water evenly at each layer below a depth of 10 cm in the ridge (Figure 5a-c), but were convergent in taking up water from the 10-20 cm soil layer in the valley (Figure 5d-f).

Interspecific Niche Overlap
The indices of water use niche overlap in specific pairwise species among the studied species were higher in the dry season than in the rainy season, and better water availability strengthened niche overlap between species in the valley in contrast to the water-limited ridge site. For example, more species niche overlap of water source use was found between R. pseudoacacia and Z. jujuba, with PS values over 90% in the dry season and up to 80% in the rainy season, and higher overlap in the valley than in the ridge. However, the PS values of V. negundo with R. pseudoacacia and Z. jujuba varied in different periods of the growing season. V. negundo showed high niche overlap with R. pseudoacacia and Z. jujuba in the dry season, but shifted to niche segregation with R. pseudoacacia and Z. jujuba in the rainy season (Table 4). During the rainy season, V. negundo mainly derived up to 60% and 40% of the water from surface soil above 10 cm in July and August, respectively, and showed no significant difference between sites (Figure 5a,d). However, similar to the dry season, R. pseudoacacia and Z. jujuba both continued to utilize water stably from the soil layer of 40-50 cm in July and August (Figure 5b,c,e,f). In September, with the decrease in soil water, V. negundo started to absorb deeper soil water at a depth of 30-40 cm whereas R. pseudoacacia and Z. jujuba explored shallower soil water of 20-30 cm ( Figure 5).

Interspecific Niche Overlap
The indices of water use niche overlap in specific pairwise species among the studied species were higher in the dry season than in the rainy season, and better water availability strengthened niche overlap between species in the valley in contrast to the water-limited ridge site. For example, more species niche overlap of water source use was found between R. pseudoacacia and Z. jujuba, with PS values over 90% in the dry season and up to 80% in the rainy season, and higher overlap in the valley than in the ridge. However, the PS values of V. negundo with R. pseudoacacia and Z. jujuba varied in different periods of the growing season. V. negundo showed high niche overlap with R. pseudoacacia and Z. jujuba in the dry season, but shifted to niche segregation with R. pseudoacacia and Z. jujuba in the rainy season (Table 4). Table 4. Niche overlap between species in the dry season (i = 10, 5 soil layers × 2 months), in the rainy season (i = 15, 5 soil layers × 3 months) and in the whole season (i = 25, 5 soil layers × 5 months).

Soil Water Content and Its Isotopic Composition
The isotopic composition of soil water is mainly affected by evaporation and osmotic processes, as well as the mixing of new and old rainwater [70,71]. In our study, the δ 2 H and δ 18 O values of soil water presented more negatively in the rainy season than in the dry season, indicating isotopic depletion in the rainy season. Our results agree with Deng et al. [72] who found that more isotopic composition enrichment occurred in the dry season than in the wet season. In soil strata, the soil water isotopic values of δ 2 H and δ 18 O were higher at surface (0-10 cm) than other soil layers because water evaporation strongly occurs in the surface soil [30,70,73]. Additionally, due to water transfer and accumulation in lowlying areas through surface runoff and infiltration after rainfall, and being less subject to wind, which greatly reduces the evaporation of soil water [56], SWC was significantly higher in the valley than in the ridge, especially in soil layers of 0-30 cm. However, even though the SWC is higher in the valley, the isotopic values of subsoil layers (30-50 cm) were normally enriched compared to the ridge during the dry season ( Figure 3). Thus, we speculate that this result may be due to occasionally low-intensity rainfall, which could not easily infiltrate soil layer below 30 cm, or the greater tree height of R. pseudoacacia creating a different canopy density which influenced the isotopic composition [30].

Seasonal Variation of Water Source Use in Different Species
The seasonal variation of the isotopic composition of water in xylem can reflect the difference of water sources and water use strategies among species [74]. V. negundo used a water source majorly from a shallow (0-10 cm) soil layer in July and August, yet shifted to use soil water below 20 cm in the dry season. The result is consistent with the study of Zhao et al. [75], who found that the main water use source of the shrub V. negundo was from the soil layers above 30 cm. In fact, the plasticity of soil water use in different periods of the growing season is closely associated with plant root traits and root distribution in the soil profile. Many plants in arid and semi-arid ecosystems possess a functionally dimorphic root system, which has a zone of lateral roots to predominantly acquire water from surface soil during wet seasons, while deriving water through the tap root from deep soil layers when surface water is unavailable in the dry seasons [65,76]. The dimorphic root species normally switch to absorbing water mainly from shallow to deep soil layers in the dry season [5,65]. Accordingly, plants with dimorphic root systems are more competitive than those solely depending on surface soil water in water-limited environments [36]. As studies have reported that V. negundo has a dimorphic root, with a set of lateral roots radiating from the main root crown and one or more deeply penetrating tap (sinker) roots [66], we also observed in our sites that the extensive network of lateral roots was predominantly concentrated in the top 30 cm soil layer with the highest root surface area, allowing V. negundo to monopolize the shallow moisture. The flexible water use strategy explains the widespread distribution of the shrub V. negundo in semiarid northern China. R. pseudoacacia and Z. jujuba both have a coarse, deeper tap root to fix plants and a series of lateral fine roots growing from a deep tap root. The deeper root systems of these two tree species enable their capacity to use water in the subsoil at 40-50 cm regardless of dry or rainy season. Besides, R. pseudoacacia and Z. jujuba can also take up water from soil layers below 10-20 cm in the dry season (particularly evident in June, Figure 5e,f) because, on the one hand, water was redistributed by plant roots lifting deep soil water to the shallower layer and using it [77][78][79][80]; on the other hand, water acquisition of fine roots in a shallow layer increased as a result of less water or low frequency of a small rainfall event in the dry season penetrating the deep soil. For R. pseudoacacia and Z. jujuba, there were no significant differences of isotopic values and seasonal variation, which implied that these two plants (arbor and semi-arbor) can stably use water from subsoil till to a depth of 40-50 cm. The similar water source use pattern could lead to water competition between them [21,22,81,82].
Normally, the proportion of soil water absorbed by plants increases with the rise of SWC [48,52,83], which can explain our result that the proportions of the major sources of water uptake were higher in the valley than in the ridge ( Figure 5). Compared to the differentiated root distribution of R. pseudoacacia (relatively shallower) and Z. jujuba (deeper) in a drier habitat, the root distribution of Z. jujuba became shallower due to improved soil and water conditions, and overlapped with R. pseudoacacia, which also exacerbated water competition in the valley. However, high-heterogeneity SWC between microtopographic sites did not change the seasonal dynamics of water source use patterns of the three species, perhaps because the difference in SWC did not reach the threshold of changing the species' water use strategy [29,84,85].

Niche Overlap in Water Source Use and Interspecies Relationship
According to niche theory, the niche segregation of water resource use can promote species coexistence through better allocation and utilization of common resources [28,86,87]. The water use strategy of species is the result of environmentally induced phenotypic plasticity, adaptive evolution or both [88]. V. negundo and Z. jujuba can stably coexist in the native communities of rocky mountainous areas due to water use niche segregation. However, when R. pseudoacacia afforested into the native shrub community, it intensified the water use overlap with shrubs, especially Z. jujuba ( Table 4). Differentiation of water use strategy exists in shrub species V. negundo and arbor/semi-arbor species of R. pseudoacacia and Z. jujuba due to the different plant growth form and root distribution feature. With the seasonal variation of water source use of V. negundo, the status of water use niche overlap changed. Water use sources among three species showed intensive overlap in the dry season but relieved overlap in the rainy season. In other words, V. negundo appeared to show niche segregation with R. pseudoacacia and Z. jujuba in the rainy season rather than in the dry season (Table 4). This result is contrary to some studies which found that co-occurring species showed ecological niche segregation of using water in the dry season to avoid drought stress, but showed niche overlap in the rainy season because of taking seasonal rainfall dominance to meet growth demand [48,86,87]. This result could be explained by the dimorphic rooting feature of V. negundo. In the dry season, V. negundo greatly used the tap root system to take up water in deep layer subsoil, which thereby enhanced the water competition with deep root species of R. pseudoacacia and Z. jujuba. In contrast, V. negundo shifted to a lower lateral root to predominantly use surface water, in contrast to the other two species, which still used deeper water in the rainy season.
After afforestation of R. pseudoacacia into the shrub community, the interspecific relationship between R. pseudoacacia and Z. jujuba is among the most important to shape the constructed community due to intensified niche overlap of water source use. In the ridge site, despite both having deeper root systems, Z. jujuba distributed a deeper root and could absorb deeper water than R. pseudoacacia [89][90][91], and thus had lesser niche overlap with the latter. However, at the valley site, the root of Z. jujuba became shallower in relatively improved water and soil conditions [82,89] and thus overlapped with R. pseudoacacia. Consequently, niche of water source use displayed more overlap in the valley than in the ridge. Moreover, shaded by clonally sprouting R. pseudoacacia in the valley, the light-demanding and drought-tolerant Z. jujuba grew better in the ridge than in the valley. On the contrary, R. pseudoacacia is a clonal plant with strong root sprouting [92], which leads to the competitive advantage of R. pseudoacacia over Z. jujuba in the valley. This partly explained why R. pseudoacacia grew better in the valley, whereas Z. jujuba performed well in the ridge (Table 1). Therefore, R. pseudoacacia and Z. jujuba used a similar water source and had higher niche overlap of water use sources. If R. pseudoacacia is used for afforestation, the site should be in the lower slope with benign soil conditions. If planted in the upper slope or ridge, R. pseudoacacia will intensify water competition and the plantation cannot persist in the long term.

Conclusions
The stable isotope-based MixSIAR model was applied to determine the spatiotemporal patterns of the water source use of planted trees R. pseudoacacia with native semi-arbor Z. jujuba and shrub V. negundo. The seasonal variation of soil water availability was similar, but it was significantly lower in the ridge than in the valley. The isotopic composition of the shallow soil water was more enriched. R. pseudoacacia and Z. jujuba can stably derive water from subsoil at a depth of 40-50cm in both dry and rainy seasons, indicating higher water competition between them. However, V. negundo exhibited a flexible water uptake pattern which progressively switched its water source from deep to shallow soil layers as soil water availability improved from dry to rainy season. For this reason, the interspecific relationship of water source use between V. negundo vs. R. pseudoacacia and V. negundo vs. Z. jujuba was different between seasons, which indicates higher water competition in the dry season, but greater water use segregation in the rainy season. The different life form and root distribution of species contributed to various water use strategies, which importantly laid the foundation for the coexistence of the co-occurring species. Although introducing R. pseudoacacia accelerated the establishment of forest plantation and improved soil nutrient conditions, it also intensified competition for water with native species. Our study on seasonal water use source at different microsites provides an insight into the species interaction and stand selection for R. pseudoacacia afforestation in the native shrub community in rocky mountainous areas. R. pseudoacacia is better planted in the valley to avoid intense water competition and control soil erosion.
Author Contributions: Conceptualization, W.Z., P.S. and W.L.; Funding acquisition and project administration, P.S. and J.C.; W.Z. was responsible for the data collection; N.Z. helped with the statistical analysis; W.Z. and S.G. implemented the study site and investigated the vegetation stations; W.Z. wrote the original draft; P.S., J.C. and W.L. conducted the writing-review and editing. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: The data in this study are available on request from the corresponding authors.