Hydraulic Constraints to Whole-Tree Water Use and Respiration in Young Cryptomeria Trees under Competition

Although extensive studies have focused on carbon and water balance from aboveground measurements, the link between the belowground and aboveground processes deserves greater attention. In this context, the aim of this work was to assess the bi-directional feedback between whole-plant respiration and transpiration. The study was performed on 25 saplings of Sugi (Japanese cedar, Cryptomeria japonica D. Don), including dominant and suppressed individuals (total fresh weight ranging between 0.2 and 8.0 kg). During one week, the integrated water use (WU) was determined using the Deuterium dilution method. After this, the trees were uprooted and the root, stem, and leaf respiration were measured using incubation chambers and CO2 infrared sensors. The stem and root respiration followed a power response to mass (power exponent b < 1), implying a decline in mass-specific respiration with size. Conversely, the leaf respiration followed a near-linear increase with size (power exponent b ≈ 1), but was negatively affected by the stem density, indicating the hydraulic limitations of the leaf metabolism. The water use followed a power response with the tree size (b < 1), showing a decline in the transpiration per leaf mass with the tree size, but was also negatively correlated with the stem density. Our results indicate that dominant trees are more efficient in the use of water, and highlight the role of hydraulic limitations to leaf metabolism in suppressed trees.


Introduction
One of the main challenges in ecology is to predict the response of forests to global change, and in particular, the interaction between CO 2 fertilization, increased temperatures, and drought [1][2][3][4][5][6][7][8].Although extensive studies have been devoted to carbon and water balance from aboveground measurements at different scales, the interaction between the belowground and aboveground processes demand greater attention [9,10], in particular, the trade-offs and synergies between water use (WU) and whole-tree carbon metabolism and redistribution [11][12][13][14][15].For example, the root and soil respiration depends on the export of assimilates from the leaves to the roots [11,12,15], but is also modulated by the local soil water content [13,14].Furthermore, when water is limited, plants tend to increase the carbon allocation to the roots [16,17], but the overall carbon fixation may be limited because of the reduced leaf area and drought-induced photosynthesis limitations [11,12,18,19].Optimality theory postulates that plants should allocate more biomass to the shoot when carbon fixation (i.e., energy) is the limiting factor, as well as to the roots under water-or nutrient-limited conditions [20][21][22][23].The trade-off between the belowground demand for photosyntates and the aboveground growth is particularly critical in young competing trees, where water and light limitations often impose opposing constraints to tree development [24][25][26][27][28].In this context, characterizing the dominance-dependent association between the carbon allocation patterns and whole-tree water uptake capacity would improve our mechanistic understanding of forest regeneration dynamics in response to global change [25,29].
So far, the changes in carbon allocation have been mainly assessed on a dry-mass basis, and only recently has the relevance of metabolic partitioning among different plant organs been highlighted [30,31].The metabolic rate (i.e., respiration rate, R) of individual organisms' scales with body size (M), can be described as the simple power function of body mass, as follows: where a is the normalization constant (antilog of the intercept in a log-log plot), and b is the scaling exponent (slope in a log-log plot) [30, [32][33][34][35].The scaling exponent is often less than unity, thus indicating that the mass-specific metabolic rates decrease with size, and the theoretical models have predicted a universal scaling exponent of 3/4 [32][33][34].However, empirical evidence in plants suggest that the allometric scaling of the metabolic rates may differ among plant organs, because of their particular functional requirements, as well as among individuals with different resource limitations [9,30,36].
In this context, the aim of this study is to explore the size scaling of respiration in different plant organs, as well as the whole-tree water use for young, competing individuals of Sugi (Japanese cedar, Cryptomeria japonica D. Don).Sugi is a long-living, fast-growing species adapted to temperate, maritime climates, but suffering seasonally high water demanding conditions [27,37,38].As a tall tree (>50 m), Sugi has developed a large stem and leaf storage capacity in order to overcome the hydraulic limitations to tree growth, maintaining high transpiration rates along the canopy [39][40][41][42].Considering the trade-off between light interception and water securing as a major constrain for the distribution of resources to the aboveground and belowground tissues, we hypothesize a progressive increase in the root carbon allocation from small, slow-growing, suppressed individuals (i.e., more light-limited) to large, fast-growing, dominant trees (i.e., more water-demanding).Bearing in mind the aforementioned adaptations to maintain hydraulic function in tall trees, we also postulate that the decreasing shoot/root metabolic ratio with tree size should translate into a near-linear (b ≈ 1) scaling of the whole-tree water use.Our first hypothesis will be assessed not only with regard to the biomass ratios, but also in terms of metabolic (respiration) rates, as a better proxy for the energy partition among tissues.To test our second hypothesis, we will assess the mass-and metabolic-dependent scaling of whole-tree water use, as estimated for the same individuals, using a Deuterium dilution technique.

