Responses to defoliation of Robinia pseudoacacia L. and Sophora japonica L. are soil water condition dependent

Defoliation significantly affected biomass allocation of Robinia pseudoacacia L. and Sophora japonica L., but leaf physiology readjusted to control levels at the end of the experiment. Considering carbon or sink limitation and relative height growth rate, defoliated R. pseudoacacia grew faster than S. japonica under well-watered conditions, while defoliated S. japonica and R. pseudoacacia had similar performance under drought conditions. Climate change may result in increases of both drought intensity and insect survival, thereby affecting both exotic and native trees in warm temperate forests. In this study, we examined the interaction effects of defoliation and drought on an exotic species Robinia pseudoacacia and a native species Sophora japonica in a warm temperate area, to provide a theoretical basis for predicting the distribution and dynamics of the two species under future climate change. In a greenhouse, both species were exposed to three soil moisture (75%, 55%, and 35% of field capacity) and three defoliation treatments (no defoliation, 50% defoliation, and 100% defoliation). Leaf physiology, biomass, and non-structural carbohydrate were determined. Leaf physiology of defoliated trees did not differ from controls trees, but defoliated seedlings allocated relatively more resources to the leaves at the end of the experiment. In well-watered conditions, defoliated R. pseudoacacia was not carbon or sink limited and defoliated S. japonica was carbon limited, while defoliated individuals of the two species were sink limited under drought. Defoliated R. pseudoacacia grow more rapidly than S. japonica in well-watered conditions. Defoliated R. pseudoacacia had a similar growth rate to S. japonica in drought. Defoliation clearly affects biomass allocation of the two species, but not leaf physiology. Considering the carbon or sink limitation, the growth of S. japonica and R. pseudoacacia may be limited by future global climate change scenarios.