Plant Material and Study Site
The study was performed in the nursery field of the Faculty of Agriculture, Yamagata University (Tsuruoka campus, 38 • 43.98 N, 139 • 49.39 E, 17 m.a.s.l.).The monthly mean temperature in Tsuruoka ranges from 0.3 • C in January to 25.2 • C in August (Annual mean: 12.3 • C), and annual precipitation amounts to 1948 mm.For the experiment, we used 7 to 10 year old saplings of Cryptomeria japonica D. Don, growing in a 2 × 8 m rectangular plot under high competition (ca.60,000 trees ha −1 ), selecting 25 individuals representative of the range of tree size and dominance status observed in the plot (tree height ranging from 0.75 to 2.81 m).Selecting 25 individuals representative of the range of tree size and dominance status observed in the plot (tree height ranging from 0.75 to 2.81 m; full dataset is available as Supplementary Materials, Table S1).

Quantification of Water Transport and Storage through Deuterium-Dilution Method
A deuterium labelling experiment was carried out between the 16 and 22 September 2014, as described in the literature [43].Briefly, a known amount of Deuterium oxide (D 2 O 99.9 atom% D; Sigma-Aldrich, Saint Louis, MO, USA) was injected into the base of the tree trunk, followed by the collection of transpired water (condensed inside plastic bags, located at 2/3 of tree height) during the subsequent days.
The weather was sunny to moderately cloudy, with scattered light rains from 18 to 19 September (the environmental conditions and sampling schedule are summarized in Figure 1).The tracer was injected on 16 September, between 10:12 (tree 1) and 12:20 (tree 25), starting immediately with the first water collection.Following the literature [44], the amount injected was proportional to the basal area of the tree (112-1605 mm 2 ), ranging between 0.89 and 9.99 g D 2 O.During the first three days, the water was collected twice (morning and afternoon), plus two additional afternoon water collections on 19 and 22 September.Each water collection lasted between two and four hours, registering the start and end of the sampling period individually for each tree.The light environment in the canopy of each individual tree was characterized through the measurements of the Photosynthetic photon flux density (PPFD, µmol m −2 s −1 ) at 2/3 of the tree height (LI-190 quantum sensor, LI-COR Inc., Lincoln, NB, USA).The underlying principle of the Deuterium-labelling method is that the highly transpiring trees will show a faster labeling of the transpired water, but with less amplitude and for a shorter period of time, due to the dilution of the D 2 O injected over a greater pool of the transpired water (Figure 2).The integration over time of the amplitude of the D 2 O label, relative to the tracer amount, was used to estimate time-integrated tree water use (WU, g H 2 O day −1 ) [45,46], as follows: where M is the mass of the tracer injected (in g), C i is the D 2 O concentration at time i (g D 2 O g −1 H 2 O), and ∆t i is the duration of the time interval i (days) between water collection times, as previously recorded for each tree.The time response over time of the D 2 O tracer was used to determine three water transport and storage attributes, velocity, half-life, and residence time [47].The tracer velocity (Vel), a surrogate for the linear sap velocity, was calculated by dividing the distance between the injection and collection point (2/3 of tree height), by the time required to reach 10% of the maximum tracer concentration, determined from a linear regression between the time and standardized values during the upward section of the curve.The tracer half-life (t 1 2 ) and tracer residence time (RT) were estimated from an exponential decay function fitted to the downward section of the curve, as follows: where C 0 and C stand for tracer concentrations at time zero (i.e., peak maximum) and time t, respectively, and λ is the decay constant.From this equation, t 1 2 and RT were solved as the time needed to reach half of the maximum and pre-injection background values (D/H atom% = 0.015), respectively.

Isotope Analyses
The isotopic composition of the water samples was analyzed using the H 2 /H 2 O equilibration method in a Gas Bench (Thermo Fisher Scientific, Waltham, MA, USA) attached to a Delta V isotope ratio mass spectrometer (Thermo Fisher Scientific, Schwerte, Germany), at the Faculty of Environmental and Earth Science, Hokkaido University.The isotopic composition was expressed relative to an international standard (VSMOW, Vienna Standard Mean Ocean Water) in per mil notation (‰), as follows: where R sample and R vsmow are the isotope ratio of water (D/H) for the sample and for the Vienna Standard Mean Ocean Water (VSMOW), respectively.Analytical error, determined as the standard deviation of the internal standards, was less than 2‰.To avoid a peak saturation because of the high D/H ratios in labelled water, the isotope ratios were determined on 1:10 dilutions of the original samples, using water with a known isotopic composition.

Measurement of Respiration Rates and Biometric Variables
At the end of the labelling experiment, the trees were uprooted and the respiration rates were determined separately for the stem, roots, and leaves, following the literature [30].The measurements were taken during three consecutive days (24-26 September), and within 20 min of the tree uprooting, keeping the plant material wrapped in wet paper to prevent transpiration.The different tissues were enclosed in custom-made incubation chambers, including forced air circulation with an electric fan, with a volume adapted to the size of each tree part (2.6 to 120 L).The increment rate of the CO 2 concentrations within the chamber was then measured every five seconds using an infrared CO 2 analyzer (GMP343, Vaisala, Helsinki, Finland), and subsequently normalized to 20 • C, assuming Q 10 = 2 [30].For each individual, the fresh weight of each part, tree height, and length of the photosynthetic section of the stem were measured.The stem basal area (BA, mm 2 ) and basal area increments (BAI, mm 2 year −1 ) were calculated as the area of an ellipse, with the maximum and minimum diameters within each tree ring as the two axis.As a result of the small size of the trees, we determined the BA and BAI at 10% of tree height, instead of the standard breast height.A subsample of leaves, coarse roots, and a section of the base and the apex of the stem were weighted fresh, and then dried for 48 h at 60 • C to determine water concentration (WC, % fresh weight).For the two stem sections, the stem density (kg m −3 ) was estimated as the quotient between the dry weight and the volume of a truncated cone (including bark).The total dry weight for the root, stem, and leaves was estimated from the fresh weight measurements and WC.Finally, we calculated the mass-specific respiration and mass-specific water use on a fresh-weight and dry-weight basis.

Data Analysis
Following the literature [30], we assessed the mass-scaling relationships for respiration and water use using a simple-power function on the log-log coordinates of the log-transformed version of Equation (1) (N = 25).We used an analysis of variance (ANOVA) to assess the effect of the tree ages on the size-related variables, and the principal component analysis (PCA) and linear regression to assess the association between the variables.The means are shown together with standard errors of the mean.The normality of variables was tested with the Shapiro-Wilks test.The size-related variables, respiration rates, WU, and Vel were log-transformed to satisfy the normality assumptions in the statistical analyses.The rest of variables were used without transformation.

Tree Morphometric Variables
The summary statistics of the main morphometric variables are shown in Table 1.The sampled trees showed a 3.7-fold range in tree height (0.75-2.81 m), and a 37-fold range in the whole-plant fresh weight (0.21-7.98 kg).The stem size showed the largest relative differences (52-and 58-fold, in weight and volume, respectively), compared with the root and leaf weight (43-and 36-fold variation, respectively).Comparatively, the range of variation in the recent BAI was much larger than the range of variation in the total BA (34-and 16-fold, respectively, see Table 1).The proportion of BA formed during the last three years ranged from 30 to 75% of the total BA, and explained 45% of the variability in the total shoot biomass (y = 3.78x 2.47 , R 2 = 0.45, p < 0.001).Conversely, we did not find significant differences among the age classes in any of the size-related variables studied (p = 0.214-0.588),indicating that the size differences were largely determined by the growth rates.The distribution of the fresh weight among the different plant organs varied with the tree size (Figure 3a).The changes in the mass allocation with the tree size did not follow a linear trend; we found the largest differences among the smaller trees (ca.<2 kg), and was best depicted on the log-log coordinates (Figure 3b).We found a significant increase with the total fresh weight in the stem/root ratio (y = 1.016x 0.206 , R 2 = 0.35, p = 0.001), and a significant decrease in the leaf/stem ratio (y = 2.129x −0.222 , R 2 = 0.47, p < 0.001), but no significant changes in the leaf/root ratio (p = 0.499).Conversely, the shoot/root ratio did not change significantly with the tree size (3.27 ± 0.126; p = 0.204).These results indicate that the leaf mass was adjusted to the root size, coupling the water demand with the water uptake, whereas the conducting tissue (stem) increased with the size, at the expense of the water uptake capacity (root size).