Introduction
Climate change models predict that summer in East Asia may become warmer with lower frequency and intensity of precipitation (Dai, 2011;Yu et al. 2014). These changes may cause higher rates of evapotranspiration and lower soil moisture content. The annual average rainfall in northern China has shown a significant downward trend during the period 1951-2010 (Yu et al. 2014). Moreover, the increased annual average temperature is likely to increase insect damage to plants (Netherer and Schopf 2010). Increased temperature enhances the reproductive potential of insects and reduces the time of hibernation (Dale and Frank 2017). As a result, drought and insect defoliation caused by increased temperature are likely to occur simultaneously in northern China under global climate change. Simultaneous drought and insect defoliation affect exotic and native species in the warm temperate forests of northern China (Lombardero et al. 2012). A better understanding of the combined effect of drought and defoliation on exotic and native tree species will be beneficial to understanding the distribution and dynamics of alien species in warm temperate forests under future global climate change.
In trees, photosynthates can be stored as nonstructural carbohydrates (NSC; soluble sugar and starch) that act as a buffer to reconcile the temporal asynchrony between carbon supply and demand (Gricar et al. 2019). NSC fuels respiratory metabolism and supports the buildup of structural biomass. NSC are often used to help distinguish between carbon and sink limitation under stresses or disturbances (Wiley et al. 2017). In general, if growth is carbon limited, NSC will be used for metabolism and maintaining growth, so that eventually NSC will decrease (Barry et al. 2011). For example, when plants are subjected to defoliation NSC will significantly decrease, indicating that the growth of the plant is experiencing carbon limitation (Wiley et al. 2017). If the stored NSC increase or are unchanged but growth is limited under stresses (e.g., drought), plants are considered sink limited (Millard et al. 2007, Piper 2020. In this case, plants are carbon sufficient, and growth would be limited by stress factors. However, long-term drought can result in decreased NSC, since plants need to use the stored carbon to survive and maintain their metabolic processes (Bréda et al. 2006, Galiano et al. 2011). Growth will be carbon-limited in such conditions. Distinguishing carbon and sink limitation is important, because some ecological factors that increase carbon availability will only enhance tree growth under carbonlimited conditions (Wiley et al. 2017).
Defoliation always affects the physiology and growth of plants (Barry et al. 2011;Quentin et al. 2012). After defoliation, plants maintain metabolism by consuming stored NSC (Jacquet et al. 2014), resulting in carbon limitation. Then, plants may increase the net photosynthetic rate of the remaining leaves (Quentin et al. 2012) and decrease the net photosynthetic rate of the bark (Eyles et al. 2009b), thereby increasing carbon absorption to reallocate stored NSC to maintain metabolism and product new organs (Weber et al. 2018). Therefore, after a regrowth period, the growth of defoliated plants may no longer be carbon limited. By this time, defoliated plants with an increased NSC level could eventually increase total biomass, reaching or even surpassing the biomass levels of nondefoliated plants (Eyles et al. 2009a;Barry et al. 2011). Moreover, the degree of growth after defoliation depends on the defoliation intensity, species identity, and environmental context (Chen et al. 2017, Jacquet et al. 2014, Quentin et al. 2012. There are complex interactions between drought and defoliation treatments (Itter et al. 2018, Teskey et al. 2015. The three primary scenarios that follow have been proposed. (i) Defoliation could exacerbate the adverse effects of drought. Defoliated trees may deplete NSC reserves (Maguire and Kobe 2015) and decrease the carbohydrate transport from leaves to roots, thereby limiting root growth capacity and soil water absorption under drought conditions (Teskey et al. 2015). (ii) Defoliation may mitigate the negative impact of drought on plants. Due to the reduction of foliar area, defoliation may reduce transpiration and thereby mitigate the negative impact of drought (Puri et al. 2015, Itter et al. 2018). (iii) Defoliation and drought may have no interactive effects on Populus tremuloides Michx., which may be relevant to the size and age of a tree (Chen et al. 2018). In addition, the interaction of drought and defoliation on a tree also depends on the intensity of drought or defoliation (Jactel et al. 2012). Hence, the effects of drought on the growth of defoliated trees remain debatable.
Robinia pseudoacacia was introduced to China from North America and has become one of the dominant trees in warm temperate forests of China (Wang and Zhou 2000). Since the 1980s in northern China, approximately 8000 ha of R. pseudoacacia has been planted for reforestation due to its capacity for rapid growth (Wang et al. 2020, Zhang andXing 2009). Sophora japonica is a widespread native tree in northern China with significant ecological (Kwon et al. 2016) and medicinal value (Zhang et al. 2014). S. japonica belongs to the same subfamily (Papilionoideae) as R. pseudoacacia. With juicy and palatable leaves, R. pseudoacacia and S. japonica are frequently defoliated by insects (Kwon et al. 2016. Therefore, the exotic R. pseudoacacia and native S. japonica will be exposed to simultaneous insect defoliation damage and drought in the context of climate change. We conducted a greenhouse experiment with three defoliation and three soil moisture treatments to (1) determine whether the growth of R. pseudoacacia and S. japonica is carbon or sink limited under different treatments; (2) reveal the responses of the two species to defoliation; and (3) reveal the responses of defoliated trees to different soil moisture conditions.

Study site and plant materials
This study was conducted at the Fanggan Research Station of Shandong University (36°26′ N, 117°27′ E) in the Central Mountain Region of Shandong Province, China. The site has a warm temperate monsoon climate with an average temperature of 13 ± 1°C and annual precipitation of 700 ± 100 mm (Zhang et al. 2006). The experiment was performed in a greenhouse with a steel pipe frame covered by a plastic film to avoid natural precipitation. The mean air temperature of the greenhouse during the experiment was 29.6°C in the daytime and 20.8°C at night. The greenhouse was kept well ventilated by rolling up the plastic film at the sides during the experiment.
Seeds of S. japonica and R. pseudoacacia were collected in the garden of Dacheng Seed Company (Jinan, China) in autumn of 2014. The seeds were stored at 4°C in a refrigerator. On April 29, 2015, after soaking in distilled water for 24 h, the seeds were transferred to a growth chamber for germination. After 10 days, the seeds started to germinate. Then, we chose healthy and uniformly germinated seedlings to be sown in plastic pots (25 × 21 cm, height × diameter). Each pot was filled with 4.9 kg of air-dried loam and 1.6 kg of sand that remove debris and stones beforehand. The same substrate was used by the growth chamber for germination. Weeds and insects were manually controlled.
Throughout the experiment, the pots received compensatory irrigation by weighing daily at 18:00 to maintain a constant soil moisture level. The weights of pots together with soil of W1, W2, and W3 treatments were supplemented by water to 7.6 kg, 7.3 kg, and 7 kg, respectively. The corresponding soil water contents are: W1 = 16.9%, W2 = 12.4%, W3 = 7.9% according to the formula as follows: W t : the weight of pots together with soil before or after watering during the experiment, W 0 : the weight of pots together with soil before the experiment began.
The weights of pots before watering were recorded at the eighth day, the twenty-second day, the thirty-sixth day, and the fiftieth day after defoliation. The soil water content before watering in different watering treatments was listed in Fig. 1 (using the above formula).
Leaves (excluding apical buds) were removed using scissors from the crown apex downwards in defoliation treatments. Total leaf area removed from individuals was determined in 100% defoliation treatment. For each seedling subjected to 50% defoliation treatment, the shoot was first separated equally into upper and lower parts. The leaves of the upper part were removed and scanned to measure the leaf area. After calculation, the removed leaf area ratio in 50% defoliation treatment was approximately 48% of the total leaf area for each seedling.
Five replicates were included for R. pseudoacacia and S. japonica at each treatment (nine combinations of treatments). As such, there were 90 pots of R. pseudoacacia and S. japonica in total. Each pot had one seedling. Approximately 20cm height of seedlings of both species were subjected to the specified experimental treatments for approximately 8 weeks, from July 12, 2015, to September 6, 2015.

Leaf traits
Gas exchange parameters were measured using a GFS-3000 Portable Photosynthesis System (Walz, Effeltrich, Germany). An LED light source mounted on the GFS-3000 Portable Photosynthesis System was set at a photosynthetic photon flux density of 1200 μmol m −2 s −1 . The CO 2 concentration in the leaf chamber was set as 400 ppm using a CO 2 cylinder, and the temperature in the leaf chamber was set at 28°C to provide a relatively stable environment. Maximum net photosynthetic rate (A max ), stomatal conductance (G s ), and transpiration rate (E) were logged when the system showed stable readings. Gas exchange parameters were measured between 8:30 and 12:00 on sunny days from prior to defoliation, 2, 5, and 8 weeks after defoliation (from July 8-11, July 20-26, August 16-23, and September 1-6, respectively). The three youngest fully developed leaves of three seedlings from each treatment were sampled for measurement.
Chlorophyll a and chlorophyll b concentrations were determined using a 722S visible light spectrophotometer (Leng Guang, Inc., Shanghai, China) after the leaves were extracted in 95% ethanol at end of experiment according to the method of Lichtenthaler and Wellburn (1983). Four replicates of both species in each treatment were sampled, and the concentrations were calculated on a fresh weight basis. The chlorophyll a and b concentrations (chl a+b) and the ratio of chlorophyll a and b (chl a/b) were also calculated.

Growth and biomass partitioning
At the end of the experiment, the seedlings were harvested and divided into root, stem, and leaf. Every part of five replicates for both species in each treatment was weighed after oven-drying at 80°C for 48h. Total biomass (TB) and biomass allocation were calculated as follows: where RB is root biomass, SB is stem biomass, and LB is leaf biomass.

Non-structural carbohydrates
Five replicates for each species in each treatment were used to measure carbohydrate concentration. Leaf, stem, and root of the seedlings were dried at 80°C for 48 h, weighed, and ground to powder using a grinding miller. Following the method of Palacio et al. (2007), 50 mg of dried tissue was prepared for extracting the soluble sugar (SS) in 10 ml 80% (v/v) ethanol at 100°C in a water bath. After centrifugation at 4000 rpm for 10min, the supernatant was separated for SS measurement. The remainder of the sample was digested to glucose with 9.2 mol L −1 perchloric acid to determine the starch (St) concentration. The concentrations of SS and St were determined using a phenol-sulfuric acid colorimetric assay (DuBois et al. 1956) as modified by Buysse and Merckx (1993). NSC concentration was the sum of SS and St.

The relative growth rate in the height
We measured the plant heights at both the start and end of each treatment. The relative growth rate in the height (RGR H , cm day −1 ) was calculated as follows ): Fig. 1 The soil water content (%) before receiving compensatory irrigation in different watering treatments. a the eighth day after defoliation; b the twenty-second day after defoliation; c the thirty-sixth day after defoliation; d the fiftieth day after defoliation. Values are means ± SE (n = 3). W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. W, water treatments; D, defoliation treatments where H 2 and H 1 indicate the seedling height at the end and the beginning of drought × defoliation treatments, and t indicates the time duration (56 days).

Statistical analysis
Three-way analysis of variance (ANOVA) was applied to test the effects of species identity, soil moisture, and defoliation treatments for every performance parameter of R. pseudoacacia and S. japonica. One-way ANOVA and Duncan's multiple range tests were performed to detect significant differences between/among treatments within each species. Repeated measures ANOVA was used to examine the effects of drought, defoliation, and time on gas exchange parameters. Before applying ANOVAs, data were tested for variance normality and homogeneity and were transformed to increase normality and homogeneity of variances if necessary. To test species differences specifically, independent sample t-tests were performed between every parameter of R. pseudoacacia and S. japonica under each water regime × defoliation treatment separately. Dataset used for analysis can be found in the figshare repository (Li et al. 2021). All of the statistical analyses were conducted using IBM SPSS Statistics 22.0 (SPSS Inc., Chicago, IL, USA). Figures were drawn using Origin 9.0 software (OriginLab Co., MA, USA).

Leaf traits
Defoliation had a significant effect on most gas exchange parameters after 2 weeks (Table 1). Defoliation increased A max of R. pseudoacacia in W2 and A max of S. japonica in both W1 and W2 (Fig. 2b, d and e). Moreover, A max of both species began to decline 2 weeks after defoliation, and A max of both species recovered to the control levels 5 weeks after defoliation. Drought had a significant effect on A max at 5 weeks after defoliation (Table 1). Two weeks after defoliation, defoliation increased E of S. japonica and Gs of both species in both W1 and W2 (Appendix Fig. 8 and 9). Five weeks after defoliation, most E and Gs both species recovered to the control levels.
Eight weeks after defoliation, drought continued to have a clear effect on A max , E and G s of the two species (Table 1). Most gas exchange parameters in W1 were clearly higher than those in W3 in any defoliation treatment (Fig. 3a, b, and c). Moreover, A max of R. pseudoacacia was higher than that of S. japonica in all soil moisture and defoliation, except W1D0 (Table 1 and Fig. 3a). Species identity, soil moisture, and defoliation had significant effects on chl a/b, and species identity had a significant effect on chl a+b (Table 1). Chl a/b of R. pseudoacacia in W1D1 and W1D2 were higher than those of S. japonica ( Fig. 3d and Appendix Fig. 10).

Non-structural carbohydrate concentration
Leaf NSC concentrations of both species were unaffected by defoliation or drought. Leaf NSC concentration of R. pseudoacacia was clearly higher than that of S. japonica in all drought and defoliation treatments, except W1D1 and W1D2 (Table 1 and Fig. 4a). Stem NSC concentration was affected by species identity, drought, defoliation, drought × defoliation (the interaction between drought and defoliation), and drought × defoliation × species identity (the interaction among drought, defoliation, and species identity), respectively (Table 1). Under W1, stem NSC concentration of S. japonica remarkably decreased with defoliation intensity (Fig. 4b). However, stem NSC concentration of R. pseudoacacia was unaffected by defoliation. Root NSC concentration was influenced by drought, defoliation, and species identity × drought (Table 1). Root NSC concentration of R. pseudoacacia was unaffected by defoliation and drought (Fig.  4c). Root NSC concentration of S. japonica in W2 and W3 was clearly higher than that in W1 (Fig. 4c).

Growth and biomass allocation
The species identity, soil moisture, defoliation, and species identity × soil moisture clearly affected leaf biomass of both species (Table 1). In W1 and W3, LB of R. pseudoacacia showed no significant differences among the three defoliation treatments; in W2, LB of R. pseudoacacia in D0 was clearly higher than that in D2 (Fig. 5a). LB of S. japonica in D0 was higher than that in D2 under W1 and W2 treatments (Fig. 5a). And LB of R. pseudoacacia was higher than that of S. japonica in W1 and W3, except W1D0.
TB of the two species was significantly affected by species identity, soil moisture, and defoliation (Table 1). In W1, TB of R. pseudoacacia showed no significant differences among the three defoliation treatments; in W2 and W3, TB of R. pseudoacacia in D0 was clearly higher than that in D2 (Fig. 5b). TB of S. japonica in W1 and W2 decreased with defoliation intensity (Fig. 5b). TB of R. pseudoacacia was higher than that of S. japonica in W1 and W3 (Table 1 and Fig. 5b).
The soil moisture and species identity clearly affected most of the biomass allocation of both species (Table 1). Drought clearly restricted LMR of the two species (Fig.  6a). RMR and R/S of S. japonica were lower in W1 than that in W2 and W3, while those of R. pseudoacacia was unaffected by drought (Fig. 6 c and d).

The relative growth rate in the height
RGR H was affected by drought, defoliation, and species identity (Table 1). In W1, RGR H of S. japonica was lower than that of R. pseudoacacia. In W2 and W3, RGR H of R. pseudoacacia has no significant difference from that of S. japonica (Fig. 7).

Carbon or sink limitation after defoliation and drought
Our study showed that leaf, stem, and root NSC concentrations of R. pseudoacacia were not affected by defoliation or drought (Fig. 4), indicating that defoliated R. pseudoacacia had recovered the carbon concentration of each organ to the control level, consistent with studies on Nothofagus pumilio and Pinus sylvestris (Piper et al. 2015, Susiluoto et al. 2010. Moreover, TB of defoliated R. pseudoacacia was not significantly different from that of undefoliated trees under high soil moisture conditions, suggesting that the growth of defoliated R. pseudoacacia in high moisture conditions is not carbon or sink limited. However, stem NSC concentration of defoliated S. japonica under high soil moisture conditions was remarkably lower than in control (Fig. 4b). Combined with the lower TB of defoliated S. japonica than that of undefoliated one, we inferred that defoliated S. japonica did not have enough carbohydrate stored to support growth, which is consistent with the definition of carbon limitation (Wiley et al. 2017). Hence, Leaf, stem, and root NSC concentrations of two species in drought had no significant difference with those of seedlings in high soil moisture conditions, indicating that seedlings with or without defoliation in drought had enough stored carbohydrate to support growth (Piper 2020). Moreover, TB of both defoliated species in drought were lower than those of undefoliated trees in the high moisture treatments, meaning that growth of defoliated trees was inhibited by drought and defoliation. In drought, defoliated R. pseudoacacia and S. japonica had enough stored carbohydrate, but could not recover their full growth, suggesting that drought inhibited the process of plant synthesis of biomass using available carbohydrates. That was in line with the characteristics of sink limitation (Wiley et al. 2017).

R. pseudoacacia and S. japonica responses after defoliation
In our study, seedlings of defoliated R. pseudoacacia and S. japonica had a transient upregulation of A max at 2 weeks after defoliation (Table 1 and Fig. 2). The results are consistent with studies in Eucalyptus globulus and Pinus sylvestris (Quentin et al. 2010, Schönbeck et al. 2020). The reason for stronger photosynthesis could be that defoliated seedlings had insufficient stores of NSC in organs to maintain growth (Wang et al. 2020). Pinkard et al. (2011a) showed that lower NSC concentration in plant organs may stimulate A max producing more carbohydrates. Therefore, seedlings increased Gs (Table 1 and Appendix Fig. 9), which was strongly correlated with the increase in A max of defoliated Eucalyptus globulus (Quentin et al. 2012) to increase the capacity of photosynthesis at leaf level for increased uptake of carbohydrate.
A max and E of the two species were not significantly different between defoliation treatments at 8 weeks after defoliation (Table 1, Fig. 2, and Appendix Fig. 8), indicating that these parameters of both species had grown to the control levels. Our results are consistent with the study on Eucalyptus globulus (Pinkard et al. 2011b). In our study, this result may be related to unchanged leaf NSC concentration under defoliation treatments (Fig.  4a). According to the positive relationship of stored NSC and A max (Pinkard et al. 2011a), we inferred that unchanged leaf NSC concentration may have no positive feedback on A max (Iqbal et al. 2012). In addition, we observed that chl a+b was unchanged by defoliation (Table  1 and Appendix Fig. 10). Pinkard et al. (2011b) showed that A max and chlorophyll content usually had a strong Fig. 2 A max of R. pseudoacacia and S. japonica grown in different water and defoliation treatments over the experimental period (mean ± SE, n = 3). Different lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05 Li et al. Annals of Forest Science (2022) 79:18 positive relationship. Hence, the unchanged chlorophyll content may also have contributed to the unchanged A max at 8 weeks after defoliation. Leaf NSC concentrations of the two species were unchanged by defoliation, while stem and root NSC concentrations of the two species decreased with defoliation intensity, showing NSC concentrations of the two defoliated species did not completely recover. Importantly, this result also suggested that defoliated seedlings invested more NSC to leaf than to stem and root. Some studies have found that producing new organs requires the consumption of carbohydrates from the nearest storage pool (Dietze et al. 2014). Thus, higher leaf NSC concentration may mean more new leaves of seedlings that can improve the carbon absorption of the whole seedling and fix more carbon in order to grow.
With regard to biomass allocation, we found that both species were unaffected by defoliation. Indeed, when calculating LB in the defoliation treatments, we neglected the removed leaf biomass at the beginning of defoliation. Therefore, LB of defoliation treatments may be lower than that of control. After a regrowth period, LMR of defoliated trees was unchanged, indicating that they may allocate more resources to the leaves. Our finding is consistent with another study showing that defoliation may increase the leaf or aboveground biomass ratio to increase carbon assimilation of the whole tree (Eyles et al. 2009a).

Defoliated R. pseudoacacia and S. japonica responses under different soil moisture conditions
There were no significant interaction effects between defoliation and drought treatments on most of the parameters (except for E and G s at 5 weeks after defoliation and stem NSC concentration at 8 weeks after defoliation) in this study. However, some studies have identified antagonistic or synergistic effects of biotic and abiotic stresses on growth in different woody plants, such as in Fig. 3 A max , E, G s , and chl a/b (panels a-d) of R. pseudoacacia and S. japonica grown in different water and defoliation treatments (mean ± SE, n(A max , E, and G s ) = 3, n(Chl a/b) = 5) at 8 weeks after defoliation. Different uppercase and lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test, respectively. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05 Li et al. Annals of Forest Science (2022) 79:18 Eucalytus globulus and Quercus rubra (Kolb et al. 1999, Quentin et al. 2012. In our study, R. pseudoacacia and S. japonica belong to the Leguminosae, which may result in both species had similar interaction effects between defoliation and drought conditions. Therefore, we inferred that our results could be attributed to species' identities or family specific. Stem NSC concentration of S. japonica decreased with defoliation under well-watered conditions, while that of R. pseudoacacia was unchanged by defoliation. This may have resulted from the higher A max of R. pseudoacacia compared to S. japonica (Fig. 3a) causing relatively more leaf NSC of R. pseudoacacia to be allocated to stems. However, stem NSC concentrations of R. pseudoacacia and S. japonica were unaffected by defoliation under drought suggesting that NSC is crucial to maintaining hydraulic functioning in trees under water stress (Jacquet et al. 2014).
Defoliated R. pseudoacacia and S. japonica had different RGR H under different water conditions: RGR H of defoliated R. pseudoacacia was higher than that of S. japonica in high moisture conditions, while R. pseudoacacia and S. japonica did not differ significantly in RGR H in drought. The results suggested that the exotic R. pseudoacacia was able to take advantage of increased water resources and enhance its growth after defoliation. Caramaschi et al. (2016) also reported that exotic species after defoliation had a greater ability to grow than native species. This is one of the characteristics of exotic species, i.e., that they often exhibit a greater response than native species when resource availability increases (Leishman et al. 2010, Luo et al. 2015. However, compared with native species, most exotic species cannot capitalize and maintain greater performance in a stressful environment (Dawson et al. 2012), which is supported in our study. In addition, A max and chl a/b of the Fig. 4 Comparisons of NSC concentration in leaf (a), stem (b), and root (c) of R. pseudoacacia and S. japonica from different water and defoliation treatments (mean ± SE, n = 5) at the end of the experiment. Different uppercase and lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test, respectively. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05  5 Comparisons of LB (a) and TB (b) of R. pseudoacacia and S. japonica grown in different water and defoliation treatments (mean ± SE, n = 5) at 8 weeks after defoliation. Different uppercase and lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test, respectively. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05. Fig. 6 Comparisons of biomass allocation of R. pseudoacacia and S. japonica grown in different water regimes (mean±SE, n = 5) at the end of experiment. a LMR, b SMR, c RMR, and d R/S. Different letters denote significant differences in each species at p ≤ 0.05 by Duncan's test. W1, W2, and W3: 75%, 55%, and 35% of field capacity exotic R. pseudoacacia was higher than that of the native S. japonica at 8 weeks after defoliation under wellwatered conditions (Table 1 and Fig. 3 a and d), which may be a key reason for stronger RGR H of R. pseudoacacia.
In addition, we found that leaf NSC concentration of R. pseudoacacia had no significant difference with that of S. japonica in high soil moisture conditions. Combined with the high A max of R. pseudoacacia in wellwatered conditions, these results indicated that R. pseudoacacia used more carbon to supply the growth of new leaves (Dietze et al. 2014). Meanwhile, more new leaf production of R. pseudoacacia can improve carbon fixation of whole seedlings, and this may be another reason for the stronger RGR H of R. pseudoacacia under wellwatered conditions.
Interestingly, RGR H of S. japonica and R. pseudoacacia did not have a significant difference in drought (Fig. 7). S. japonica in drought conditions had the higher R/S, RMR, and root NSC concentration than that in the wellwatered conditions, suggesting that S. japonica allocates more resources to roots in drought conditions. This could improve root osmotic pressure (higher root NSC concentration means higher solute concentration), taking up more water and alleviating the negative effects of drought (Guo et al. 2019, Jacquet et al. 2014. Comparing the unchanged R/S, RMR, and root NSC concentration of R. pseudoacacia in drought, we inferred that R. pseudoacacia had weaker drought tolerance than S. japonica. Although R. pseudoacacia had higher carbon absorption capacity at leaf level (A max ) and leaf biomass than those of S. japonica under drought conditions (Fig.  3a), weak drought tolerance of R. pseudoacacia limited its ability to grow after defoliation. This may be the main reason why there is no significant difference in the RGR H of the two species under drought treatments.

Conclusion
Our study showed that A max was upregulated at 2 weeks after defoliation, indicating that both species had a transient increase in carbon fixation at the leaf level. At the end of the experiment, the leaf physiology of both defoliated species had recovered to the levels of undefoliated trees. Defoliated seedlings allocated more resources to the leaves to increase carbohydrate production. In wellwatered conditions, R. pseudoacacia achieved clearly greater RGR H than S. japonica, while R. pseudoacacia had similar RGR H to S. japonica in drought. In wellwatered conditions, defoliated R. pseudoacacia did not show carbon or sink limitation. However, both species were sink limited in drought. In conclusion, the growth of both R. pseudoacacia and S. japonica may be limited under future drought and defoliation. Nevertheless, to fully reveal the different responses between sympatric R. pseudoacacia and S. japonica, long-term experiments and monitoring data are needed. Fig. 7 RGR H of R. pseudoacacia and S. japonica grown in different water and defoliation treatments (mean±SE, n = 5). Different uppercase and lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test, respectively. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between treatment within each species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05 Different lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05  10 Chl a+b of R. pseudoacacia and S. japonica grown in different water and defoliation treatments at the end of the experiment (mean ± SE, n = 5). W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation Fig. 9 G s of R. pseudoacacia and S. japonica grown in different water and defoliation treatments over the experimental period (mean ± SE, n = 3). Different lowercase letters denote significant differences in each species at p ≤ 0.05 by Duncan's test. W1, W2, and W3: 75%, 55%, and 35% of field capacity. D0, D1, and D2: control, 50% defoliation, and 100% defoliation. Asterisks indicate significant differences between species in each drought and defoliation treatment separately with independent sample t-test: **p ≤ 0.01, *p ≤ 0.05