Respiration Rates and Their Association with Tree Size
Overall, the respiration rates showed a smaller relative range of variation, as compared to the range in the fresh weight.The sampled trees showed a 25-fold increase in the total respiration, ranging between 0.10 and 2.49 µmol CO 2 s −1 .The lack of a linear increase in the respiration was in accordance with the theoretical size-scaling of the respiration rates, which could be described with a power function of the fresh weight (Figure 4).The shoots showed a steeper and more linear (i.e., b coefficient, closer to 1) increase in the respiration with weight (y = 0.361x 0.904 , R 2 = 0.98, p < 0.001), as compared to the roots (y = 0.224x 0.842 , R 2 = 0.91, p < 0.001) (Figure 4a).However, this was mainly determined by the steep response of the leaf respiration (y = 0.375x 0.932 , R 2 = 0.94, p < 0.001), whereas the response of the stem was similar to the root (y = 0.260x 0.847 , R 2 = 0.95, p < 0.001) (Figure 4b).Indeed, the 95% confidence interval for b in the two woody tissues was significantly lower than the unity (0.73-0.96 in the root; 0.76-0.93 in the stem), whereas this was not the case for the leaves (0.83-1.03).This indicates that the mass-specific metabolic rates in the woody tissues decline significantly with organ size, but remain nearly constant in the photosynthetic tissues.The existence of a near-linear increase with the mass in the total leaf respiration, but not in the total wood respiration, was further evidenced when comparing the mass-specific respiration rates with the organ size (Figure 5).On the log-log scale, both the shoot and root specific respiration decreased significantly with the fresh weight (y = 0.361x −0.096 , R 2 = 0.32, p = 0.025 in the shoot; y = 0.224x −0.158 , R 2 = 0.27, p = 0.064 in the root) (Figure 5a).However, this trend was significant for the stem (y = 0.260x −0.153 , R 2 = 0.37, p < 0.001), but not for the leaves (p = 0.125) (Figure 5b).Due to the differential response of respiration to the mass in the green and woody tissues, the respiratory partitioning among the plant organs followed less consistent trends than those observed for the fresh weight.We only found a weak association between the total fresh weight and the stem/root respiration ratio (y = 1.165x 0.172 , R 2 = 0.20, p = 0.020).Neither the leaf/stem respiration ratio nor the leaf/root respiration ratio changed significantly with the total fresh weight (p = 0.162 and p = 0.196, respectively).Considering the association between the dry weight and respiration, as well as the dry weight allocation among the organs, we observed nearly the same trends as for the fresh weight, although slightly weaker (data not shown).

Water Use and Water Storage
On average, the label peaked within one day after injection (Table 2; for peak plot see Figure A1 in Appendix A).The maximum D/H in the transpired water was 3-to 41-fold larger than the background (D/H atom % = 0.015), approaching pre-injection values in the last sampling time, about six days after injection (Table 2).
As it occurred with the respiration, the Vel and WU showed a smaller relative range of variation than the tree size (32-and 18-fold increase, respectively; see Table 2).The Vel and WU are surrogates for the linear sap velocity and sap flow, respectively, and were consistently associated through a power function (WU = 214.9vel 0.656 , R 2 = 0.46, p < 0.001).The leaf-specific WU (i.e., transpiration rate per unit of leaf mass) showed a much smaller range of variation (2.6-fold).The two variables associated with the water storage dynamics, t 1 2 and RT, also showed a relatively small range of variation (6-and 4-fold, respectively).Both the WU and Vel followed a power response to the tree size (Figure 6).For WU, the strongest association was found with the leaf fresh weight (y = 1162x 0.861 , R 2 = 0.95, p < 0.001), followed by the stem and root fresh weight (y = 1737x 0.687 , R 2 = 0.94, p < 0.001 and y = 2136x 0.822 , R 2 = 0.91, p < 0.001, respectively; see Figure 5a).Similar responses were found considering either the shoot or total fresh weight (y = 815x 0.807 , R 2 = 0.96, p < 0.001 and y = 653x 0.816 , R 2 = 0.95, p < 0.001, respectively).Although the b coefficient was significantly larger in the leaves and roots than in the stem, in all cases, the 95% confidence interval for b was significantly lower than the unity (0.78-0.95 in leaves; 0.61-0.76 in stem; 0.71-0.93 in roots).In other words, the increase in the transpiration with size was not linear, being less pronounced in the larger individuals.Accordingly, the leaf-specific WU showed a negative association with the leaf fresh weight (y = 1162x −0.139 , R 2 = 0.33, p = 0.003) (i.e., transpiration rates per unit of leaf mass were smaller in the largest tree).The Vel showed the strongest association with the root fresh weight (y = 13.67x 0.616 , R 2 = 0.48, p < 0.001), followed by the leaf and stem fresh weight (y = 8.32x 0.580 , R 2 = 0.40, p < 0.001 and y = 10.76x 0.452 , R 2 = 0.38, p < 0.001, respectively; Figure 5b).The WC and stem density showed relatively small variability across individuals, with a 1.2-and 1.7-fold range, respectively (Figure 7).However, we observed significant differences in the WC among the plant tissues (p < 0.001; leaf > root > stem apex > stem base), as well as a significant decrease in the stem density from the base to the apex (p = 0.014).The WC in the roots and stem tissues increased with tree fresh weight (y = 0.577 + 0.04698 × log(x), R 2 = 0.38, p = 0.001, in the root, y = 0.505 + 0.05009 × log(x), R 2 = 0.51, p < 0.001 in the stem base, and y = 0.555 + 0.06326 × log(x), R 2 = 0.72, p < 0.001, in the stem apex; see Figure 7a).Conversely, the leaf WC did not change significantly with the tree size (p = 0.995).The stem base density decreased significantly with the tree fresh weight (y = 642 − 102.4 × log(x), R 2 = 0.34, p = 0.002; see Figure 7b), but this trend was not significant for the stem top density (p = 0.242).
Due to the differential response of the tissue the WC, changes in the relative contribution to the water storage per tissue (i.e., total water per plant organ) were slightly stronger than in the fresh and dry weight allocation.We found a significant increase in the stem/root water storage ratio with total fresh weight (y = 0.933x 0.217 , R 2 = 0.37, p = 0.001), and a significant decrease in the leaf/stem ratio (y = 2.627x −0.270 , R 2 = 0.53, p < 0.001), but no significant changes in the leaf/root ratio (p = 0.158).Unlike WU and Vel, the two variables associated with water storage dynamics (t 1 2 and RT) were not associated with the tree size.Among the different morphometric variables tested, we only found a weak positive association between the stem base density and both t 1 2 (log(y) = −0.809+ 0.001099x, R 2 = 0.16, p = 0.048) and RT (y = −1.05+ 0.00607x, R 2 = 0.18, p = 0.035).Conversely, the stem base density was negatively associated with WU and Vel (log(y) = 4.684 − 0.00289x, R 2 = 0.37, p = 0.001 and log(y) = 2.56 − 0.00280x, R 2 = 0.33, p = 0.003, respectively).This suggests that the larger water retention times and limited water transport in the small trees was partly due to the reduced water conductivity in the dense stems.Associations between the dry weight and hydraulic traits followed very similar patterns to those found for the fresh weight, but slightly weaker (data not shown).

Association between Hydraulic, Metabolic, and Morphometric Traits
A principal component analysis, including the most relevant measured traits, revealed consistent associations between the respiration rate and the hydraulic traits (WU, WC, RT, and Vel, see Figure 8).The first axis, explaining 51% of the total variability, was mainly driven by size, and showed a tight association between the total fresh weight and total respiration, as well as the WU and stem WC, but negatively correlated with the leaf to stem ratio, stem density, and specific respiration in woody tissues (roots and stem).Along this axis, we also found negative associations between the tree size and leaf-specific WU and RT, and weak positive associations with the stem to root ratio and Vel.The second axis, explaining 13% of the total variability, was dominated by the tight association between the leaf-specific respiration and leaf-specific WU, and opposed to the stem density and RT.Along this axis, we also found positive associations between the leaf-specific traits and both the tracer velocity and stem specific respiration, and negative assosciations with stem to root ratio.Looking in detail at the association between hydraulic traits and respiration, we found that respiration rates followed a power response to WU, being stronger and steeper (b > 1) in stem than in the leaves and roots (b ≈ 1; Table 3).Vel was best correlated with root and leaves respiration rates, but the relationship was steeper for stem and leaves than for the roots (Table 3).Respiration was positively associated with tissue WC in the stem and roots (R2 = 0.61 and 0.36, respectively), but not in the leaves (Table 3).Conversely, the respiration in all of the tissues was negatively correlated with the stem base density, but the association was stronger in the leaves (R 2 = 0.41) than in the stem and roots (R 2 = 0.32 and 0.26, respectively; see Table 3).

Biomass Allocation Patterns Differ in Dominant and Suppressed Trees
In the highly competitive conditions of our study, the lack of age effects implies that the observed trends in the size and metabolic activity were driven by a dominance gradient (i.e., from small, slow-growing, suppressed trees, to large, fast-growing, dominant trees), rather than by developmental changes.As a consequence, the differences in the stem BA were largely determined by the differences in the recent growth (30-75% of the total BA was formed during the last three years, and this ratio was independent from the tree age).The large competition is further confirmed by the comparatively lower size of the measured individuals than in the other studies at the same age [42,48], and the lower growth rates than predicted by allometric equations (e.g., [37,49]).
Over the studied range, we did not find significant changes in the shoot/root ratio with the tree size, suggesting that this ratio is not affected by the dominance status in the juvenile stands.The shoot to root ratios in our study (2.3-4.5) were in the range of those reported in the literature for young Sugi trees (2.9-3.6 in <1-year seedlings [27]; 3.1-3.7 in <40 age classes [49]), but smaller than in the older trees (9.2 in <50 age class [49]).Contrasting with the invariance of the shoot/root ratio, there was a substantial change in the allocation pattern among the woody tissues (stem/root ratio), resulting in larger stem growth rates in the dominant individuals, but maintaining a constant root to leaf ratio.This implies a tight regulation of the leaf mass to adapt to the root size, by means of self-pruning, and translates into a lower leaf/stem ratio in the larger, dominant individuals than in the smaller, suppressed individuals.
Overall, our results contradict our initial hypothesis, indicating that the dominant trees do not allocate more carbon to the roots than to the suppressed trees.In more energy-limited than water-limited environments, as in our study, the root size does not limit the stem growth of the dominant trees.In this context, light interception and height growth becomes the main priority, at the expense of the root development [22,50,51].However, the root size may still constrain the canopy development (i.e., leaf mass), as a strategy to couple the water losses with the root uptake capacity, thus avoiding hydraulic limitations in tall trees.These trends may not be restricted to young trees under competition; Lim et al. [49] also reported a constant leaf/root ratio from <20 to <50 age classes, together with declining leaf/stem ratios and increasing stem/root ratios.Our results also evidence the risk of using the shoot to root ratio as the only indicator of carbon allocation patterns, obviating the contrasting physiological function of the stem and leaves.

Metabolic Rates Tend to Compensate for Changes in Biomass Allocation
The size-dependent changes in the partition of the metabolic rates (i.e., respiration) differed from those found for the biomass allocation.The declining trend in the leaves/stem ratio was counteracted by the differential response to the size changes of the leaves and stem respiration.Whereas the specific respiration of the woody tissues (stem and roots) significantly declined with size, the leaves maintained nearly constant specific respiration rates, hence maintaining their relative contribution to shoot respiration.It is noteworthy that only the scaling of root respiration was not significantly larger than the b factor of the 3/4 predicted by the West, Brown, and Enquist (WBE) plant-scaling model [33], but even in this case, its confidence interval was centered at 6/7 (0.73-0.96).This is in agreement, however, with previous empirical evidence in young trees, showing b values lower than 1, but generally higher than 3/4 [30,52].
Although the stem/root respiration rates increased with the tree size, this trend was much weaker than expected from the changes in the biomass ratio.This suggests that, despite the declining biomass, the metabolic allocation to the roots was comparatively larger in the dominant trees.We also found a larger size-independent variability in the root respiration than in any other tissue, potentially linked to changes in the proportion of the fine and coarse roots among the individuals, as the root respiration tends to be negatively associated with the root diameter (see e.g., [53]).Although a decline in the proportion of the fine roots with size has been reported when comparing the young and mature stands of Sugi [50], this trend might be less consistent in a dominance-driven size gradient.In this regard, a further assessment of the biomass and metabolic partition within the root system is needed in order to understand the observed patterns in root respiration.
In the case of the leaves, b did not differ significantly from the unity, thus supporting the assumption that canopy metabolism results from the addition of N individual leaves with invariant metabolic rates [54].Bearing in mind that the leaf respiration is largely substrate-limited, and thus dependent on leaf carbon fixation [55,56], our results could be interpreted as a lack of response of the taller trees to their higher light availability, apparently contradicting the observed trends in the growth rate.However, the leaf mass per area is known to increase with the tree height and canopy openness [40], and thus our results would be compatible with higher photosynthetic rates per leaf area in the larger trees.In this regard, it would be worth expanding the empirical assessment of the tree metabolic scaling by considering the differential role of mass and area expansion.

Transpiration Rates Show a Negative Allometric Response to Tree Size
Contrary to our hypothesis, the transpiration rates in the young competing trees did not increase linearly with the tree size, with the lower leaf-specific transpiration in the larger trees, thus contradicting the assumption of the invariant metabolic rates in the leaves [54].Although the b coefficients in all of the tissues were closer to the linearity than the respiration in the woody tissues, they were less linear than the leaf respiration rates.Besides, we found a positive association between size and WC in the root and stem, indicating improved water distribution and larger water storage capacity in the large trees.However, although the leaves showed a variation in WC, it was not related to the tree size, suggesting that the WC in the leaves should be physiologically maintained within a certain threshold, at the expense of the stem and root water storage.In this regard, Kumagai et al. [57] found significant stem water storage effects, which could explain the different response to the transpiration rates in the WC at the base and the top of the stem.On the other hand, Himeno et al. [42] found that the leaves in young trees show significant changes in WC, and constitute a relevant storage pool for daily variations in water demand.Because of the longer time-scale considered in our study, we could not assess the daily changes in leaf WC, but this could be behind the scattering found in this variable.
The observed scaling factor b was compatible with the predicted 3/4 scaling, according to the WBE model [33], but the fitted values were larger than in previous empirical studies (for the shoots, b = 0.81 in our study, and b = 0.74 in the literature [58]).Unlike in our study, Meinzer et al. [58] included large, adult trees in their models, which, in agreement with the trends observed in our study, would exacerbate the limitation of transpiration with the increasing size.Our results indicate that the scaling factor for transpiration in young trees is closer to linearity than in the adult trees, similar to what has been reported for respiration [30,52].According to our findings, the leaf mass was the most critical factor determining the size-scaling of transpiration; therefore, the progressive reduction in the leaf/stem ratio could contribute to the negative allometry in the association between the aboveground biomass and transpiration, particularly in the larger trees.This trend would be reinforced by the significant decline in water use per leaf mass with the tree size, although, as for respiration, this trend could be mainly driven by an increase in the leaf mass per area [40].In this regard, despite methodological differences, our scaling models successfully predicted the tree water use in similarly young trees without competition (Himeno et al. 2017 [42]), but substantially overestimated the water use in much larger trees (Kumagai et al. 2009 [57]) (See Appendix B, Figure A2).
At this point, it should be noted that the larger trees tend to receive more direct sunlight, and hence may have a larger integrated evaporative demand [40].Thus, the empirically determined scaling factors include both the physiological scaling and canopy-height environmental gradients.In our study, although we found a positive association between the average daily maximum PPFD at 2/3 of the tree height and WU (y = 0.0005336x 2.155 , R 2 = 0.76, p < 0.001), it was much weaker than the associations found between the fresh weight and WU (R 2 = 0.91 − 0.95, p < 0.001).Furthermore, after standardizing the WU by the PPFD values (WU/PPFD), we still found a strong positive association with the fresh weight, particularly strong for the leaves (R 2 = 0.77 − 0.84, p < 0.001; see Appendix B, Table A1).Overall, this suggests that the observed association with the PPFD was mainly an indirect effect of the correlation between the PPFD and tree size.In any case, although removing the effect of the canopy gradients may be useful for the empirical validation of the theoretical models, they cannot be obviated when the aim is to attain a realistic estimate of the size-scaling under natural conditions.

Hydraulic Limitations in Suppressed Trees May Reduce Their Water-Use Efficiency
Considering that dominant trees show a larger biomass production (BAI), and constant (i.e., not declining) mass-specific leaf metabolic rates, but show a decline in the leaf-specific water use, the dominant trees would be more efficient in their use of water than the suppressed individuals.In our study, we found a significant decline in the stem density with the tree size.The stem density tends to be inversely correlated with the conduit diameter, and hence can be associated with the stem water conductivity [59][60][61].This would explain the negative association between the stem density and WU and Vel, indicating a greater potential for the water transport per unit of sapwood area in the large, dominant trees.The distribution of the variables in the PCA revealed that the negative effect of the stem density on the WU and Vel was largely associated with a direct size effect.Therefore, the smaller trees were also the most limited by stem density, despite the larger demand for water transport in the dominant trees.From these results, we can conclude that the self-pruning was effective in adapting the water demand to the hydraulic capacity of the tree.This interpretation is supported by the steeper increase in the water concentrations with the size in the upper part of the stem, which suggests a greater ability to mobilize water resources along the stem in the larger trees.
Interestingly, we also found a negative association between the stem density and leaf-specific respiration rates, highlighting that the hydraulic limitations affected not only the transpiration but also the leaf metabolic rates.According to the strength of the association with the stem density and its distribution along the PCA axis, the stem density more strongly limited the leaf metabolic activity than the transpiration rates, contributing to the overall reduction in the water-use efficiency for the smaller trees.This might be a general trend, at least in the energy-limited environments, coinciding with the findings by Fardusi et al. [22], who presented a global positive intraspecific association between water-use efficiency and growth, particularly strong in the boreal conifers.

Conclusions
Contrary to our first hypothesis, the belowground to aboveground carbon allocation did not change with tree size and dominance status, but we found significant changes in the carbon allocation among the tree organs.On the one hand, the increase in the stem/root ratio with size, highlights that when the soil water is not limiting, the trees promote height growth and water transport through the stem, at the expense of root development.On the other hand, the leaf canopy development in larger trees was limited, coupling with the root size; hence, modulating the total tree water use to fit the water uptake capacity.
Our results partially supported our second hypothesis, namely, that the size-scaling factor of whole-tree transpiration was larger than the predictions of mass-scaling models, but still below unity.Interestingly, this was partly due to a decrease in the leaf-specific transpiration in larger trees, contrasting with the invariance in the leaf metabolic rates.
Comparing the scaling of water use and carbon metabolism, we can conclude that larger trees tend to increase the water-use efficiency, maintaining leaf metabolic rates while reducing the water use.According to the associations found with stem density, this can be attributed to larger hydraulic limitations in smaller, suppressed trees.Future directions of research should strengthen the focus on the interaction between water use and carbon metabolism at the whole-tree level, and consider the contrasting response of different plant organs in determining whole-plant metabolic scaling.

Figure 2 .
Figure 2. Deuterium-labelling method.(A,B) Injection of the D 2 O tracer in the trunk; (C) D 2 O mixing with the pool of water transpired; (D) sampling of the transpired water for isotope analyses; (E) scheme of the D 2 O-labelling peak and the main curve parameters that can be obtained.X-axis-time; Y-axis-D 2 O concentration above reference values in transpired water; Inj.-time of label injection; T1-time to reach 10% of the label; A-peak area; D2-exponential decay of the label.

Figure 3 .
Figure 3. (a) Fresh weight of the different plant organs for each individual, and in relation to total fresh weight; (b) Log-log plot of the association between total fresh weight and the ratio between different plant organs (on a fresh-weight basis).

Figure 4 .
Figure 4. Log-log plot of the relationship between fresh weight and respiration rates (a) for the root and the shoot; and (b) for the stem and leaves.

Figure 5 .
Figure 5. Log-log plot of the relationship between fresh weight and mass-specific respiration rates, on a fresh weight basis (a) for the root and the shoot; and (b) for the stem and leaves.

Figure 6 .
Figure 6.Relationship between log-transformed fresh weight of the different tree parts and (a) time-integrated whole-plant water use; and (b) tracer velocity.

Figure 7 .
Figure 7. Relationship between log-transformed total fresh weight and (a) water concentration (% of fresh weight) in the different tissues; and (b) stem density (dry weight per unit of stem volume).

Figure 8 .
Figure 8. Bi-plot of the principal component analysis including metabolic, hydraulic, and morphometric traits.RT, Vel, tracer residence time, and velocity; SD, stem density; L/S, S/R, leaf to stem and stem to root fresh weight ratio; R-sR, S-spR, L-spR, root, stem and leaf specific respiration; tree-FW, tree-R, tree-WU, total fresh weight, respiration and time-integrated water use; L-sWU, leaf-specific water use; S-WC, Stem water concentration.Symbol sizes are proportional to tree height.

Table 1 .
Mean and standard error, range, and median for tree morphometric variables (N = 25).
1 BAI-3y-basal area increment in the last three years; FW-fresh weight.

Table 2 .
Mean and standard error, range, and median for water use and storage dynamics (N = 25).
1Leaf-specific WU-WU per unit of leaf fresh weight.

Table 3 .
Summary statistics and coefficient estimates for the association between the most relevant hydraulic traits and organ respiration (N = 25).